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Summary 

The  research  has  identified,  characterized,  and  quantified  various  important  domains  of 
behavior  in  the  nonlinear  distortion  and  disintegration  of  injected  liquid  fuel  streams.  Var¬ 
ious  liquid-stream  configurations  resulting  from  fuel  injectors  have  been  analyzed:  conical, 
annular,  and  planar  streams  with  and  without  swirl;  twin-fluid  and  single-fluid  atomizers. 
Linear  and  nonlinear  theories  of  distortion  and  disintegration  have  been  developed  and  have 
predicted  initial  stream  break-up  characteristics.  Distinct  regimes  of  ligament  break-up  and 
cellular  break-up  have  been  determined.  Modulations  of  both  liquid  streams  and  gas  streams 
have  been  studied  as  means  of  active  control.  The  characteristics  of  two-dimensional  cap¬ 
illary  wave  phenomena  have  been  determined.  The  effects  of  impacting  gas  jets  have  been 
compared  with  the  Kelvin-Helmholtz  effect  of  parallel  jets.  Pulsed  gas  jets  have  been  shown 
to  be  more  effective  than  gas  jets  as  a  break-up  mechanism.  The  importance  of  the  rate  of 
forced  stretching  of  the  liquid  stream  compared  to  the  rate  of  disturbance  propagation  in  the 
stream  (i.e.,  capillary  wave  velocity  or  characteristic  viscous  velocity)  has  been  quantified; 
various  regimes  for  the  forced  stretching  have  been  identified  and  characterized. 

I.  Introduction 

The  research  program  focused  on  an  extension  of  previous  studies  on  planar  and  axisym- 
metric  swirling  and  non-swirling  annular  liquid  sheets.  Effects  of  a  surrounding  gas-phase  on 
the  distortion  and  disintegration  of  the  injected  liquid  stream  has  beeen  analyzed  for  a  prac¬ 
tical  two-dimensional  twin-fluid  atomizer  configuration.  Directly  (through  a  liquid-phase) 
or  indirectly  (via  a  gas-jet)  modulated  sinuous  (anti-symmetric)  and  dilational  (symmetric) 
film  distortion  have  been  analyzed  and  evaluated  with  respect  to  the  atomization  efficiency. 

Active  control  of  liquid  stream  distintegration  has  been  analyzed  for  three-dimensionally 
modulated  swirling  conical  liquid  films.  The  latter  configuration  is  of  practical  relevance  in 
state-of-the-art  injector  systems  of  combustion  engines.  In  this  context,  a  linear  analysis  of 
three-dimensionally  distorting  swirling  annular  liquid  films  has  been  performed  in  order  to 
provide  guidance  for  the  studies  of  nonlinearly  distorting  swirling  annular  and  conical  sheets 
with  three-dimensional  disturbances.  Comparison  between  linear  and  nonlinear  results  for 
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the  annular  case  quantifies  the  influence  of  nonlinear  effects  on  film  distortion  and  film 
breakup. 

In  an  effort  to  understand  better  the  dynamics  of  the  larger  ligaments  first  detached 
from  a  continuous  discharging  liquid  stream,  an  analysis  of  the  dynamic  stretching  of  a 
two-dimensional  liquid  film  has  been  conducted.  This  analysis  included  the  consideration  of 
liquid  viscosity  and  long  range  intermolecular  forces  in  addition  to  surface  tension  and  liquid 
inertia. 

In  the  following  three  sections,  we  will  summarize  major  findings  from  the  described 
analyses  of:  planar  liquid  films  discharging  from  twin-fluid  atomizers;  the  capillary  stability 
of  modulated  swirling  liquid  films;  and  the  dynamic  stretching  of  a  planar  liquid  bridge.  The 
details  of  the  analyses,  results  and  literature  reviews  are  provided  in  References  1  through 
5;  since  they  are  not  yet  in  print,  copies  of  References  1,4  and  5  are  attached  as  addenda 
to  this  report.  References  6  though  13  are  conference  papers  and  a  book  article,  the  former 
preceeding  the  primary  journal  articles  (1-5). 

II.  Planar  Liquid  Film  Discharging  From  Twin-Fluid  Atomizer 

The  flowfield  in  the  vicinity  of  a  twin-fluid  atomizer  has  been  analyzed  numerically.  See 
Reference  1.  Considered  is  a  thin  two-dimensional  inviscid  incompressible  liquid  film  dis¬ 
charging  from  the  atomizer  centerline  surrounded  on  both  sides  by  gas  jets  that  impact 
symmetrically  (in  phase)  or  antisymmetrically  (out  of  phase)  onto  the  discharging  film.  The 
value  of  the  ambient  pressure  is  manifested  through  the  density  ratio.  Initial  film  distor¬ 
tion  is  enforced  actively  by:  1)  modulation  of  the  impacting  gas  jets  (indirect  or  gas-phase 
modulation)  with  gas  momentum  components  parallel  and  normal  to  the  liquid  stream,  2) 
modulation  in  the  liquid-phase  at  the  atomizer  exit  (direct  or  liquid-phase  modulation)  or 
3)  direct  liquid  film  modulation  with  superposition  of  continuous  gas  streams.  The  inves¬ 
tigation  focused  on  gas-phase  modulated  films.  Effects  of  different  flow  parameters  on  film 
breakup  characteristics  were  studied  by  varying  one  of  the  flow  parameters  from  a  given 
base-case  configuration. 

For  direct  (i.e.,  gas-phase)  modulated  films,  energy  input  at  the  gas-inlet  ports  was  found 
to  be  inadequate  in  identifying  effective  film  rupture  conditions,  the  latter  being  characterized 
by  short  break-up  length  and  times  at  minimum  energy  input  into  the  system.  For  the 
parameter  domain  surrounding  the  considered  base  case,  film  rupture  was  obtained  more 
effectively  by  sinuous  forcing  and  subsequent  nonlinear  sinuous-dilational  mode  coupling 
rather  than  by  dilational  forcing  at  the  same  average  energy  flux  into  the  system.  Also,  it 
was  observed  that  intermediate  gas-jet  pulsing  (without  underlying  continuous  gas  flow)  can 
yield  the  same  or  larger  amplitude  disturbances  than  obtained  through  wave  growth  from 
continuous  jets  at  the  same  overall  energy  input  into  the  system. 

Direct  modulation  of  liquid-phase  flow-parameters,  i.e.  axial  or  transverse  film  velocity, 
at  the  nozzle  exit  is  more  effective  than  gas-jet  modulation,  the  latter  transferring  only  a 
portion  of  the  modulation  energy  onto  the  film.  Computational  results  also  indicate  that 
for  direct  liquid-phase  forcing,  dilational  modulation  (i.e.  modulation  yielding  dilational 
film  distortion)  is  more  effective  in  causing  film  rupture  than  sinuous  modulations,  while, 
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as  noted  earlier,  the  opposite  is  true  for  gas-jet  modulated  liquid  films.  Also,  for  liquid- 
phase  modulated  films,  admission  of  ambient  gas  streams  significantly  reduces  sinuous  mode 
wavenumbers  whereas  influence  on  wavelengths  remains  small  for  dilationally  distorting  films. 

For  gas-phase  modulated  films,  variations  of  Weber  number,  gas-jet-to-liquid-film  mo¬ 
mentum  ratio,  gas-to-liquid  density  ratio  and  pulse  period  were  considered  and  their  effect 
on  film  distortion  was  analyzed.  An  increase  in  Weber  number  from  its  base-case  value 
(i.e.  from  10  to  25)  resulted  in  stretched  band-like  films  and  a  delay  in  film  rupture  due  to 
reduced  sinuous-dilational  mode  coupling.  On  the  other  hand,  reduction  in  Weber  number 
(i.e.  from  10  to  5)  resulted  in  smaller  disturbance  amplitudes  and  larger  break-up  time  and 
length  despite  stronger  nonlinear  mode  coupling  manifested  by  the  contraction  of  the  film 
into  fluid  cylinders  connected  by  thinner  fluid  films.  Transverse  deflection  of  the  film  center- 
line  for  a  five-fold  increase  in  density  ratio  was  similar  to  the  one  observed  for  the  increased 
Weber  number  case;  however,  fluid  blob  formation  was  still  observed. 

A  five-fold  increase  in  vertical  or  parallel  gas-jet-to-liquid-film-momentum  ratio  resulted 
in  “immediate”  film  rupture  caused  by  the  dynamics  of  the  impacting  gas  jets  rather  than  by 
continuous  growth  of  film  disturbances  downstream  with  energy  transfer  from  the  adjacent 
gas-streams. 

Variation  in  pulse  period  affects  film  distortion  in  two  ways:  It  changes  the  amplitude  of 
the  initial  film  disturbances  generated  by  the  pulsed  gas-jets,  since  increasing  the  pulse  period 
provides  more  time  for  transverse  film  movement  (due  to  one-sided  gas-jet  impact)  before 
being  counter-acted  by  an  opposed  gas  jet.  On  the  other  hand,  it  alters  the  wavelengths 
of  the  disturbances  generated  on  the  film  and  therefore  the  gas-to-liquid  energy  transfer 
downstream  from  Kelvin-Helmholtz-type  wave  growth. 

At  constant  gas-to-liquid  jet  momentum  ratio,  higher  ambient  pressure  conditions,  i.e.  higher 
gas-to-liquid  density  ratios,  will  affect  film  distortion  and  disintegration  only  through  changes 
in  the  film  dynamics  or  interface  dynamics.  If  the  ambient  pressure  or  gas-to-liquid  density 
ratio  is  increased  while  keeping  the  gas-  and  liquid-phase  velocities  constant,  film  distortion 
is  determined  by  the  interface  dynamics  directly  but  also  indirectly  through  the  changed 
gas-phase  momentum  influx  which  will  affect  the  force  balance  at  the  interfaces.  With  in¬ 
creasing  density  ratio,  larger  sections  of  the  discharging  liquid  film  are  being  deflected  by 
the  impacting  gas  jets  and  film  rupture  time  is  greatly  reduced.  Rupture  or  rather  film 
tearing  occurs  before  the  development  of  a  wave  structure  on  the  discharging  film.  A  similar 
behavior  is  observed  if  gas  and  liquid  injection  velocities  are  kept  constant  rather  than  the 
momentum  ratio.  The  density-ratio  effects  described  above  are  identical  in  the  limit  where 
the  gas-jet-momentum  influx  is  reduced  to  zero.  In  that  case,  i.e.  for  a  liquid  film  moving 
through  a  quiescent  gas,  linear  (and  also  nonlinear)  theory  shows  that  (temporal)  growth 
rates  of  unstable  sinuous  (and  dilational)  waves  present  on  the  liquid  film  increase  with 
increasing  density  ratio  at  fixed  liquid  Weber  number. 

Details  of  the  described  analysis  and  of  the  results  are  provided  in  the  first  addendum 
(Reference  1)  to  this  report. 
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III.  Capillary  Stability  of  Modulated  Swirling  Liquid  Films 

Linear  and  nonlinear  analyses  of  modulated  three-dimensionally  distorting  thin  inviscid 
free  liquid  films  discharging  into  a  gas  of  negligible  density  have  been  presented.  See  Ref¬ 
erences  2,  3  and  4.  The  nonlinear  numerical  analysis  uses  a  lubrication  model  reducing  the 
three-dimensional  problem  to  a  system  of  two-dimensional  unsteady  equations.  Linear  the¬ 
ory  for  swirling  annular  films  predicts  that  for  We  >  2  only  one  unstable  wave  is  generated 
on  the  film  due  to  its  modulation  at  the  nozzle  exit.  Depending  on  Weber  number  and  film 
radius,  linear  growth  rates  for  this  unstable  dilational  mode  wave  might  be  larger  for  non- 
axisymmetric  modes  than  for  the  corresponding  axisymmetric  case.  However,  for  the  consid¬ 
ered  parameter  range,  maximum  growth  rates  for  a  given  Weber  number  and  annular  radius 
are  still  observed  for  the  axisymmetric  case.  Modulation  of  multiple  dilational  mode  waves 
at  the  nozzle  exit  allows  for  uniform  breakdown  of  swirling  and  non-swirling  annular  films 
indicated  by  the  generation  of  uniformly  sized  liquid  volumes  which  are  uniformly  spaced 
and  connected  by  thinner  liquid  layers.  Film  topology  and  break-up  of  three-dimensionally 
modulated  conical  films  were  discussed.  Comparison  with  the  corresponding  swirling  annular 
film  shows  that  for  sinuous  film  modulation,  film  divergence  causes  the  separation  of  initially 
formed  larger  fluid  blobs  into  a  pair  of  smaller  fluid  volumes  connected  by  an  even  thinner 
liquid  layer. 

On  clockwise  swirling  conical  films,  spiraling  dilational  and  sinuous  waves  moving  in  the 
clock-wise  direction  increase  in  slope  as  the  film  thins  out  in  the  downstream  direction.  The 
slope  of  counter-clockwise  propagating  waves  decreases  with  downstream  distance.  The  de¬ 
scribed  changes  in  slope  can  be  attributed  to  a  decrease  in  swirl  velocity  with  increase  in 
annular  film  radius  downstream  that  results  from  the  conservation-of-angular-momentum 
principle.  Based  on  the  location  of  the  initial  film  rupture  points  and  the  thickness  dis¬ 
tribution  at  the  time  of  film  rupture,  characteristic  break-up  pattern  have  been  identified 
for  the  cases  with  pure  standing  or  travelling  dilational  or  sinuous  mode  waves  modulated 
at  the  nozzle.  For  mixed  standing/travelling  and  sinuous/dilational  wave  modulation,  the 
break-up  pattern  are  modified  or  mixed  versions  of  the  patterns  identified  for  the  pure  cases. 
For  the  investigated  cases  with  travelling  dilational  or  sinuous  mode  circumferential  waves, 
superposition  with  the  imposed  axial  wave  results  in  an  oblique  wave  spiraling  clock-wise 
downstream  on  the  annular  or  conical  film.  Here,  initial  film  breakup  occurs  simultane¬ 
ously  at  various  points  with  constant  downstream  distance.  The  perforations  are  expected 
to  expand  along  the  line  of  minimum  film  thickness  which  results  in  the  formation  of  liquid 
filaments  spiraling  downstream  in  the  counter-clockwise  direction  (“filament  break-up”).  On 
swirling  annular  films  with  superimposed  dilational  standing-wave  modulations  in  the  cir¬ 
cumferential  direction,  film  rupture  first  occurs  after  the  formation  of  larger  fluid  blobs  and 
upstream  just  behind  these  larger  fluid  masses.  Again,  initial  film  rupture  will  take  place 
simultaneously  at  various  locations  at  the  same  downstream  position.  Subsequent  expansion 
of  the  film  perforations  can  be  expected  to  generate  a  more  or  less  circular  pattern  within 
each  cell  formed  by  neighboring  fluid  blobs.  This  type  of  film  rupture  can  be  characterized 
as  “cellular  break-up”.  The  same  break-up  pattern  can  also  be  observed  for  non-swirling 
(pressure-stabilized)  annular  films  and  for  swirling  conical  films  with  similar  forcing  con- 
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ditions.  However,  due  to  film  divergence  the  “cells”  will  stretch  in  the  conical  case  and 
film  thickness  within  the  cells  will  be  reduced.  Therefore,  expansion  of  the  film  perforation 
is  expected  to  occur  significantly  faster  than  in  the  corresponding  annular  case.  However, 
initial  rupture  of  the  film  in  the  conical  case  is  delayed.  Analogous  to  the  dilational  case 
discussed  earlier,  initial  film  perforation  for  swirling  annular  films  with  modulated  standing 
sinuous  circumferential  waves  takes  place  just  behind  the  larger  fluid  blobs.  However,  in  this 
case  the  perforations  are  expected  to  propagate  in  the  circumferential  direction  at  a  (more 
or  less)  constant  downstream  distance,  ultimately  resulting  in  the  detachment  of  fluid  rings 
with  pronounced  thickness  fluctuations  in  the  circumferential  direction  (“ring  break-up”). 
The  similar  standing  wave  sinuous  mode  modulation  imposed  onto  the  conical  film  geometry 
yields  a  “filament  break-up”  pattern.  The  filaments  are  spiraling  clock- wise  in  the  down¬ 
stream  direction.  In  contrast  to  the  single  wave  modulation,  breakup  of  these  filaments  will 
be  greatly  influenced  by  the  existing  non-uniform  mass  distribution  along  the  filaments. 

Details  of  the  three-dimensional  analysis  and  the  preceeding  axisymmetric  analysis  of 
swirling  annular  and  conical  films  are  provided  in  References  2,  3  and  4.  Reference  4  is  not 
yet  in  print  so  a  copy  is  provided  as  second  addendum  to  this  report. 

IV.  Dynamic  Stretching  of  A  Planar  Liquid  Bridge 

A  thin  incompressible  viscous  planar  free  liquid  film  in  a  void  and  under  zero  gravity  was 
analyzed  by  means  of  the  previously  developed  reduced-dimension  (lubrication)  approach 
extended  to  include  liquid  viscosity  and  long-range  intermolecular  forces.  Linear  analysis 
focused  on  films  with  harmonic  modulations  in  the  axial  film  velocity  enforced  at  the  ends 
of  the  planar  bridge.  Effect  of  changes  in  the  problem  parameters  on  the  overall  distortion 
characteristics  of  the  film  were  discussed.  Nonlinear  film  distortion  and  break-up  was  in¬ 
vestigated  for  the  case  of  temporally  increasing  velocity  at  the  end  of  the  film  resulting  in 
continuous  film  stretching  eventually  leading  to  film  rupture.  Implementation  of  the  em¬ 
ployed  numerical  model  was  validated  for  the  linear  limit  by  comparison  with  the  analytical 
linear  solutions  and  for  harmonically  modulated  film-end  velocities. 

The  importance  of  the  film  stretching  rate  in  comparison  to  either  the  capillary  wave 
velocity  or  the  characteristic  viscous  velocity  has  been  demonstrated.  If  the  stretching  rate 
is  large,  the  distortion  signals  are  slow  in  reaching  large  distances  from  the  forced  end; 
so  distortion  is  confined  to  a  smaller  portion  of  the  sheet.  For  small  stretching  rates,  the 
distortion  is  significant  over  a  larger  portion  of  the  film. 

For  films  with  an  initial  length-to-thickness  ratio  of  0(10)  and  with  a  time-scale  for 
film-end  acceleration  (from  zero  to  maximum  pull  velocity)  comparable  to  the  propagation 
time  of  capillary  waves  from  film  end  through  film  center,  various  distinct  film  topologies 
were  observed,  depending  on  Weber  number  and  Reynolds  number.  Here,  film  topology  is 
typically  characterized  by  three  distinct  regions,  i.e.  a  film  wedge  forming  at  the  pulling 
end(s),  the  film  center  region  and  a  transition  region.  The  size  and  shape  of  these  regions 
greatly  depends  on  the  particular  case  under  investigation. 

For  cases  with  length-to-thickness  ratio  of  0(100)  information  on  the  acceleration  of  the 
film  ends  does  not  reach  the  film  center  region.  Here,  evolution  of  film  topology  is  similar 
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for  all  the  considered  Weber  number  and  Reynolds  number  combinations;  with  a  more  or 
less  undisturbed  film  center  region  smoothly  transitioning  into  a  narrow  liquid  wedge  at  the 
pulling  end.  For  cases  with  low  film-end  acceleration,  film  rupture  is  significantly  delayed 
resulting  (in  all  the  considered  cases)  in  a  significantly  stretched  film  prior  to  rupture  with 
or  without  significant  amounts  of  fluid  remaining  in  the  film  center  region. 

Nonlinear  film  distortion  has  also  been  investigated  for  the  continuously  compressed 
planar  film  bridge,  illustrating  the  relevance  of  the  film  bridge  analysis  for  the  contracting  and 
stretching  of  a  free  planar  film  such  as  those  found  in  the  atomization  process  of  liquid  fuels 
in  typical  gas  turbine  combustors.  Film  distortion  characteristics  observed  for  continuously 
compressed  planar  films  conform  with  observations  made  by  other  authors  for  the  similar 
case  of  contracting  free  liquid  films. 

Details  of  the  described  stretching-film  analysis  are  provided  in  Reference  5;  since  it  is 
not  yet  in  print,  a  copy  is  attached  as  the  third  addendum  to  this  report. 
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Abstract 

The  distortion  and  break-up  of  a  thin  planar  liquid  film  impacted  by  two 
gas  jets  while  discharging  from  a  twin-fluid  atomizer  is  studied  numerically. 

The  gas  momentum  vector  has  components  normal  and  parallel  to  the  liquid 
stream.  Viscosity  and  compressibility  are  neglected  in  both  the  liquid  phase 
and  the  gas  phase.  The  reduced-dimension  (lubrication)  approximation  is 
employed  to  describe  the  nonlinear  distortion  and  breakup  of  the  thin  film. 

The  gas-phase  dynamics  are  modelled  by  using  a  boundary-element-method 
formulation.  For  the  considered  parameter  range  and  for  a  given  energy  ex¬ 
penditure,  direct  modulation  of  liquid-phase  velocities  at  the  nozzle  exit  is 
found  to  be  more  effective  in  causing  film  rupture  than  indirect  modulation 
via  adjacent  impacting  gas  jets.  In  the  former  case,  dilational  film  modu¬ 
lation  results  in  shorter  breakup  lengths  than  sinuous  modulation.  On  the 
other  hand,  for  gas-jet  modulated  films,  sinuous  mode  forcing  is  more  effec¬ 
tive  than  dilational  forcing  for  the  same  energy  input.  Co-flowing  gas  streams 
significantly  alter  wavelengths  and  amplitudes  of  film  disturbances  generated 
by  direct  film  modulation.  Large  ratios  of  gas-jet  momentum  to  liquid-film 

i 

momentum  result  in  “immediate”  film  rupture  in  response  to  the  dynamics  i 


of  the  impacting  gas  jets,  whereas  for  lower  ratios  films  disintegration  occurs 
further  downstream  after  continuous  growth  of  the  initial  disturbances.  Film 
distortion  is  characterized  by  the  formation  of  fluid  blobs  or  long  band-like 
films  depending  on  Weber  number  values  and  density  ratio. 


Typeset  using  REVT^K 


I.  INTRODUCTION 

In  the  past,  free  liquid,  films  have  not  been  analyzed  only  for  their  intriguing  scientific 
substance  but  also  due  to  their  relevance  for  a  wide  variety  of  technological  applications, 
including  the  formation  of  liquid  sprays  from  fluid  films  which  discharge  from  so-called  at¬ 
omizer  nozzles.1  The  manner  in  which  these  films  disintegrate  into  droplets  depends  upon 
the  type  and  geometry  of  the  atomizer  as  well  as  on  the  operating  conditions.  One  com¬ 
mon  type  of  atomizer  is  described  as  the  twin-fluid  atomizer,  where  the  interfaces  of  the 
discharging  liquid  film  are  subjected  to  neighboring  gas  streams.  Here,  the  main  cause  of 
film  instability  is  interaction  of  the  film  with  the  surrounding  gas  flow  either  through  gas-jet 
impact  or  continuous  growth  (amplification)  of  film  disturbances  via  energy  input  from  the 
neighboring  gas  streams. 

Subsequently,  we  review  previous  work  on  discharging  free  liquid  films  in  an  ambient 
gas  phase;  some  references  are  also  provided  considering  theoretical  or  numerical  analyses 
of  periodically  disturbed  infinite  films  (also  with  the  effect  of  an  ambient  gas  considered). 
In  subsection  A,  we  review  analyses  dealing  with  films  under  steady  discharge  and  ambient 
conditions  (or  boundary  conditions).  Subsection  B  addresses  previous  work  on  discharging 
liquid  films  with  modulations  enforced  onto  the  liquid  phase.  In  subsection  C,  we  describe 
the  flow  configuration  considered  within  the  present  work. 

A.  Steady-State  Discharging  Films 

Until  the  recent  past,  steady-state  operational  conditions  were  imposed  onto  the  pre¬ 
scribed  atomization  systems.  Accordingly,  the  majority  of  previous  experimental  and  an¬ 
alytical  work  focused  on  the  aerodynamic  stability  of  these  films  under  steady  free-stream 
and  boundary  (e.g.,  inflow  and  outflow)  conditions. 

Arai  &  Hashimoto2,3  and  Hashimoto  &  Suzuki4  presented  experimental  results  on  the 
three-dimensional  breakdown  of  a  planar  liquid  sheet  in  a  high-speed  co-current  uniform 
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gas  stream  under  steady-state  conditions;  more  extensive  analyses  of  a  similar  configuration 
were  later  conducted  by  Mansour  &  Chigier,5,6  by  Stapper  et  al.7,a  and  discussed  by  Feman- . 
des  et  al.9.  In  Ref.  9,  the  authors  illustrate  that  effectiveness  of  atomization  depends  not 
only  on  gas-to-liquid  momentum  ratio  but  also  on  film  thickness,  which  is  attributed  to  the 
observation  that,  for  effective  film  disintegration,  the  characteristic  length  scale  of  pertur¬ 
bations  should  be  of  the  order  of  the  film  thickness.  The  authors  also  observe  that,  for  high 
liquid  momentum,  the  effect  of  viscous  shear  on  film  breakup  is  confined  to  a  thin  region 
near  the  gas-liquid  interface.  In  this  case,  atomization  can  be  improved  by  superposition 
of  artificial  film  modulations  which  can  cause  disturbances  of  the  inner  core.  Even  though 
external  liquid  and  gas-phase  modulation  via  piezoelectric  transducers  or  loudspeakers  is 
described  in  Ref.  9,  no  results  have  been  presented  under  such  forcing  conditions. 

Within  their  work,  Aral  &  Hashimoto3  also  presented  a  perturbation  analysis  of  the 
corresponding  inviscid  three-dimensional  boundary-value  problem.  A  linear  analysis  based 
on  the  Orr-Sommerfeld  equation  was  also  attempted  by  Hashimoto  &  Suzuki.4  Ibrahim  &; 
Akpan10  presented  a  three-dimensional  linear  analysis  of  a  periodically  disturbed  planar 
viscous  liquid  sheet  moving  in  an  inviscid  gas  medium  and  with  axial  and  transverse  sheet 
disturbances. 

Kawano  et  al.11  investigated  annular  liquid  films  discharging  from  a  two-fluid  atomizer 
under  various  steady-state  conditions.  The  authors  also  provided  a  linear  temporal  analysis 
for  inviscid  incompressible  annular  films  subject  to  unequal  gas  velocities  on  both  sides  of 
the  annular  film.  A  similar,  experimental  analysis  was  reported  later  by  Berthoumieu  et 
al ,12’13 

Few  nonlinear  analyses  of  planar  sheets  with  surrounding  gas  streams  have  been  reported. 
Rangel  and  Sirignano14  and  later  Lozano  et  a/.17,18  employed  discrete-vortex  methods  to 
describe  the  nonlinear  distortion  and  rupture  of  two-  or  three-dimensionally  distorting  liquid 
films  in  a  co-flowing  gas  stream. 

The  nonlinear  calculations  conducted  in  Ref.  14  revealed  that:  1)  oscillating  (stable), 
sinuous  modes  exist  at  gas-to-liquid  density  ratios  of  order  1,  and  2)  sinuous  distortion  may 
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result  in  ligaments  interspaced  by  half  of  a  wavelength,  whereas  dilational  distortion  may 
result  in  ligaments  interspaced  by  one  wavelength.  These  break-up  features  were  already 
described  by  Dombrowski  &  Hooper15  and  identified  later  by  Jazayeri  &  Li16. 

In  Ref.  17,  aerodynamic  growth  of  antisymmetric  (sinusoidal)  streamwise  sheet  distur¬ 
bances  with  superimposed  symmetric  (dilational)  or  antisymmetric  disturbances  in  the  trans¬ 
verse  direction  was  studied,  in  order  to  better  understand  the  presence  of  transverse  and 
longitudinal  filaments  observed  in  liquid  sheet  air-assisted  atomization  experiments.  For 
both  initially  sinusoidal  or  dilational  transverse  sheet  disturbances,  flattening  of  the  sheet  in 
the  transverse  direction  is  reported.  The  formation  of  streamwise  tubular  sections  with  point 
edges  is  predicted  when  both  air/liquid  interfaces  touch  each  other  at  points*  interspaced  by 
one  wavelength  in  a  plane  perpendicular  to  the  main  flow  (i.e.  streamwise)  direction.  As 
illustrated,  and  due  to  the  fact  that  gas  and  liquid  velocities  are  assumed  to  be  initially 
oriented  in  the  streamwise  direction,  the  transverse  perturbation  grows  at  a  slower  rate  than 
the  streamwise  one  whose  initial  growth  is  exponential  (by  linear  theory).  Accordingly,  the 
transverse  perturbation  growth  rate  might  be  insufficient  to  generate  streamwise  filaments. 
In  this  context,  oblique  waves  have  been  studied18  and  were  found  to  contribute  to  the  evolu¬ 
tion  of  a  finite  amplitude  transverse  wave  starting  from  infinitesimal  perturbations.  Surface 
tension  effects  and  edge  effects  were  addressed  very  briefly  in  Ref.  18. 

B.  Modulated  Discharging  Films 

More  recently,  atomization  systems  are  being  developed  that  aim  at  an  active  control 
of  the  film  disintegration  or  overall  spray  formation  process  in  order  to  optimize  certain 
parameters  within  in  the  overall  system  (e.g.,  NOx  reduction  in  spray  combustion  of  fossil 
fuels19).  Several  experimental  and  theoretical  analyses  have  addressed  this  interest. 

The  spatial  stability  of  thin  liquid  sheets  subject  to  forced  vibrations  applied  to  the  nozzle 
was  first  investigated  by  Hagerty  &  Shea21  and  later  by  Crapper,  Dombrowski  &  Pyott.20  The 
latter  authors  used  the  dispersion  relation  derived  by  Squire,22  in  order  to  analyze  the  spatial 
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(aerodynamic)  growth  of  sinuous  disturbances.  Comparison  of  the  theoretical  results  with 
experimental  observations  revealed  considerable  discrepancies:  1)  unstable  waves  were  found 
at  frequencies  where  linear  inviscid  theory  (based  on  spontaneous  growth  of  infinitesimally 
small  perturbations)  predicts  stable  waves,  and  2)  the  wave  growth  saturates  (i.e.,  the  wave 
amplitude  eventually  no  longer  grows)  as  the  waves  propagate  away  from  the  orifice,  even 
for  wavelengths  which  are  unstable  by  the  employed  linear  theory.  The  latter  phenomenon 
was  attributed  to  the  formation  of  vortices  cast  off  from  the  wave  crests  and  moving  into 
the  troughs. 

After  Crapper  et  al.,20  Asare,  Takahashi  &  Hoffman23  also  investigated  the  stability  of 
liquid  sheets,  harmonically  forced  at  the  nozzle  (sinuous  disturbances  only).  Their  experi¬ 
mental  results  (for  thin  and  thick  sheets  and  variable  air  pressures)  agreed  well  with  linear 
theory21’22  if  the  disturbance  amplitudes  were  small  and  if  the  observations  were  made  near 
the  nozzle  exit.  As  a  result  of  nonlinear  effects,  deviations  between  experimental  observa¬ 
tions  and  linear  theory  were  more  pronounced  at  larger  downstream  positions  and  for  higher 
forcing  amplitudes.  Similar  to  the  observations  made  by  Crapper  et  al.,20  the  wave  envelope 
amplitude  appeared  to  saturate  at  large  falling  distances,  i.e.  where  the  sheet  jets  begin 
to  break  up.  Asare  et  al.2Z  also  provided  a  simplified  trajectory  theory  based  on  the  force 
balance  for  a  fluid  particle  in  the  crest  of  a  propagating  wave.  The  initial  wave  envelope 
amplitudes  predicted  by  this  simplified  theory  agreed  well  with  experimental  results  and  the 
employed  linear  theory. 

The  questions  on  the  discrepancies  between  linear  theory  and  experimental  results  raised 
in  Refs.  20  and  23  were  subsequently  addressed  in  Refs.  24  and  25.  Crapper,  Dombrowski 
&  Jepson24  presented  a  linear  analysis  for  both  sinuous  and  dilational  waves  on  thin  viscous 
liquid  sheets  in  a  viscous  gas  flow  (with  focus  on  sinusoidal  waves  with  wavelengths  much 
longer  than  the  sheet-thickness) .  They  found  that  liquid  viscosity  has  no  effect  on  the  initial 
wave  growth.  Furthermore,  gas  viscosity  was  found  to  extend  the  region  of  instability  over 
a  greater  range  of  frequencies  than  that  given  by  the  inviscid  theory;  which,  as  stated  in 
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Ref.  24,  might  explain  the  discrepancy  with  respect  to  this  issue  found  in  Ref.  20.  (Wrap¬ 
per  et  al.24  also  observed  that  by  their  linear  viscous  theory  (i.e.  viscous  gas  and  liquid 
phases),  a  maximum  growth  rate  does  not  always  exist,  which  (as  stated)  indicates  that,  for 
these  situations, -the  experimentally  observed  waves  cannot  be  the  fastest  growing  as  hith¬ 
erto  suggested  (for  example  by  Squire22)  but  must  be  imposed  on  the  sheet  by  some  other 
means,  e.g.  by  a  natural  frequency  of  the  apparatus.  The  importance  of  nozzle  vibrations  in 
practical  applications  was  subsequently  demonstrated  by  Crapper  &;  Dombrowski.25  Capil¬ 
lary  waves  on  planar  liquid  sheets  emanating  from  a  nozzle  or  atomizer  were  also  observed 
and  investigated  experimentally  by  Hashimoto  &  Suzuki4  and  held  responsible  by  Clark 
&  Dombrowski2®  for  the  prescribed  discrepancy  between  their  theoretical  (infinite-sheet) 
analysis  and  experimental  results  for  fan  sheets. 

More  recently,  an  experimental  study  on  modulated  discharging  films  was  also  presented 
by  Chung  et  al.27  The  authors  provided  an  extensive  study  on  disintegrating  conical  liq¬ 
uid  films  subject  to  piezoelectric  transducer  modulations  inside  the  nozzle.  The  effect  of 
transducer  modulation  on  liquid  film  disintegration  was  studied  for  different  fluid  viscosi¬ 
ties,  driving  frequencies  and  input  perturbation  powers.  An  optimum  driving  frequency  was 
identified  at  which  liquid  sheet  breakup  length  is  shortest.  Increased  liquid  viscosity  was 
found  to  render  liquid  modulation  less  effective  delaying  liquid  film  disintegration.  Higher 
input  modulation  power  enhanced  liquid  disintegration.  The  authors  suggest  that  break-up 
length  is  a  function  of  cosh-1  of  input  power. 

The  linear  analysis  of  a  sinuous  and/or  dilationally  distorting  viscous  liquid  film  in  a 
co-flowing  ambient  inviscid  gas  has  been  presented  recently  by  Mitra  et  al.2S  For  large 
Weber  numbers  (and  by  using  Gaster’s  relation29),  the  authors  extend  their  analysis  to 
liquid  films  discharging  from  a  nozzle  or  atomizer.  Film  breakup  lengths  and  interface 
profiles  are  reported  for  various  combinations  of  Reynolds  number,  Weber  number,  gas- 
to-liquid  density  ratio,  gas-to-liquid  velocity  ratios  and  phase  angle  between  sinuous  and 
dilational  mode  disturbances;  always  assuming  the  presence  of  sinuous  and  dilational  mode 
waves  with  maximum  growth  rates.  Computational  results  are  compared  with  experimental 
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observations  made  by  Jazayeri  and  Li,30  indicating  the  inadequateness  of  their  linear  theory 
for  large  amplitude  disturbances,  i.e.  when  nonlinear  effect  become  important.  In  fact, 
film  breakup  as  predicted  by  the  described  linear  analysis  solely  depends  on  the  dilational 
mode  and  is  independent  of  any  sinuous  mode  disturbances.  However,  large  amplitude 
sinuous  mode  disturbances  have  been  shown  in  Ref.  37,  as  well  as  by  Jazayeri  and  Li16  to 
cause  significant  dilational  film  distortion  due  to  nonlinear  sinuous-dilational  mode  coupling. 
Consequently,  any  predictions  for  film  break-up  length  based  on  linear  theory  has  to  be 
considered  with  great  caution. 

Similar  to  the  described  linear  analysis  for  viscous  films,  Jazayeri  and  Li31  also  extended 
their  nonlinear  analysis  of  Ref.  16  on  inviscid  liquid  films.  As  in  Ref.  30,  the  authors 
used  Gaster’s  relation29  to  transform  the  results  of  their  temporal  analysis  into  a  spatial 
analysis  describing  the  nonlinear  distortion  of  a  discharging  liquid  film  into  a  quiescent 
ambient  gas.  Based  on  this  transformation,  the  authors  determine  nonlinear  film  break¬ 
up  lengths  under  various  operating  conditions.  Analogous  to  Ref.  30,  the  results  from 
the  perturbation  expansion  of  Ref.  31  for  the  spatially  developing  film  always  assume  the 
presence  of  only  the  dominant  wavenumber  at  the  nozzle.  The  results  obtained  for  the 
spatially  developing  film  by  using  the  perturbation  analysis  of  Ref.  31  are  only  valid  for  large 
Weber  numbers  and  liquid  films  that  discharge  into  an  overall  unrestricted  quiescent  ambient 
gas.  Practical  twin-fluid  atomizers  do  not  conform  to  this  assumption.  In  context  with  the 
large  Weber  number  assumption,  it  is  worthwhile  to  notice  that  dilational  capillary  waves  on 
thin  films  are  dispersive  in  nature  with  increase  in  wave  velocity  for  decreasing  disturbance 
wavelength.37  Consequently,  for  finite  Weber  number  values,  propagation  velocities  of  shorter 
wavelength  disturbances  will  always  deviate  from  the  film  discharge  velocity.  Also,  results 
for  spatially  developing  films  obtained  from  a  temporal  analysis  via  transformation  using 
Gaster’s  relation29  do  not  consider  the  physical  constraints  imposed  at  the  nozzle  exit  in  a 
practical  situation. 

The  nonlinear  analysis  presented  here  differs  from  the  analysis  presented  in  Ref.  31.  The 
employed  nonlinear  model  is  based  on  the  assumption  of  thin  films  and  is  exact  in  the  limit 
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of  an  infinite  value  for  the  ratio  of  disturbance  wavelength  to  undisturbed  film  thickness.37 
Considered  is  the  initial-and-boundary  value  problem  of  a  semi-infinite  sinuous  and/or  dila- 
tionally  distorting  inviscid  liquid  film  discharging  from  a  practical  twin-fluid  atomizer  (with 
adjacent  gas-stream)  into  a  planar  channel.  Liquid  film  modulation  is  enforced  by  imposing 
forcing  conditions  onto  the  discharging  gas  streams  which  impact  the  liquid  film  from  both 
sides  at  a  specified  angle.  Direct  modulation  of  the  discharging  film  via  modulation  of  the 
liquid-phase  velocity  at  the  nozzle  exit  is  also  considered.  A  more  detailed  description  of 
the  employed  model  is  given  below. 

It  should  be  noted  here,  that  the  active  control  and  break-up  of  cylindrical  liquid  jets  via 
liquid  phase  modulation  has  a  longer  history  than  the  analysis  of  modulated  liquid  films.  The 
first  detailed  experimental  and  (linear)  analytical  studies  of  modulated  discharging  liquid  jets 
were  conducted  by  Pimbley,32  Pimbley  and  Lee33  and  Bogy34  motivated  by  its  application 
in  ink-jet  printing.  More  recently,  Hilbing  and  Heister35  analyzed  modulated  discharging 
nonlinear  liquid  jets  in  context  with  the  likelihood  of  droplet  recombination  downstream 

j 

from  the  initial  jet  pinch-off  point.  Modulation  of  microjets  has  currently  received  attention 
in  context  with  micro-dispensing  systems  for  pharmaceutical  products  and  in  context  with 
the  formation  of  precision  powders  having  very  uniform  particle  size.  For  a  more  complete 
review  of  previous  theoretical  work  on  free  liquid  films  and  jets,  the  reader  is  referred  to 
Ref.  43. 


C.  New  Configuration  under  Investigation 

The  present  study  is  a  continuation  of  previous  work  by  the  authors.36  In  particular, 
the  present  analysis  considers  the  general  sinuous  and/or  dilational  nonlinear  evolution  of  a 
thin  semi-infinite  planar  liquid  film  downstream  from  a  twin-fluid  atomizer  and  under  the 
influence  of  capillary  and  aerodynamic  effects.  The  numerical  model  employs  the  reduced- 
dimension  approach  of  Ref.  37  to  describe  the  thin  planar  sheet  and  a  Boundary-Element 
Method  (BEM)  formulation38  for  the  inviscid  incompressible  gas  streams  impacting  onto  the 
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film  from  both  sides  of  the  injector  centerline. 

The  unsteady  Bernoulli  equation  in  combination  with  the  boundary-element  method 
is  used  to  determine  the  instantaneous  gas  pressure  at  the  liquid-gas  interfaces  which  is 
needed  within  the  lubrication  equations  governing  the  liquid  phase.  The  use  of  a  discrete 
boundary-element  method  for  the  gas  phase  allows  the  consideration  of  practical  applications 
where  liquid  streams  are  injected  into  a  gaseous  flowfield  with  its  own  physical  constraints 
or  boundary  conditions.  The  latter  is  of  particular  importance  with  regard  to  the  compli¬ 
cated  flowfields  within  and/or  around  fuel  injection  elements  or  atomizers  used  for  spray 
combustion  purposes,  e.g.  prefilming  airblast  atomizers. 

The  particular  injector  configuration  considered  in  this  work  is  shown  in  Fig.l. 

A  thin  liquid  film  is  injected  from  a  two-dimensional  slit-nozzle  located  in  the  symmetry 
plane  of  a  two-dimensional  twin-fluid  atomizer.  The  liquid  film  is  impacted  on  both  sides  by 
gas-jets  discharging  from  two  inlet  ports  angled  with  respect  to  the  film  discharge  plane  at 
an  angle  a.  Both  liquid  film  and  gas  streams  are  injected  into  a  two-dimensional  channel. 
The  width  of  the  gas-jet  inlet  ports  is  s,  the  channel  width  is  H  and  the  width  of  the  liquid- 
phase  nozzle  is  h.  For  the  analysis,  the  coordinate  system  is  fixed  to  the  nozzle  exit  with  the 
y— axis  located  within  the  film  discharge  plane  and  the  x— axis  coinciding  with  the  symmetry 
plane  of  the  atomizer.  For  the  spatially  developing  semi-infinite  liquid  film,  capillary  waves 
generated  by  the  impacting  gas  jet  do  not  reach  the  outflow  boundary  of  the  computational 
domain  during  the  computation. 

The  specification  of  liquid-phase  boundary  conditions  at  the  nozzle  exit  follows  previous 
work  by  the  authors  on  discharging  semi-infinite  films  without  gas-phase  effects.37  In  Ref. 
37,  liquid  film  modulation  was  imposed  only  locally  at  the  nozzle  (x  =  0).  The  number  of 
boundary  conditions  specified  at  the  nozzle  was  chosen  according  to  linear  theory  and  by 
employing  the  Sommerfeld  radiation  condition.  The  latter  condition  excluded  any  wavenum¬ 
ber  solution  with  negative  group  velocity  that  resulted  from  the  specific  forcing  frequency. 
Within  the  present  investigation,  the  presence  of  a  non-zero  density  incompressible  ambient 
gas  and/or  impacting  gas  jets  downstream  from  the  nozzle  exit  plane  allows  generally  for 
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capillary  waves  with  negative  group  velocities  effectively  resulting  in  upstream  energy  trans¬ 
port  on  the  film.  However,  within  this  work  liquid-film  injection  velocities  are  assumed  high 
enough  so  that  upstream  energy  transport  and  upstream  wave  propagation  across  the  nozzle 
exit  plane  can  be  neglected,  eliminating  the  need  to  specify  additional  reflective  boundary 
condition  at  this  boundary.  Nevertheless,  it  should  be  recalled  that  the  incompressibility 
assumption  within  the  gas-phase  implies  that  the  speed  of  sound  is  much  larger  than  the 
the  wave  speed  of  any  capillary  waves  observed  on  the  liquid.  So  signals  may  propagate 
upstream  (and  downstream)  through  the  gas. 

.  Certain  boundary  conditions  have  been  imposed  within  the  gas  phase.  Gas-phase  velocity 
components  normal  to  solid-walls  have  been  set  to  zero  (i.e.  along  the  walls  of  the  discharge 
channel  and  the  atomizer).  This  implies  that,  in  the  boundary-element  analysis,  the  gradient 
of  the  velocity  potential  normal  to  these  boundaries  is  zero.  Parallel  outflow  conditions  have 
been  specified  within  the  gas-phase.  In  particular,  the  value  of  the  velocity  potential  has 
been  prescribed  as  <j>  —  0  at  the  nozzle  exit  on  both  sides  of  the  liquid  film. 

Gas-jet  injection  at  the  two  gas-jet  inlet  ports  has  been  specified,  in  general,  according 
to 

oo  t  i 

%w(*)  =  Y  % o  sin2  [tt(—  -  (2 n  +  nfl;ij2))]  for  (2 n  +  ne. ij2)  <  —  <  (2n  +  1  +  ne. x,2) 

n=0  1P  MP 

and  ua;x,2(t)  =  0  ,  otherwise.  (1) 

Here  Tp  denotes  the  pulse  period  and  subscripts  1  and  2  refer  to  the  lower  and  upper  gas 
jets,  respectively.  The  prescribed  forcing  function  indicates  that  a  gas-jet  pulse  of  period 
Tp  at  either  inlet  ports  is  followed  by  an  equal-time  period  where  there  is  no  pulse  for  that 
particular  port.  Parameters  ng. ij2  describe  any  time  delay  or  phase-shift  in  the  pulsing  of 
the  gas  jets  at  the  lower  (1)  and  upper  (2)  injection  ports.  In  this  analysis  ng; i  =  0  always, 
whereas  ng.2  ranged  from  0  to  1.  Note  that  gas  jet  pulsation  with  ng ^  =  ng$  =  0  will  always 
result  in  dilational  film  deformation  whereas  ngt i  =  0  and  ngt2  =  1  will  generate  initially 
(predominantly)  sinuous  film  distortions. 
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II.  GOVERNING  EQUATIONS 


With  the  prescribed  boundary  conditions,  the  solution  to  Laplace’s  equation  within  the 
gas-phase  is  determined  if  the  instantaneous  location  of  the  liquid-gas  interface  and  the 
velocity  component  normal  to  this  interface  is  known.  More  precisely,  using  the  discrete 
boundary-element  method,  the  solution  for  the  gas-phase  velocity  potential  at  a  given 
point  within  or  on  the  boundaries  of  the  gas-phase  domains  in  Fig.  1  is  determined  by  defining 
contributions  of  sources  from  all  nodes  along  the  boundary  of  each  domain.  Note  that  for 
dilational,  i.e.  symmetric,  film  disturbances  a  boundary-element  solution  has  to  be  obtained 
only  for  gas-phase  region  on  one  side  of  the  liquid  film.  The  flowfield  on  the  opposite  side 
of  the  film  is  then  merely  a  mirror  image  of  the  calculated  one.  However,  for  more  general 
film  disturbances,  boundary-element  solutions  have  to  be  computed  simultaneously  for  both 
gas-phase  regions,  i.e.  above  and  below  the  thin  film. 

The  BEM  approach  adopted  here  to  solve  the  gas-phase  flow  field  follows  the  method 
proposed  by  Brebbia.38 

The  integral  representation  of  Laplace’s  equation  for  the  gas-phase  velocity  potential 
4>g.  1,2  may  be  written  as 


+  JT 


dr  =  o 


where  is  the  potential  at  a  point  fj  in  the  gas-phase  region  below  (1)  or  above  (2) 

the  liquid  film,  r1)2  denotes  the  boundary  of  the  particular  domain,  q\^  =  d4>g^^/dn  is  the 
gradient  of  <f>g- i(2  on  normal  to  the  particular  boundary,  a  is  a  constant  (for  a  given 
node),  and  G  is  the  free-space  Green’s  function  corresponding  to  the  governing  equation. 
Since  the  above  equation  involves  an  integration  only  around  the  boundary  Ti  and/or  T2, 
we  need  not  discretize  the  entire  gas-phase  domain.  It  is  presumed  that  either  <f>g-i,2  or  <?i)2 
is  specified  at  each  ‘node’  on  each  boundary  while  the  other  quantity  is  returned  as  part  of 
the  solution.  Details  regarding  the  BEM  solution  procedure  can  be  found  in  the  appendix  of 
Ref.  39,  together  with  validation  simulations  of  the  employed  method  for  two  steady-state 
problems. 
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Coupling  between  the  gas  and  liquid  phase  is  described  through  kinematic  and  dynamic 
boundary  conditions  at  the  phase  interface.  Kinematic  conditions  are  represented  by  the 
instantaneous  interface  locations  and  the  normal  velocity  of  the  interfaces;  this  affects  the 
solution  Tor- the  velocity  potential  in  the  gas  phase  above- and  below  the  thin  film.  On  the 
other  hand,  the  flowfield  in  the  gas  phase  impacts  the  liquid  film  distortion  through  the 
dynamic  condition  at  the  interfaces,  requiring  that  the  pressure  inside  the  liquid  balances 
with  the  pressure  on  the  gas  side  of  the  interface  combined  with  the  capillary  pressure. 

The  equations  governing  the  nonlinear  film  distortion  are  Eqs.  (2.10)  through  (2.16)  of 
Ref.  37,  which  can  be  combined  to  yield  jg? 

l+^w)=°  <2> 


du  du  =  1  d(pg,+  +Pg-) 

dt  dx  2  p  dx 


+Z.I±  U  +  f-IPy  _  #y 

2fi\dx  2  dr?  u+  ‘zdx* 


fl?ll 

dx2\  j 


fOL  ,JZL  _  lPg,+  ~Pg,- 
dt  dx  p  y 


+py  1 +  }dx2  +  2  dx2 


where  /+  and  /_  are  given  by 


_  dy  _dy 

V=ft+Udi 


f  =  [i  +  (!Hl\2  i  (?i\2 

J±  [dx  dx  dx  4\dx 


with  additional  term  appearing  on  the  right-hand-side  of  Eqs.  (3,  4)  due  to  the  nonzero  gas 


pressures  p9t+  and  pg-  along  the  upper  (+)  and  lower  (— )  interfaces.  In  the  above  equations 
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y  is  the  sheet  thickness,  y  is  the  sheet-centerline  location  in  the  transverse  direction,  and  u 
and  v  denote  the  axial  and  transverse  film  velocity  averaged  over  the  film  thickness.  The 
underbraced  terms  in  the  original  Eqs.  (3,  4)  identify  nonlinear  coupling  terms  between 
sinuous  and  dilational  film  disturbances  (or  capillary  waves)  for  the  general  nonlinear  case 
resulting  from  capillary  forces.  Note  that  in  the  nonlinear  case  }±  depends  on  both  y  and 
y  (see  Eq.  6).  Mode  coupling  is  also  influenced  by  the  surrounding  gas  flow  represented  in 
Eqs.  (3)  and  (4)  by  the  gas  pressures  p3)+  and  pg-  along  the  gas-liquid  interfaces.  Eqs.  (2) 
through  (5)  are  valid  for  thin  inviscid  incompressible  liquid  films.  Here  “thin”  implies  that 
the  film  thickness  h  is  small  compared  to  the  wavelength  of  any  disturbance  appearing  on 
the  film  interfaces.  A  concise  derivation  of  Eqs.  (2  -  5)  has  been  presented  in  Ref.  39  and  in 
Ref.  37,  as  well  as  in  Ref.  40  for  the  more  general  case  of  thin  annular  films.  For  dilationally 
distorting  films,  we  have  y  =  0,  v  =  0  and  p9t+  =  pg-  so  that  Eqs.  (4)  and  (5)  are  identical 
to  zero  with  f±  =  [1  +  (dy/dx)2/ 4]-3/2  in  Eq.  (3).  Eqs.  (2,  3)  and  Eqs.  (4,  5)  are  decoupled 
by  linear  analysis,  governing  linear  dilational  mode  disturbances  (i.e.  u  and  y)  and  linear 
sinuous  mode  disturbances  (i.e.  v  and  y),  respectively. 

The  unsteady  Bernoulli  equation  provides  a  relationship  between  the  pressure  in  the  gas 
phase  pgii,2  and  the  velocity  potentials  <f>g;ip  governed  by  V2<pgiip  =  0: 


1,2 

dt 


+  2  (^s; i.s)2  -  “  (Pg;  1,2 


(7) 


where  p°.12  =  fctn(i)  is  the  stagnation  pressure  in  the  gas  below  (subscript  1)  and  above 
(subscript  2)  the  liquid  film.  The  gas  densities  on  both  sides  ,  p9ii  and  p9t 2,  were  assumed  to 
be  equal;  the  static  pressure  at  the  outflow  boundary  has  been  prescribed  at  a  constant  value 
(referenced  at  0  Pa).  Consequently  p°;12  =  Eq-  (7)  since  =  0-  Evaluation 

of  Eq.  (7)  at  the  liquid-gas  interfaces  governs  the  gas  pressure  p9t±.  The  adjustment  of 
P°g-,  1,2ft)  is  employed  in  order  to  prevent  transverse  film  deflection  at  the  outflow  boundary, 
particularly  for  non-symmetric,  i.e.  mixed  or  sinuous  mode  film  pulsing  due  to  the  non-zero 
pressure  difference  across  the  undisturbed  film.  Note  that  properly  prescribed  variations 
in  the  outflow  static  pressure  or  total  pressure  in  both  gas  regions  could  be  used  to  model 
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the  effects  of  combustion  instabilities  (i.e.  acoustic  resonant  oscillations)  on  the  breakup 
of  the  injected  liquid  fuel  film.  Such  an  analysis  has  already  been  presented  by  Hilbing  et 
al.41  for  finite  length  liquid  jets  using  a  boundary-element  method  to  describe  both  liquid 
and  gas  phases.  A  similar  analysis  for  the  case  of  semi-infinite  thin  planar  films  would 
require  modifications  of  the  specified  outflow  conditions  and  possibly  inflow  conditions  for 
the  discharging  liquid  film. 

The  kinematic  interface  condition  governing  the  motion  of  the  phase  interface  is  implicit 
in  Eq.  (2)  which  governs  mass  continuity  for  the  liquid  phase. 

Prom  Eq.  (7)  we  observe  that,  in  order  to  evaluate  the  local  gas  pressure,  d^g-i^/dt  has 
to  be  known  along  the  interface.  Here,  the  local  values  of  this  time-derivative  are  obtained 
explicitly  using  previous  values  of  at  the  interface  and  at  the  same  downstream  location. 
This  corresponds  to  the  analysis  presented  by  Spangler,  Hilbing  and  Heister,42  who  analyzed 
the  two-dimensional  planar  sheet  problem  by  using  a  BEM  for  the  solution  of  Laplace’s 
equation  in  both  gas  and  liquid  phases.  Those  authors  tracked  the  motion  of  the  nodes  at 
the  liquid-gas  interface  only  in  the  direction  vertically  to  the  liquid  stream  or  jet. 

Eqs.  (2)  through  (7)  have  to  be  integrated  in  time  which  also  includes  the  simultaneous 
solution  of  Laplace’s  equation  for  (f>g] i  and  (f)g-2  at  each  time  step.  The  procedure  begins 
with  the  solution  of  Laplace’s  equation  for  <f>g-,i,2  by  using  the  current  interface  locations 
and  the  velocity  normal  to  the  interfaces,  i.e.  d<j>g;ifl/ dn.  The  solution  gives  the  values 
of  the  velocity  potentials  <j>g.x  and  03;2  along  the  two  interfaces.  This  information  permits 
the  updating  of  the  gas  pressures  p3)±  using  Eq.  (7).  The  newly  determined  gas  pressure 
pan  now  be  used  within  Eqs.  (3)  and  (4)  which  in  combination  with  Eqs.  (2,  5)  is  solved 
to  update  the  interface  locations  (including  the  locations  for  the  BEM  surface  nodes)  and 
the  velocity  components  normal  to  the  interfaces.  By  repeating  this  procedure  at  each  time 
step,  the  film  shape  can  be  determined  at  all  times  prior  to  droplet  pinch-off.  The  prescribed 
solution  procedure,  illustrated  schematically  in  Fig.  2  for  the  dilational  case,  closely  follows 
the  procedure  employed  by  Spangler,  Hilbing  and  Heister42  who  analyzed  the  nonlinear 
evolution  of  an  axisymmetric  liquid  jet  by  using  a  boundary-element  method  (BEM)  for  the 
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solution  of  Laplace’s  equation  in  both  gas  and  liquid  phases  and  by  employing  the  unsteady 
Bernoulli  equation  in  both  phases  along  the  interface. 

Note  that  the  solution  procedures  employed  for  the  gas-phase  (discrete  BEM)  and  the 
liquid  phase  (reduced-dimension  analysis)  solve  for  nondimensional  quantities.  Both  solution 
methods  use  the  thickness  of  the  film  as  characteristic  length  for  the  nondimensionalization. 
The  characteristic  velocity  for  the  reduced-dimension  analysis  of  the  liquid  film  is  the  capil¬ 
lary  velocity  u3  =  ^cr/(2pih)  which  eliminates  the  Weber  number  from  the  nondimensional 
forms  of  Eqs.  (3,  4).  u3  is  the  velocity  at  which  small  amplitude  sinuous  capillary  waves 
travel  along  a  planar  liquid  film.37  For  the  BEM  solution  of  the  gas-phase  velocity  poten¬ 
tial,  and  the  unsteady  Bernoulli  equation  in  the  gas  phase,  the  maximum  gas-jet  injection 
velocity  is  used  as  the  characteristic  velocity. 

With  the  prescribed  non-dimensionalization,  the  following  (nondimensional)  flow  pa¬ 
rameters  are  to  be  considered:  (1)  The  ratio  of  maximum  gas  jet  momentum  to  liquid  film 
momentum  M  =  pgu2, 0s/(piu2ti),  (2)  the  nondimensional  square  of  the  film  injection  veloc¬ 
ity  or  liquid-phase  Weber  number  Wej  =  piufh/a,  (3)  the  gas-to-liquid  density  ratio  pg/pu 
and  three  geometric  parameters,  i.e.  the  gas  jet  injection  angle  a,  as  well  as  s/h  and  H/h 
denoting  the  ratios  of  gas-jet  inlet  port  width  or  discharge  channel  width  to  the  thickness 
of  the  undisturbed  liquid  film,  respectively.  Within  the  present  analysis  the  parameters  M 
and  a  have  been  replaced  by  M||  =  M  cos  a  and  Mj_  =  M  sin  a  representing  the  maximum 
values  of  the  gas  jet  components  parallel  and  perpendicular  to  the  undisturbed  liquid  film. 
The  former  is  relevant  for  the  initial  generation  of  film  disturbances,  while  the  latter  can 
cause  amplification  of  an  existing  disturbance  through  Kelvin-Helmholtz  wave  growth. 

From  Eqs.  (3,  4)  we  see  that,  within  the  considered  inviscid  model,  the  gas-phase  influ¬ 
ences  the  liquid-phase  dynamics  solely  through  the  static  pressure  at  the  interfaces  pg<±.  Its 
nondimensional  value  p*  ±  =  p9l±/(pi  u2s)  within  the  nondimensional  form  of  Eqs.  (3,  4)  can 
be  expressed  in  terms  of  the  nondimensional  pressure  pft±  =  p9t±/ (0.5  pgu2.0)  obtained  from 
the  solution  of  the  unsteady  Bernoulli  equation  along  the  interface, 
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(8) 


_*  =h8i  ?i® 

Psi±  2  Pi  U* 


The  nondimensional  parameter  of  interest  in  this  expression  is  the  product  formed  by  the 
gas-to-liquid  density  ratio  pgf  pi  and  the  ratio  between  the  maximum  gas  jet  velocity  ugto  and 
the  capillary  velocity  us  on  the  liquid  film.  For  practical  applications  this  product  is  of  the 
order  one.  In  particular,  for  a  water  sheet  with  a  thickness  of  0.2mm  at  20°  C  (i.e.  a  =  0.0727 
N/m  and  pi  =  998  kg/m3),  and  a  co-flowing  gas  stream  with  pg  —  1.2  kg/m3  moving  at  15 
m/s,  the  prescribed  nondimensional  value  is  pgv^.0/ (piuf)  =  1.5.  pgu2gQJ (piuf)  =  2M  Wej  h/ s 
represents  the  ratio  of  aerodynamic  forces  exerted  on  the  liquid-gas  interface  to  capillary 
forces.  As  a  parameter  relevant  to  the  gas/liquid  field  governing  equations,  it  does  not 
contain  any  geometrical  parameters  with  exception  of  the  film  thickness.  Due  to  the  different 
choices  made  for  the  characteristic  velocities  in  the  nondimensional  equations  for  the  gas 
and  liquid  phases,  ug;o/u3  is  another  nondimensional  parameter  relevant  to  the  problem. 
However,  since  u9-q/u3  can  be  expressed  in  terms  of  We;,  M,  s/h  and  pg/pi,  the  density 
ratio  pglpi  has  been  chosen  as  the  second  (geometry  independent)  nondimensional  quantity 
relevant  to  the  field  equations.  Note  that,  if  u9]q  ^  u„,  the  same  nondimensional  time-steps 
for  the  numerical  procedures  employed  to  solve  the  liquid-  and  gas-phase  equations  implies 
different  physical  time  steps.  Consequently,  and  since  the  capillary  velocity  us  is  in  general 
smaller  than  the  velocity  ug- o,  subiteration  for  the  solution  of  the  liquid-phase  equations  is 
needed  (if  the  same  nondimensional  time-step  is  used)  or  the  nondimensional  time-step  for 
the  liquid-phase  analysis  has  to  be  adjusted  to  the  one  employed  for  the  BEM  solution  of 
the  gas  phase.  Note  that  parameters  pgJ pi  and  pgu^.Q/ (piuf)  are  essential  to  the  problem  or 
field  equations,  while  other  nondimensional  parameters  such  as  a,  s/h  and  H/h  depend  on 
the  particular  atomizer  configuration. 


Model  Validation  and  Numerical  Accuracy 

The  various  components  of  the  numerical  model  employed  within  this  analysis  and  dis¬ 
cussed  above  have  been  benchmarked  independently.  Implementation  of  the  boundary- 
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element  method  has  been  validated  by  predicting  the  flow  across  a  cylinder  and  a  NACA 
airfoil  and  by  comparing  the  computed  pressure  distributions  with  analytical  and  experimen¬ 
tal  results,  respectively.  Here,  the  pressure  across  the  different  shapes  has  been  evaluated 
by  using- the  steady-state  Bernoulli  equation.- The  liquid-phase  model  (without  gas-phase 
coupling)  has  been  analyzed  in  detail  within  previous  work,  including  an  analysis  on  the 
_accuracy  of  the  employed  thin  film  approximation  or  long  wavelength  assumption  close 
to  the  point  of  film  pinch-off  where  short  wavelength  contributions  can  be  present.  The 
latter  has  been  accomplished  by  comparison  of  the  computed  results  for  an  infinite  period¬ 
ically  disturbed  liquid  film  with  solutions  obtained  from  a  two-dimensional  discrete-vortex 
method.37’43  Accuracy  in  the  evaluation  of  the  unsteady  term  within  Eq.  (7)  has  been  estab¬ 
lished  independently  by  comparing  the  pressure  distribution  along  the  surface  of  an  object 
placed  into  a  gas  stream  with  the  pressure  distribution  predicted  for  the  same  object  now 
moving  with  the  same  relative  velocity  in  the  same  but  quiescent  gas.  With  respect  to  the 
pressure  evaluation  within  the  gas-phase,  the  latter  case  is  unsteady  (due  to  the  moving 
boundary)  and  d<f>g/&t  ^  0,  whereas  in  the  former  (steady)  case,  d(f>g/dt  —  0. 

Note  that  the  angled-derivative  scheme45  employed  here  and  in  Ref.  40  for  the  solution 
of  the  liquid  film  equations  is  second-order  accurate  in  time  and  space.  Numerical  solutions 
for  the  pure  convection  equation  showed  that  this  scheme  is  superior  to  the  Lax-Wendroff 
Method  with  Richtmyer  splitting  previously  employed  in  Ref.  44,  as  it  produces  significantly 
less  numerical  diffusion  for  this  test  case.  The  latter  feature  is  essential  if  details  of  the  capil¬ 
lary  wave  propagation  along  the  liquid  film  are  to  be  predicted  accurately.  More  information 
on  the  numerical  method  can  be  found  in  Ref.  45. 

III.  RESULTS 


Base  Case  Analysis 

Fig.  3  shows  a  dilationally  distorting  liquid  film  subject  to  the  forcing  conditions  imposed 
onto  the  gas  jet  inlet  velocities  tts;i  and  ug. 2  according  to  Eq.  (1)  at  the  lower  (1)  and  upper 
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(2)  gas-inlet  ports,  respectively.  Here,  both  lower  and  upper  gas  jet  pulsations  are  in  phase 
so  that  ngy\  =  ng>2  =  0.  This  pulsing  sequence,  i.e.  one  gas  pulse  of  time  period  Tp 
followed  by  a  ‘resting-period’  Tp  without  gas  flow  is  subsequently  described  as  the  base- 
case  pulsation.  Problem  parameters  for  this  case,  subsequently  referred  to  as  base-case 
conditions  are  M||  =  Mx  =  1,  p  =  pg/pi  =  0.0012,  We;  =  10  and  s/h  =  H/h  =  10.  For  all 
results  presented  in  this  analysis,  grid-size  and  time-step  independence  has  been  verified  by 
successive  time-step  reduction  and  mesh  refinement  (i.e.,  addition  of  computational  nodes 
within  the  boundary-element  analysis).  ~  •  .y  • 

Fig.  4  shows  the  dilationally  distorting  liquid  film  for  the  similar  case  but  a  modified 
gas-jet  modulation.  Here,  gas  jet  pulses  of  period  Tp  follow  immediately  upon  each' other 
without  a  ‘resting-period.’  As  observed  from  Fig.  3,  the  energy  provided  by  the  impacting 
gas  jets  and  base-case  pulsation  is  insufficient  to  break  the  film,  the  latter  being  stabilized 
by  surface  tension.  As  more  energy  is  transferred  onto  a  given  length  of  the  moving  film, 
this  stabilizing  effect  is  overcome  and  film  rupture  is  predicted  at  a  downstream  location  of 
x/h  =  106.5  and  at  a  nondimensional  time  t*  —  25. 

Assuming  base-case  parameter  conditions  as  well  as  base-case  forcing  conditions  for  both 
upper  and  lower  gas  jets,  but  allowing  for  a  phase  shift  or  time  shift  of  period  Tp  between  both 
jet  pulses,  i.e.  ngti  =  0  and  ngy2  =  1,  one  obtains  an  initially  sinuous  distorting  film  as  shown 
in  Fig.  5.  Note  that  changing  the  pulse  shape  from  a  sin2-function  to  a  (1  —  cos)  distribution 
while  keeping  the  pulse  period  Tp  and  the  overall  pulse  energy  constant  does  not  significantly 
alter  the  resulting  film  distortion.  Recall  that,  for  all  the  simulations  presented  here,  plug 
flow  has  been  assumed  at  both  gas  jet  outlets.  The  effect  of  start-up  conditions  has  been 
studied  for  the  sinuous  case  of  Fig.  5.  Fig.  6  displays  the  result  for  the  case  where  the  first 
gas-jet  pulse  (entering  the  discharge  channel  from  the  lower  port)  was  specified  according  to 
u'g ;1  =  (1  —  exp  [— t/Ts\)  ug-i(t)  with  u3;1  from  Eq.  (1)  and  Ta  =  2.5.  As  observed  from  Fig.  6, 
the  downstream  location  where  film  rupture  occurs  does  not  change  significantly  although 
break-up  times  do  vary.  For  the  case  shown  in  Fig.  6,  the  time  until  film  rupture  occurs 
increased  from  approximately  t*  =  16  to  t*  —  21.  However,  film  pinch-off  is  still  predicted 
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at  x/h  «  45,  although  the  location  of  the  pinch-off  within  the  waveform  differs.  Variation 
of  ne#  between  0  and  1  (with  n5;1  =  0)  produces  mixed-mode  film  deformations  (see  Fig.  7). 
Due  to  the  complexity  in  the  resulting  interface  topology,  the  present  work  focused  on  the 
analysis  of  dilational  and  sinuous  gas  jet  pulsations  only. 

For  the  cases  shown  in  Figs.  5  and  6,  strong  nonlinear  coupling  between  sinuous  and 
dilational  modes  leads  to  film  breakup  close  to  the  nozzle  exit.  Note  that,  on  average, 
the  same  amount  of  energy  is  added  to  the  system  as  previously  for  the  case  with  in- 
phase  pulsing  of  both  jets  (see  Fig.  3).  However  with  the  given  forcing  frequency  and 
atomizer  configuration,  film  rupture  is  more  readily  obtained  indirectly  by  sinuous  forcing 
and  nonlinear  sinuous-dilational  mode  coupling  rather  than  by  forcing  film  pinch-off  directly 
by  dilational  film  modulation. 

Comparison  of  Figs.  4  and  5  also  shows  that  sinuous  gas-jet  pulsation  under  base-case 
conditions  (and  with  base-case  forcing)  yields  even  shorter  breakup  times  and  lengths  than 
the  dilational  case  with  modified  forcing  conditions,  i.e.  without  ‘resting-period;  ’  despite  the 
fact  that,  for  the  modified  dilational  case  of  Fig.  4,  energy  input  into  the  system  per  unit 
time  is,  on  average,  twice  that  for  the  sinuous  case  illustrated  in  Fig.  5.  Clearly,  the  greater 
amount  of  energy  provided  at  the  gas-jet  inlet  ports  will  not  necessarily  result  in  shorter  film 
break-up  times  and  lengths.  More  importantly,  time  history  and  absolute  values  of  energy 
transfer  rate  from  gas-  to  liquid-phase  will  determine  the  effectiveness  with  which  the  gas 
jets  cause  large  amplitude  dilational  film  distortions,  which  ultimately  cause  film  rupture. 

Sinuous  (i.e.  out-of-phase)  gas-jet  modulation  will  result  in  larger  differences  between 
the  instantaneous  gas-phase  pressures  on  both  sides  of  the  liquid  film.  Consequently,  more 
energy  is  transferred  onto  the  liquid  film  than  under  dilational  (i.e.  in  phase)  forcing  con¬ 
ditions  for  the  two  gas  jets.  If  nonlinear  sinuous-dilational  mode  coupling  is  strong  (e.g., 
for  low  We;  values),  sinuous  gas-jet  forcing  conditions  will  consequently  result  in  faster  film 
break-up  than  dilational  gas-jet  forcing.  The  importance  of  nonlinear  sinuous-dilational 
mode  coupling  as  the  determining  mechanism  for  film  break-up  has  already  been  discussed 
for  annular  liquid  films  discharging  into  a  surrounding  void.40 
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Weber  Number  Effect 

Under  conditions  similar  to  those  in  Fig.  5  but  larger  Weber  number  We;,  film  resistance 
to  transverse  distortion  or  deflection  is  reduced  (see  Fig.  8  for  We*  =  25);  however,  at  larger 
-  We/- values,  nonlinear  sinuous-dilatienal-mode  coupling -remains  weak  and  the  film  does  not 
rupture  before  it  hits  the  walls  of  the  discharge  channel.  Although  impact  of  the  film  with 
Jrhe  wall  should  result  in  rupture,  the  mechanism  and  the  results  will  differ  from  those  studied 
here.  For  our  purposes,  we  will  consider  impact  with  the  walls  to  be  undesirable  although 
it  might  have  useful  consequences  in  some  technological  applications. 

Reduced  mode  coupling  is  manifested  by  stretched-band-like  film  distortion  rather  than 
the  formation  of  fluid  blobs  (cylinders)  as  found  for  We;  =  10  (see  Fig.  5)  and  as  illustrated 
also  in  Fig.  9  for  the  case  with  We/  =  5.  In  the  latter  case,  the  smaller  Weber  number 
(e.g.  larger  surface  tension  at  same  discharge  velocity  and  film  thickness)  stabilizes  the  film 
by  reducing  transverse  film  oscillations  and  consequently  also  gas-phase  pressure  variations 
along  the  film  in  the  downstream  direction.  This  results  in  larger  break-up  lengths  and 
times  for  the  smaller  Weber  number  case  despite  the  fact  that,  for  a  given  sinuous  wave 
(i.e.  wavelength  and  disturbance  amplitude),  sinuous-dilational  mode  coupling  is  stronger 
for  the  smaller  Weber  number  (i.e.  We/  =  5)  as  illustrated  by  the  underbraced  terms  in 
Eqs.  (3,4). 

Previous  analyses  of  planar  films  in  a  zero  ambient  gas  and  with  film  modulations  en¬ 
forced  locally  onto  the  discharging  film  at  the  nozzle  exit  showed  that,  at  a  given  imposed 
forcing  frequency,  up  to  two  sinuous  mode  wavenumbers  are  observed  downstream  form  the 
nozzle.  In  the  parameter  range  of  interest  (i.e.  large  Weber  number  and  low  density  ratio) 
both  of  these  wavenumbers  vary  only  little  in  magnitude  so  that  the  resulting  film  distortion 
shows  a  characteristic  beat  behavior.37  Such  a  beat  behavior  could  not  be  identified  within 
the  present  configuration  where  film  modulation  is  enforced  via  impacting  and/or  co-flowing 
gas  streams.  Note  however  that,  in  the  present  analysis,  gas  jets  of  a  given  finite  width 
impact  the  moving  film  with  a  specified  pulsation  pattern  and  at  a  given  frequency.  In  other 
words,  the  various  disturbances  generated  by  these  jets  and  their  specific  wavelengths  and 
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propagation  velocities  not  only  depend  on  the  imposed  gas-jet  modulation  frequency  but  also 
on  the  gas-jet  geometry  or  gas-jet  width;  this  is  in  contrast  to  the  analysis  of  Ref.  37.  For 
direct  film  modulation  at  the  nozzle  and  zero  density  ambient  gas  (in  Ref.  37),  disturbances 
of  various  wavelengths  result  naturally  from  the  harmonic  modulation  of  the  flow  variables 
at  the  nozzle  exit.  The  resulting  wavelengths  are  affected  by  the  imposed  forcing  frequency 
and  the  liquid  Weber  number  (based  on  film  thickness  and  film  discharge  velocity)  only. 

This  raises  the  question,  if  film  rupture  via  gas-jet  modulation  can  be  optimized  by 
co-ordinating  gas-jet  width  and  frequency  of  gas-jet  pulses  (at  a  given  nondimensional  film 
discharge  velocity  or  Weber  number  We*).  The  difficulty  here  is  that  any  modification  of 
the  pulsing  frequency  will  not  only  alter  film  distortion  directly  by  altering  the  film  length 
exposed  to  the  impacting  gas  jet  during  the  pulse  but  also  indirectly  by  modifying  the  overall 
gas-phase  flowfield. 

With  respect  to  this  study  of  Weber  number  effects  on  liquid  film  distortion,  it  should  be 
noted  that  an  increase  in  Weber  number  can  be  interpreted  as  a  decrease  in  surface  tension 
at  the  same  value  of  h  and  ui  or  as  an  increase  of  itj  at  constant  surface  tension  and  h- value. 
In  the  latter  case,  the  reference  time  within  the  solutions  shown  in  Fig.  5  and  8  are  the  same, 
resulting  in  the  same  dimensional  time  frame.  However,  in  the  former  case,  tre/  =  /re//itre/ 
and  uref  =  ^Jcx/(2pih)]  therefore  an  increase  in  the  surface  tension  coefficient  a  yields  an 
increase  in  tref  so  that  for  the  same  nondimensional  time,  the  elapsed  physical  time  is  larger 
in  Fig.  8  than  in  Fig.  5.  For  constant  liquid  velocity  ut,  we  have  uref  ~  We~lj/2,  so  that 
the  dimensional  time  scale  in  Fig.  8  (with  Wej  =  25)  is  stretched  by  a  factor  of  y/23  with 
respect  to  the  dimensional  time  in  Fig.  5  (with  Wej  =  10). 

Fig.  10  shows  dilational  film  deformation  for  the  same  parameter  configuration  as  used  in 
Fig.  8,  only  now  i  =  ngi 2  =  0,  whereas  before  ngt  1  =  0  but  ng^  =  1.  Figs.  8  and  10  suggest 
that,  for  larger  Weber  numbers  (and  film  modulation  via  gas  jet  pulsing),  film  rupture  can 
be  effectively  controlled  by  dilational  gas-jet  pulsations,  particularly  if  constraints  of  the  flow 
geometry  are  to  be  considered.  In  addition,  for  large  Weber  number  flows,  sinuous-dilational 
mode  energy  transfer  becomes  increasingly  ineffective  so  that  film  rupture  lengths  for  similar 
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dilational  and  sinuous  mode  forcing  conditions  become  comparable. 

In  what  follows,  the  discussion  will  focus  on  parametric  studies  for  the  sinuous  mode 
case  of  Fig.  5.  Sinuous  mode  film  distortion  has  previously  been  considered  to  be  the  most 
-relevant  mode  of  film  disintegration  for  practical  atomization  systems.  This  does  not  imply 
that  the  sinuous  mode  will  be  more  effective  in  causing  film  rupture  than  the  dilational  mode 
jinder  every  possible  operating  condition  (e.g.,  for  high-pressure  or  large  density  ratio  pg/pi 
applications).  In  the  past,  experimental  studies  have  focused  on  the  sinuous  mode  of  film 
disintegration,  while  numerical  models  considered  (with  some  exceptions)  the  analysis  of  the 
dilational  mode  only.  Based  on  this  observation,  the  present  analysis  is  intended  to  provide 
some  insight  to  film  distortion  under  the  mode  least  considered  within  previous  numerical 
or  analytical  work  on  nonlinear  liquid  film  distortion  and  disintegration. 

Gas  Momentum  Effect 

Several  variations  and  combinations  of  gas-jet  inlet  momentum  My  and  Mx  have  been 
considered  for  the  case  shown  in  Fig.  5  where  My  =  M_l  =  1.  The  result  obtained  for 
Mj_  =  0.5  is  illustrated  in  Fig.  11.  According  to  the  reduced  vertical  momentum,  the 
effect  of  the  impacting  gas  jet  onto  the  film  is  reduced,  resulting  in  smaller  initial  sinuous 
film  disturbances.  Film  break-up  after  downstream  growth  of  these  disturbances  (possibly 
enhanced  by  capillary  and/or  Kelvin-Helmholtz  instabilities)  is  delayed.  A  result  similar  to 
the  one  shown  in  Fig.  11  was  obtained  for  the  case  where  M  =  ^/Mjj  +  Mj_  is  kept  constant 
while  a  is  reduced  from  45  degrees  (see  Fig.  5)  to  30  degrees  (see  Fig.  12).  Mx  >  0.5 
and  My  >  1  in  these  cases.  However,  despite  the  larger  values  of  Mx  and  My,  maximum 
amplitudes  at  a  given  time  and  downstream  location  remain  somewhat  smaller.  Note  that, 
the  physical  relevance  of  the  short  wavelength  disturbances  observed  in  Fig.  12  at  t*  =  27 
for  x/h  >  75  as  predicted  by  the  employed  thin  film  model  is  uncertain.  The  accuracy  and 
limitations  of  the  current  liquid-phase  model  with  respect  to  short  wavelength  contributions, 
particularly  close  to  the  point  of  film  breakup,  has  been  discussed  in  detail  in  Ref.  44. 

If  a  is  decreased  further  (i.e.  to  a  =  15°),  gas-jet  normal  momentum  is  further  reduced 
in  favor  of  an  increase  in  gas  momentum  parallel  to  the  undisturbed  film.  Fig.  13  shows  that, 
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as  a  result,  amplitudes  of  sinuous  mode  film  disturbances  are  somewhat  reduced  compared 
to  those  found  in  Fig.  12  at  the  same  time  and  downstream  location.  Results  obtained  for 
a  =  75°  varied  only  slightly  from  those  computed  for  a  =  45°  shown  in  Fig.  5. 

~  Figs.  14  and  15  illustrate  cases  where  the  vertical  or  parallel  component  of  the  gas-jet 
momentum  was  increased  from  1  to  5  (at  both  inlet  ports)  while  keeping  all  other  parameters 
the  same.  In  both  instances,  film  pinch-off  occurs  close  to  the  nozzle  exit  as  a  result  of  the 
pressure  distribution  generated  by  the  impacting  gas  jets  rather  than  through  continuous 
wave  growth,  the  latter  being  relevant  for  the  case  shown  in  Fig.  11  and  13,  for  example. 

Note  that  continuous  wave  growth  in  Figs.  11  and  13  does  not  imply  monotonic  growth. 
The  monotonic  growth  can  be  expected  for  steady  gas-phase  boundary  conditions  (i.e.  zero 
or  constant  velocity  at  gas  inlet  ports).  For  gas-jet  pulsations  without  underlying  base 
flow,  relative  velocities  between  gas  jets  and  liquid  film  vary  significantly  over  time  causing 
disturbance  amplitudes  to  growth  sporadically. 

Density  Effects 

Changes  in  the  density  ratio  pgj pi  from  the  case  shown  in  Fig.  5  have  also  been  analyzed. 
Fig.  16  shows  the  result  obtained  for  the  similar  parameter  set  as  in  Fig.  5  but  with  a  five-fold 
increase  in  pg/ pi .  Due  to  the  relative  increase  of  gas-phase  inertia,  transverse  film  deflection 
is  now  larger  (at  the  same  gas-jet  forcing  conditions).  Film  distortion  as  illustrated  in 
Fig.  16  resembles  the  prediction  made  in  Fig.  8  for  an  increase  in  liquid  Weber  number  Wei. 
However,  in  the  present  case,  mode  coupling  is  not  altered  and  the  film  ruptures  shortly 
before  impacting  on  the  channel  walls.  Also,  for  the  larger  density-ratio  case,  film  deflection 
extends  over  a  larger  domain  than  shown  in  Fig.  8. 

Pulse  Period 

The  effect  of  variations  in  pulse  period  Tp  on  film  distortion  for  the  sinuous  mode  of 
Fig.  5  will  now  be  summarized:  1)  For  large  Tp- values  (and  constant  jet  width),  the  trans¬ 
verse  displacement  of  the  film  becomes  large  before  the  counter-acting  gas  jet  causes  the 
deceleration  and  reversal  of  the  displacement.  Here,  transverse  film  motion  caused  by  an 
impacting  gas-jet  pulse  might  result  in  film  impingement  onto  the  channel  walls  (before  re- 
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versal  of  the  film  dislocation  due  to  gas-jet  impact  onto  the  film  from  the  other  side).  2)  For 
small  Tp-values  (see  Fig.  17),  transverse  film  distortion  remains  small  (due  to  more  rapid 
appearance  of  the  opposed  gas  jets)  with  only  shorter  wavelength  disturbances  appearing 
on  the  film.  Note  that  the  gas  jet  width  is  s/h  =  10  and  gas-jet  impact  changes  from 
one  side  of  the  film  to  the  opposite  side  at  frequency  1/TP.  Therefore,  even  though  one 
gas  jet  generates  an  instantaneous  film  disturbance  of  approximate  length  l  =  (s/h)  cos  or, 
the  only  part  of  this  disturbance  not  compensated  by  the  action  of  the  opposed  jet  has  a 
length  of  A  =  uj  Tp.  In  order  to  achieve  an  effective  initial  sinuous  mode  film  distortion  with 
gas-jets  phased  according  to  n^i  =  0  and  —  1,  the  pulsing  period  Tp  should  be  longer 
than  the  time  needed  for  a  liquid  particle  to  travel  across  the  length  of  the  gas  inlet  ports, 
Le.  Tp>  s/[h  cosattj).  However,  as  mentioned  above,  for  very  large  Tp  values  the  film  will 
tend  to  impact  onto  the  walls  of  the  discharge  channel.  In  the  present  configuration,  gas-jet 
impact  is  confined  to  the  vicinity  of  the  atomizer,  whereas  further  downstream,  amplitude 
growth  of  the  liquid  film  disturbance  will  be  increasingly  influenced  by  Kelvin-Helmholtz 
effects.  Therefore,  the  pulse  period  providing  optimum  wave  growth  will  depend  on  1)  the  - 
generation  of  an  optimum  (to  be  defined)  initial  film  disturbances  in  the  vicinity  of  the  at¬ 
omizer  (i.e.,  large  amplitude  and  large  initial  dilational  mode  contribution  or  film  straining), 
and  2)  amplification  of  the  initial  disturbances  through  (sporadic)  Kelvin-Helmholtz  wave 
growth,  where  the  wavelength  of  the  disturbance  with  optimum  growth  rate  will  depend 
on  the  given  parameter  configuration,  e.g.  ui ,  uSii/2(t),  We;,  pg/pi ,  and  the  dynamics  of  the 
unsteady  gas  streams. 

Gas-Phase  Continuous  Flow 

The  previous  discussion  focused  on  liquid  film  distortion  and  disintegration  resulting  from 
time-periodic  modulation  of  the  gas-phase  velocity  on  both  sides  of  the  liquid  film.  Forcing 
conditions  onto  the  initially  quiescent  gas  phase  were  imposed  at  the  gas-jet  outlets  according 
to  Eq.  1.  Accordingly,  gas-phase  injection  remained  confined  to  pulse  events  of  time-period 
Tp.  Between  these  injection  events,  gas-phase  velocities  on  both  sides  of  the  liquid  film  were 
determined  by  the  dynamics  of  the  deforming  gas-liquid  interfaces.  In  this  paragraph,  gas- 
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jet  pulses  are  superimposed  onto  continuous  gas  streams  discharging  from  the  two  gas-inlet 
ports.  Figs.  18,  19  and  20  illustrate  3  cases  with  time-dependent  forcing  conditions  similar 
to  the  sinuous  mode  case  of  Fig.  5  now  superimposed  onto  steady  gas  streams  entering  from 
the  two  gas  inlet  ports  at  constant  velocity  ug;3.~  Steady  inlet  velocities  are  assumed  to  be 
the  same  at  both  gas-inlet  ports  in  order  to  prevent  continuous  transverse  film  deflection, 
the  latter  being  limited  by  the  discharge  chamber  walls.  Pulse  period  Tp  in  Figs.  18  through 
20  is  Tp  =  5, 8.75  and  20,  respectively;  gas-stream  steady  base-flow  velocity  and  amplitude 
of  the  gas-phase  velocity  modulation  are  fixed  at  ug]3  =  0.25  and  ug<o  =  0.5.  Therefore,  the 
maximum  injection  velocity  reached  by  the  gas-jets  in  Figs.  5  and  18  is  the  same,  while  the 
average  energy  flux  provided  by  the  gas-jets  in  the  latter  case  is  larger,  due  to  the  underlying 
time-independent  velocity  component.  The  average  energy  flux  for  gas  inlet  ports  1  and  2 
is  given  by 


1.2  = 


n 


(9) 


with  the  mass  influx  per  unit  width  of  the  two-dimensional  gas  inlet  ports  m\p(t)  = 
pg  s  ug;ip(t).  The  integration  in  Eq.  (9)  is  performed  over  a  multiple  number  of  pulse  periods 
Tp  (n  =  even  integer).  The  time-independent  component  of  the  overall  gas-jet  momentum 
is  equally  distributed  onto  the  gas  jets  on  either  side  of  the  liquid  film;  consequently,  the 
fluctuating  part  of  the  gas-jet  momentum  or  gas  pulse,  which  is  reduced  from  its  original 
value,  will  cause  less  deflection  of  the  liquid  film  in  the  transverse  direction.  This  results 
in  an  increase  in  breakup  length  and  time  (see  Fig.  18),  even  though  the  underlying  non- 
modulated  gas-phase  velocity  allows  for  continuous  energy  transfer  from  gas-to  liquid  phase, 
i.e.  continuous  Kelvin-Helmholtz  wave  growth. 

It  is  useful,  in  this  context,  to  compare  film  distortion  at  t*  =  32  with  the  corresponding 
interface  locations  predicted  for  the  case  shown  in  Fig.  13,  i.e.  at  t*  =  27.  Note  that, 
within  the  simulations  of  Figs.  18  through  20,  gas-jet  velocities  on  both  sides  of  the  liquid 
film  were  increased  from  zero  to  the  specified  ug-3  value  within  a  nondimensional  start¬ 
up  time  of  5  units.  After  this  time,  gas-jet  modulation  was  initiated  beginning  with  the 
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lower  gas  jet.  Also,  in  the  case  shown  in  Fig.  13,  the  maximum  amplitude  of  the  vertical 
gas-jet  velocity  component  associated  with  the  generation  of  the  initial  film  disturbance  is 
Ug.fi  sin  (71-/6),  whereas  in  the  former  case  it  is  0.5  ug.Q  sin(7r/4).  Comparison  of  the  liquid 
film  location  for  the  two  cases  at  the  prescribed  times,  therefore  illustrates  that  wave  growth 
from  intermediate  gas-jet  pulsing  can  result  in  larger  overall  amplitude  disturbances  than 
^obtained  via  continuous  wave  growth  from  continuous  gas-jets,  the  latter  providing  larger 
overall  energy  input  into  the  system. 

The  effects  of  increasing  pulse  period  Tp  for  the  case  with  underlying  base  flow  can  be 
observed  from  Figs.  18,  19  and  20.  As  Tp  is  increased,  transverse  film  deflection  resulting 
from  gas-jet  impact  is  allowed  to  proceed  further  before  being  ‘compensated’  or  counter¬ 
acted  by  a  gas-pulse  from  the  opposite  side  of  the  film.  Consequently,  amplitudes  in  the 
centerline  disturbance  of  the  film  increase  as  Tp  is  increased.  On  the  other  hand,  larger  values 
result  in  longer  wavelength  sinuous-mode  film  disturbances.  However,  at  similar  values 
for  the  disturbance  amplitude,  longer  wavelength  sinuous-mode  waves  result  in  reduced 
dilational  mode  wave  amplitudes  due  to  a  reduction  in  nonlinear  sinuous-dilational  mode 
coupling  with  increases  in  the  wavelength  of  the  sinuous  mode  wave. 

As  for  most  of  the  results  presented  in  earlier  sections,  the  computational  domain  for  the 
solutions  in  Figs.  18,  19  and  20  extended  from  x/h  =  0  to  x/h  =  125,  with  prescribed  static 
pressure  in  the  gas-phase  at  the  downstream  location.  Simulations  for  the  similar  cases  but 
extended  computational  domain,  i.e.  x/h  =  250  showed  that,  in  the  former  case,  the  imposed 
pressure  condition  effectively  suppresses  long  wavelength  film  disturbances  generated  during 
the  start-up  phase  of  the  film  modulation.  Note  that,  this  has  not  been  observed  for  the  cases 
without  continuous  gas-flow  discussed  earlier.  (The  importance  of  the  start-up  condition  for 
the  case  without  continuous  gas  jets  has  already  been  discussed  earlier  for  the  case  illustrated 
in  Fig.  5.) 

Linear  theory  for  continuous  wave  growth  on  thin  films  via  the  Kelvin-Helmholtz 
mechanism43  shows  that,  for  the  base-case  parameter  set  considered  here,  the  described 
longer  waves  have  significantly  larger  growth  rates  than  the  ones  generated  by  the  impact- 
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ing  gas  jets.  This  is  illustrated  by  the  growth  rate  plot  shown  in  Fig.  21  for  harmonic 
sinuous  waves  r]±  =  t]0  exp  [i(A;x  —  art)]  on  a  thin  liquid  film  subject  to  two  co-flowing  gas 
streams  in  an  infinite  domain  (here  r)±  denotes  the  transverse  deflection  of  the  disturbed  liq¬ 
uid  interfaces).  Since  the  longer  wavelength  contributions  are  not  suppressed  on  the  longer 
computational  domain,  transient  effects  characterized  by  growth  of  the  longer  wavelength 
disturbances  generated  during  start-up  will  result  in  significant  transverse  deflection  of  the 
liquid  film  which  is  already  being  distorted  at  shorter  wavelength  due  to  the  periodically 
impacting  gas  jets.  See  Fig.  22  for  the  case  with  Tp  =  5.  A  similar  observation  has  already 
been  made  for  modulated  thin  annular  liquid  films  in  Ref.  40  where  transient  effects  caused 
film  collapse  or  fluid  blob  formation. 

The  longest  wavelength  suppressed  on  the  shorter  domain  but  present  on  the  longer 
domain  is  of  the  order  of  the  length  of  the  longer  domain;  this  implies  that,  for  the  case 
considered  here  (Figs.  20,  22),  the  nondimensional  wavenumber  kh  suppressed  on  the  shorter 
domain  is  of  0(0.01).  From  Fig.  21  we  see  that  the  corresponding  nondimensional  growth 
according  to  the  linear  analysis  of  Ref.  43  is  about  0.035  which  corresponds  to  a  growth  rate  of 
0.04  in  terms  of  the  nondimensionalization  used  in  this  work.  This  yields  a  threefold  increase 
in  the  initial  disturbance  amplitude  of  the  prescribed  long- wavelength  contribution  after  t*  = 
25,  if  in  fact  this  wavelength  is  not  suppressed.  (Compare  Figs.  20  and  22.)  From  Fig.  20,  one 
can  estimate  the  dominant  nondimensional  disturbance  wavenumber  to  be  approximately 
kh  =  0.25,  the  corresponding  growth  rate  according  to  Fig.  21  is  approximately  six  times 
smaller  than  for  the  described  larger  wavenumber  contribution.  It  should  be  noted  again 
that  nonlinear  sinuous-dilational  mode  coupling  is  reduced  as  the  wavelength  of  the  sinuous 
wave  increases.  Consequently,  even  though  the  longer  sinuous  waves  will  result  in  larger 
transverse  film  deflection,  they  will  be  less  effective  in  causing  film  rupture  or  break-up 
than  the  shorter  wavelength  sinuous  mode  waves  generated  by  the  imposed  periodic  gas-jet 
pulses. 
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Direct  Film  Modulation 

To  this  extend,  liquid  film  distortions  were  generated  by  gas  jets  impacting  onto  the 
liquid  film  from  either  sides.  An  alternative  method  to  generate  initial  film  disturbances  is  to 
impose  forcing  conditions  directly  onto  the  liquid-phase  flow  variables  at  the  slit  nozzle  where 
the  liquid  discharges  from  the  atomizer.  Practically,  such  a  modulation  can  be  achieved  for 
^example  via  piezoelectric  actuators  placed  inside  or  adjacent  to  the  liquid  feeding  channel. 
This  has  been  analyzed  for  the  case  where  gas  density  is  negligible37.  However,  now  gas 
density  and  motion  will  have  an  effect. 

There  are  two  important  characteristics  which  distinguish  the  prescribed  direct  film 
modulation  from  the  indirect  modulation  via  impacting  gas  jets: 

1)  Only  part  of  the  energy  used  to  modulate  the  gas  jets  is  transferred  into  the  liquid 
phase  resulting  in  film  distortion.  Direct  liquid-phase  modulation,  on  the  other  hand,  guar¬ 
antees  that  all  the  energy  invested  in  modulating  one  or  more  liquid-phase  flow  parameters 
will  translate  into  film  disturbances-  Indeed,  simulations  similar  to  those  shown  in  Figs.  3 
and  5  but  without  gas  jet  modulation  and  time-dependent  forcing  conditions  imposed  onto 
the  axial  or  transverse  film  velocities  u  and  v  at  the  nozzle,  i.e. 

_  00  i  t 

w(x  =  0,  i)  =  ^  wi;q  sin2  [tt(— —  2n)]  for  2n  <  —  <  (2n  + 1) 

71—0  ^ P  -Lp 

and  w(x  =  0,  i)  =  0  ,  otherwise  (10) 

where  w  =  u  or  m  —  v,  resulted  in  nearly  instantaneous  film  rupture  if  the  energy  in¬ 
vested  into  the  modulation  is  assumed  to  equal  the  energy  added  to  the  gas-phase  for  the 
corresponding  gas-phase  modulated  case  shown  in  Figs.  3  and  5,  respectively. 

2)  Direct  liquid-phase  modulation,  as  described  above,  imposes  energy  only  locally  on  the 
film,  i.e.  at  its  discharge  location  from  the  atomizer.  Indirect  film  modulation  via  impacting 
gas  jets  causes  energy  transfer  onto  the  liquid  film  over  a  finite  length  of  the  discharging 
stream.  In  fact,  a  gas-jet  pulse  discharging  from  one  of  the  gas-jet  inlets  shown  in  Fig.  1 
will  cause  energy  transfer  or  film  modulation  over  the  entire  length  of  the  disturbed  liquid 
film;  this  is  in  contrast  to  the  discussed  direct  film  modulation  which  remains  confined 
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to  the  film  discharge  plane  throughout  the  entire  simulated  time  frame.  If,  as  considered 
later,  the  discharging  film  is  modulated  directly  in  the  presence  of  continuous  gas  streams, 
energy  transfer  onto  the  film  occurs  both  locally  at  the  nozzle  exit  and  globally  through 
interaction  of  the  disturbed  liquid  film  interfaces  with  the  “co-flowing”  gas  streams.  Even 
if  the  surrounding  gas-phase  is  quiescent,  an  originally  localized  film  disturbance  will  have 
some  non-local  effect  due  to  its  interaction  with  the  gas-phase.  The  latter  has  been  assumed 
to  be  incompressible  which  implies  an  infinite  value  for  the  speed  of  sound  or  instantaneous 
propagation  of  information  throughout  the  gas-phase. 

Fig:  23  shows  the  case  of  a  directly  modulated  liquid  film  with  or  without  steady  gas- 
jets  and  with  modulation  energy  that  equals  5  %  of  the  energy  previously  added  to  the 
gas-phase  in  Fig.  5  in  order  to  generate  sinuous  film  distortions.  In  Fig.  23,  modulations 
were  enforced  onto  the  transverse  film  velocity  v  at  the  nozzle  exit  2  =  0  according  to 
Eq.  (10).  In  other  words,  the  same  time-dependent  forcing  function  was  used  as  previously 
employed  for  the  gas-phase  modulation  in  Fig.  5.  Comparison  of  Figs.  5  and  23  shows 
that  maximum  transverse  film  deflection  is  comparable  in  both  cases,  even  though  energy 
input  for  the  direct  liquid-phase  modulated  case  is  significantly  smaller  than  for  the  indirect 
gas-phase  modulated  film.  On  the  other  hand,  as  noted  earlier,  energy  transfer  onto  the 
liquid  via  gas-jet  modulation  takes  place  along  the  entire  disturbed  interface.  This  results, 
as  observed  from  Figs.  5  and  6  in  stronger  sinuous-dilational  mode  coupling  or  dilational 
film  distortion  which  ultimately  causes  film  rupture  at  approximately  x/h  =  45.  Fig.  23 
illustrates  that  mode  coupling  without  significant  gas-phase  interaction  is  not  strong  enough 
(for  the  considered  parameter  configuration)  to  cause  film  rupture,  even  though  disturbance 
amplitudes  of  the  sinuous  waves  are  comparable  in  both  cases.  The  effect  of  continuous  gas 
jets  superimposed  onto  the  liquid-phase-modulated  discharging  film  is  also  demonstrated 
in  Fig.  23.  Comparison  of  the  corresponding  results  shown  in  Fig.  23  indicates  that  the 
admission  of  continuous  gas  jets  on  both  sides  of  the  liquid  film  causes  a  50  %  reduction 
of  the  sinuous  mode  disturbance  amplitudes  generated  by  the  imposed  modulation  of  v  at 
the  nozzle.  Furthermore,  the  dominant  sinuous  mode  wavelength  is  reduced  by  a  factor  of  2 
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upon  admission  of  equal  velocity  gas  streams;  at  the  same  time,  the  amplitude  of  dilational 
mode  film  disturbances  is  increased.  Using  only  10  %  of  the  energy  admitted  in  the  gas- 
phase  modulated  case  of  Fig.  5,  direct  liquid  film  modulation  produced  film  rupture  at 
approximately  the  same  downstream  location.  However,  as  illustrated  in  Fig.  24,  transverse 
film  deflection  is  significantly  larger  in  this  case. 

-  Direct  liquid-phase  modulation  with  or  without  adjacent  gas  streams  has  also  been  con¬ 
sidered  for  the  pure  dilational  mode  of  film  distortion.  Fig.  25  shows  the  distorting  liquid 
film  for  the  same  parameter  set  as  in  Fig.  3  but  now  with  modulation  of  the  axial  velocity 
component  u  according  to  Eq.  (10)  rather  than  gas-jet  modulation  as  employed  in  Fig.  3. 
As  for  the  previously  discussed  sinuous  case,  the  amount  of  energy  used  to  modulated  the 
liquid  film  in  Fig.  25  was  only  5  %  of  the  energy  added  in  the  case  shown  in  Fig.  3.  In 
contrast  to  the  sinuous  case  of  Fig.  23,  where  film  rupture  has  not  been  observed  within  the 
simulated  time  frame,  the  film  in  Fig.  25  ruptures  early  at  t*  =  27  at  a  downstream  location 
of  approximately  xjh  =  68.  Clearly,  for  the  considered  parameter  configuration,  dilational 
modulation  is  preferred  if  film  rupture  is  achieved  by  direct  liquid-phase  modulation.  On  the 
other  hand,  effective  film  breakup  (characterized  by  small  energy  input  and  short  breakup 
length)  for  gas-phase  modulated  films  is  achieved  by  forcing  sinuous  mode  waves  rather  than 
dilational  ones. 

Comparison  of  Fig.  25  and  26  shows  that  time  for  film  rupture  from  dilational  liquid- 
phase  modulation  is  reduced  by  the  addition  of  surrounding  gas  streams;  however,  film 
breakup  length  is  not  influenced  significantly.  Also,  the  wavelengths  of  dilational  mode 
disturbances  generated  by  direct-liquid  phase  modulation  is  not  significantly  altered  by  the 
admission  of  gas  streams  on  both  sides  of  the  discharging  film. 

From  the  previous  discussion  we  find  that,  liquid  film  disintegration  via  gas-jet  modu¬ 
lation  is  less  effective  than  direct  modulation  of  liquid-phase  flow  variables.  However,  this 
statement  only  refers  to  the  initial  film  rupture  process  discussed  here.  Once  liquid  ligaments 
have  been  detached  from  the  continuous  liquid  film,  further  breakup  of  these  ligaments  can 
only  be  achieved  via  energy  transfer  from  the  surrounding  gas-phase  onto  the  liquid-phase 
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ligaments  or  droplets.  Here,  large  relative  velocities  between  gas-phase  and  discrete  liquid 
phase  is  expected  to  be  essential  to  promote  further  breakup  of  ligaments  and  droplets. 

Extended  Quantitative  Analysis 

The  previous  sections  have  focused  on  the  distortion  and  rupture  characteristics  of  dis¬ 
charging  liquid  films  impacted  by  gas  jets.  Quantitative  data  was  only  provided  for  film 
breakup  lengths  and  times  for  the  different  configurations. 

In  order  to  allow  qualitative  and  quantitative  comparison  of  the  computed  results  with 
future  experimental  observations  of  these  films,  we  present  the  results  from  a  quantitative 
analysis  of  the  time-dependent  variations  in  film  thickness  y  and  film  centerline  location  y 
at  various  positions  downstream  from  the  nozzle.  The  employed  Fourier  analysis  and  the 
consideration  of  phase  diagrams  for  local  values  of  y  and  y  follows  the  work  by  Ramos46  who 
studied  the  drawing  of  annular  liquid  jets  at  low  Reynolds  numbers. 

Figs.  27  through  30  illustrate  phase  plots  obtained  for  the  case  shown  in  Fig.  11  at 
x/h  =  12.5  and  x/h  —  25,  respectively.  Nearly  time-periodic  variations  in  film  thickness 
y(t*)  and  centerline  displacement  y(i* )  occur  close  to  the  nozzle  (Figs.  27  and  28).  The 
centerline  displacement  becomes  increasingly  non-periodic,  further  downstream  (Fig.  30). 
Deviation  from  periodic  behavior  at  x/h  =  25  is  significantly  less  pronounced  for  y(t*)  than 
for  y(t*).  At  x/h  =  12.5  the  phase  of  y(£*)  takes  a  dumbbell-like  or  hourglass-like  shape 
whereas  the  orbit  of  y(t*)  at  the  same  location  is  circular  with  higher-harmonic  contributions 
resulting  in  an  additional  inner  heart-shaped  loop.  Fourier  analysis  of  the  signals  y(t*)  and 
y(t*)  shows  that  the  dumbbell  (hour-glass)  shaped  y-orbit  is  generated  mainly  by  the  first  and 
third  harmonic  of  the  sinuous  modulation  frequency,  i.e.  fm  =  1/(2TP).  The  phase  of  y(f*) 
is  a  result  of  the  first,  second,  third  and  fourth  harmonic  in  the  film  thickness  variations, 
whereby  according  to  the  half-wavelength  thinning  of  nonlinear  sinuous  distorting  planar 
films,  the  value  of  the  first  harmonic  in  the  thickness  variation  is  twice  that  of  the  first 
harmonic  in  the  centerline  variation  of  the  film.  As  the  contribution  of  the  fourth  harmonic 
in  y  decreases  further  downstream,  the  heart-shaped  loop  in  the  y-orbit  changes  to  a  simple 
inner  loop  without  the  “kink”.  As  noted  earlier,  oscillations  in  film  thickness  at  x/h  =  25 
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axe  quasi-periodic  whereas  deviation  from  periodic  behavior  is  already  significant  for  y(t*) 
at  x/h  —  25.  However,  Fourier  analysis  of  the  y  signal  shows  that  the  energy  contained  in 
the  sinuous  film  oscillation  at  that  location  is  still  found  predominantly  in  the  first  and  third 
harmonic  of  the  modulation  frequency  fm.  It  is  noted  here  that  the  Fourier  transform  of  y 
or  y  —  h  and  y  at  a  given  downstream  location  x/h  was  taken  within  a  time  interval  after 
the  film  oscillation  at  that  particular  location  as  been  established  or,  in  other  words,  after 
the  initial  film  disturbance  has  passed  the  considered  downstream  location. 

y  and  y  phases-plots  for  the  cases  illustrated  in  Figs.  6  and  12  are  similar  to  those 
discussed  above.  For  the  reduced  Weber  number  case  of  Fig.  9,  the  inner  loop  within  the 
phase-plot  for  y  at  x/h  =  12.5  is  very  tight  and  without  “kink” .  Fourier  analysis  of  the  time- 
dependent  signal  y  at  that  location  shows  that  here,  the  fourth  harmonic  in  the  thickness 
variation  is  no  longer  important.  For  the  corresponding  phase-plot  of  y,  y  =  0  and  dy/dt  =  0 
are  no  longer  general  symmetry  lines  (as  observed  in  Fig.  28),  however  the  orbit  still  shows 
point-symmetry  around  dy/dt  =  y  =  0. 

It  is  noted  here  that  the  appearance  or  disappearance  of  the  prescribed  “kink”  or  heart- 
shaped  inner  loop  in  the  phase  plot  of  y  does  not  necessarily  constitute  the  presence  of  a 
higher  (in  this  case  fourth)  harmonic  mode.  It  can  also  be  the  result  of  a  variation  in  the 
phase-shift  between  the  lower  frequency  modes. 

In  the  presence  of  a  non-zero  constant  gas  flow  superimposed  onto  the  gas-jet  pulsation 
of  the  base  case  (see  Fig.  18)  the  hour-glass  shaped  orbit  in  the  y  phase-plot  is  still  observed. 
However,  the  amplitudes  in  y  and  dy/dt  increase  with  each  completed  orbit.  Fourier  analysis 
of  the  corresponding  time-dependent  signal  for  y  shows  that  as  for  the  lower  Weber  number 
case  discussed  above,  the  fourth  harmonic  mode  is  no  longer  relevant  here.  The  phase  plot 
of  y  at  x/h  =  12.5  does  not  exhibit  the  heart-shaped  inner  loop. 

For  the  described  gas-jet  modulated  liquid  films,  periodic  behavior  in  the  variation  of  y 
and  y  deteriorates  fast  with  increasing  downstream  distance.  This  is  in  contrast  to  a  liquid 
film  which  is  modulated  within  the  liquid  phase  itself.  For  the  sinusoidally  modulated  film 
in  Fig.  23  without  co-flowing  gas  streams,  variations  in  film  thickness  and  film  centerline 
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location  at  fixed  downstream  location  are  highly  periodic  (after  a  certain  start-up  period)  up 
to  a  downstream  distance  of  x/h  =  50.  See  Figs.  31  and  32  for  x/h  =  12.5.  However,  with 
co-flowing  surrounding  gas-streams  re-occuring  orbits  (even  close  to  the  atomizer  nozzle)  are 
no  longer  observed. 

Even  though  periodic  behavior  in  the  variation  of  y  and  y  deteriorates  with  increasing 
downstream  distance  (and  in  particular  for  gas-phase  modulated  films),  analysis  of  the 
power-spectra  for  the  oscillation  in  film  thickness  and/or  film-centerline  location  shows  that 
in  the  investigated  cases,  the  majority  of  the  film  energy  remains  confined  to  the  harmonic 
modes  present  close  to  the  nozzle,  latter  being  the  first  and  third  harmonic  for  the  film 
centerline  location  and  the  first,  second,  third  and  fourth  harmonic  for  the  film  thickness 
variation.  Accordingly,  and  in  an  effort  to  provide  as  much  significant  quantitative  data  as 
possible,  the  included  tables  list  real  and  imaginary  parts  of  the  Fourier  transforms  for  the 
various  dominant  harmonic  modes  observed  in  the  time-signals  for  the  film  thickness  y  and 
film  centerline  location  y  at  various  downstream  locations. 

Fourier  transforms  X  ( k )  for  length  N  input  sequences  x(n)  have  been  calculated  by  using 
MatLab™  according  to 

N 

X(k)  =  x{n)  exp  [—i2n(k  —  l)[(n  —  1) /AT]]  1  <  k  <  N 

n= 1 

i  N 

x(n)  =  —  X(k )  exp  [i2n(k  —  l)[(n  —  l)/iV]]  1  <  n  <  N 

™  k=l 

where  x(n)  is  given  by  x(n)  =  y(nAt)  or  x(n)  =  y(nAt )  at  a  certain  downstream  location 
and  within  a  certain  time-interval,  i.e.  nsAt  <  nAt  <  neAt,  whereby  At  is  the  time-step 
used  within  the  simulation. 

Since  the  input  sequence  x(n)  is  real,  the  above  equation  for  x(n)  can  be  rewritten  in 
terms  of  a  summation  of  sine  and  cosine  functions  with  real  coefficients 

,  .  1  /7\  /  27t(A:  —  l)(n  —  1)\  .  f  2Tr(k  -  l)(n  -  1)\  ,  . 

x(n )  =  JjT.  a(k ) cos  [— - Y - "J  +  Sm  \  - N - " )  (U) 

where 
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(12) 


a(k)  =  real[X(A:)]  ,  b(k )  =  — imag[X(A:)]  ,  1  <  n  <  N 

With  the  knowledge  of  the  transforms  X (ki)  from  Tables  1  through  4,  the  approximate 
phase  plots  for  y  and  y  at  the  various  downstream  locations  can  be  reconstructed  and 
compared  with  the  results  obtained,  for  example,  from  experimental  observations.  Clearly, 
the  energy  contained  in  modes  other  than  the  ones  described  above  and  provided  in  the 
Tables  is  small  and  justifies  the  omission  of  these  modes  in  the  present  investigation. 

Spray  Angle 

Liquid  film  atomization  and  spray  formation  is  a  three-dimensional  phenomenon  even  for 
planar  films  and  includes  the  formation  of  streamwise  vortical  structures  leading  (depend¬ 
ing  on  operational  conditions)  to  cellular  film  breakup  or  stretched  streamwise  ligament 
breakup.7,8  Besides  the  importance  of  three-dimensional  effects  for  the  spray  formation  pro¬ 
cess,  the  resulting  spray  angle  will  greatly  depend  on  the  conditions  (e.g.  velocities)  at  which 
ligaments  are  shed  continuously  from  the  continuous  liquid  film.  The  present  analysis,  which 
applies  only  up  to  the  point  in  time  where  film  rupture  first  occurs,  cannot  address  these 
issues  (which  are  vital  for  an  accurate  determination  of  the  spray  angle). 

In  their  experiment  on  liquid  films  discharging  from  a  planar  twin-fluid  atomizer,  Man- 
sour  and  Chigier6  considered  the  envelope  of  the  transversely  distorting  liquid  film  as  a 
measure  for  the  spray  angle.  However,  for  modulated  liquid  films  Asare,  Takahashi  and 
Hoffman23  showed  that  following  linear  stability  theory,  the  amplitude  of  the  envelope  grows 
exponentially  up  to  a  point  where  it  saturates  due  to  nonlinear  effects.  The  authors  also 
showed  (for  certain  operational  conditions)  good  agreement  of  their  experimental  results 
with  theoretical  predictions  from  a  simplified  trajectory  theory.  Accordingly,  a  prediction  of 
spray  angle  based  on  the  slope  of  the  envelope  at  the  point  of  film  rupture  greatly  depends 
on  the  film  rupture  point  and  cannot  be  very  reliable. 

Nevertheless,  based  on  the  envelope  growth  observed  for  the  different  configuration, 
indications  are  that  an  increase  in  spray  angle  (with  respect  to  the  base  case  of  Figs.  5  and 
6)  is  expected  with  an  increase  in  Weber  number  (see  Fig.  8)  and  with  an  increase  in  gas- 
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to-liquid  density  ratio  (see  Fig.  16).  Naturally,  the  most  significant  increase  in  spray  angle 
is  expected  for  an  increase  in  vertical  gas-jet  momentum  ratio.  See  Fig.  14,  where  liquid 
ligaments  are  expected  to  be  shed  periodically  from  the  highly  deflecting  film.  Furthermore, 
based  on  the  film  envelope,  a  reduction  in  spray  angle  can  be  expected  for  a  decrease  in 
Weber  number  (see  Fig.  9),  a  reduction  of  vertical  gas-jet  momentum  (see  Fig.  11),  and  a 
reduction  in  gas-jet  impact  angle  (see  Figs.  12  and  13).  For  liquid  films  modulated  at  the 
nozzle  exit,  Fig.  23  indicates  that  the  addition  of  constant  co-flowing  gas  streams  will  result 
in  a  decrease  in  the  spray  angle,  if  indeed  this  angle  is  based  on  the  envelope  behavior  of 
the  distorting  film. 

IV.  SUMMARIZING  REMARKS 

The  flowfield  in  the  vicinity  of  a  twin-fluid  atomizer  has  been  analyzed  numerically. 
Considered  is  a  thin  two-dimensional  inviscid  incompressible  liquid  film  discharging  from 
the  atomizer  centerline  surrounded  on  both  sides  by  gas  jets  which  impact  symmetrically 
onto  the  discharging  film.  Initial  film  distortion  is  enforced  actively  by:  1)  modulation  of 
the  impacting  gas  jets  (indirect  or  gas-phase  modulation)  with  gas  momentum  components 
parallel  and  normal  to  the  liquid  stream,  2)  modulation  in  the  liquid-phase  at  the  atomizer 
exit  (direct  or  liquid-phase  modulation)  or  3)  direct  liquid  film  modulation  with  superpo¬ 
sition  of  continuous  gas  streams.  The  investigation  focused  on  gas-phase  modulated  films. 
Effects  of  different  flow  parameters  on  film  breakup  characteristics  were  studied  by  varying 
one  of  the  flow  parameters  from  a  given  base-case  configuration. 

For  direct  (i.e.,  gas-phase)  modulated  films,  energy  input  at  the  gas-inlet  ports  was  found 
to  be  inadequate  in  identifying  effective  film  rupture  conditions,  the  latter  being  character¬ 
ized  by  short  break-up  length  and  times  at  minimum  energy  input  into  the  system.  For  the 
parameter  domain  surrounding  the  considered  base  case,  film  rupture  was  obtained  more 
effectively  by  sinuous  forcing  and  subsequent  nonlinear  sinuous-dilational  mode  coupling 
rather  than  by  dilations!  forcing  at  the  same  average  energy  flux  into  the  system.  Also,  it 
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was  observed  that  intermediate  gas-jet  pulsing  (without  underlying  continuous  gas  flow)  can 
yield  the  same  or  larger  amplitude  disturbances  than  obtained  through  continuous  distur¬ 
bance  wave  growth  from  continuous  jets  at  the  same  overall  energy  input  into  the  system. 

Direct  modulation  of  liquid  phase  flow-parameters,  i.e.  axial  or  transverse  film  velocity, 
at  the  nozzle  exit  is  more  effective  than  gas-jet  modulation,  the  latter  transferring  only  a 
portion  of  the  modulation  energy  onto  the  film.  Computational  results  also  indicate  that  for 
direct  liquid  phase  forcing,  dilational  modulation  (i.e.  modulation  yielding  dilational  film 
distortion)  is  more  effective  in  causing  film  rupture  than  sinuous  modulations,  while,  as 
noted  earlier,  the  opposite  is  true  for  gas-jet  modulated  liquid  films.  Also,  for  liquid-phase 
modulated  films  admission  of  nonzero  ambient  gas  streams  significantly  reduces  sinuous 
mode  wavenumbers  whereas  influence  on  wavelengths  remains  small  for  dilationally  distort¬ 
ing  films. 

For  gas-phase  modulated  films,  variations  of  Weber  number,  gas  jet  to  liquid  film  mo¬ 
mentum  ratio,  gas-to-liquid  density  ratio  and  pulse  period  were  considered  and  their  effect 
on  film  distortion  was  analyzed.  An  increase  in  Weber  number  from  its  base-case  value 
(i.e.  from  10  to  25)  resulted  in  stretched  band-like  films  and  a  delay  in  film  rupture  due  to 
reduced  sinuous-dilational  mode  coupling.  On  the  other  hand,  reduction  in  Weber  number 
(i.e.  from  10  to  5)  resulted  in  smaller  disturbance  amplitudes  and  larger  break-up  time  and 
length  despite  stronger  nonlinear  mode  coupling  manifested  by  the  contraction  of  the  film 
into  fluid  cylinders  connected  by  thinner  fluid  films.  Transverse  deflection  of  the  film  center- 
line  for  a  five-fold  increase  in  density  ratio  was  similar  to  the  one  observed  for  the  increased 
Weber  number  case;  however,  fluid  blob  formation  was  still  observed. 

A  five-fold  increase  in  vertical  or  parallel  ga.s-jftt.-to-1iqiiid-film-Tnomp.nt.nm  ratio  resulted 
in  “immediate”  film  rupture  caused  by  the  dynamics  of  the  impacting  gas  jets  rather  than  by 
continuous  growth  of  film  disturbances  downstream  with  energy  transfer  from  the  adjacent 
gas-streams. 

Variation  in  pulse  period  affects  film  distortion  in  two  ways:  It  changes  the  amplitude  of 
the  initial  film  disturbances  generated  by  the  pulsed  gas-jets,  since  increasing  the  pulse  period 
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provides  more  time  for  transverse  film  movement  (due  to  one-sided  gas-jet  impact)  before 
being  counter-acted  by  an  opposed  gas  jet.  On  the  other  hand,  it  alters  the  wavelengths 
of  the  disturbances  generated  on  the  film  and  therefore  the  gas-to-liquid  energy  transfer 
downstream  from  Kelvin-Helmholtz-type  wave  growth. 
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Table  1:  Complex  Fourier  coefficients  determined  for  nondimensional  fre¬ 
quencies  /  =  0,0.2, 0.4, 0.6, 0.8  of  y(t*)  and  frequencies  g  =  0,0.1, 0.3  of 
y(t*)  at  various  downstream  locations  x/h.  Data  taken  from  cases  shown  in 
Fig.  6,  Fig.  9  and  Fig.  18.  ( t*nt  denotes  the  sampling  interval.) 


case 

Fig.6 

Fig.9 

Fig.  18 

x/h 

12.5 

12.5 

12.5 

25 

tint 

11  -  20.9 

8  -  17.9 

11  -  31.5 

16.25  -  26.25 

fo 

-1.7425 

0.1243 

-3.6091 

-1.1409 

fo.2 

-2.4001  -i  3.5750 

0.6519 -i  8.5912 

-2.2346 -i  3.2005 

4.9988  +i  2.6593 

foA 

-0.3777  -i  5.7013 

-0.0560  -i  3.4911 

0.2596  -i  2.8377 ' 

-3.3047 -i  1.3497 

fo.6 

-0.3274  -i  1.8957 

-0.9075  +i  0.8650 

0.1585  -i  0.4658 

0.2757 -i  1.2830 

/o.8 

-0.6983  +i  0.4870 

-0.0686  -i  0.1412 

0.0080  +i  0.0422 

0.0847  -i  0.0845 

0o 

9.67 

-4.9318 

46.92 

-63.8480 

00.1 

-126.48  -i  66.94 

81.5998  +i  83.3471 

117.45  +i  170.43 

-11.3874  +i  80.3349 

00.3 

-10.37  -i  27.69 

-11.1639  -i  11.2677 

-16.95  +i  10.15 

0.7101  +i  22.4862 

Table  2:  Complex  Fourier  coefficients  determined  for  nondimensional  fre¬ 
quencies  /  =  0,0.2, 0.4, 0.6, 0.8  of  y{t*)  and  frequencies  g  —  0,0.1, 0.3  of 
y(t*)  at  various  downstream  locations  x/h.  Data  taken  from  cases  shown  in 
Fig.  11  and  Fig.  12.  (t*nt  denotes  the  sampling  interval.) 


case 

Fig.ll 

Fig.  12 

12.5 

25 

12.5 

25 

50 

Emh 

8-18 

11.5  -  21.5 

6-26 

11.5  -  21.5 

17-27 

fo 

-1.1274 

-0.2904 

-2.2123 

-0.3252 

-0.1413 

/o.2 

2.1313  +i  0.1244 

5.2311  +i  4.6508 

-3.0126  -i  2.7513 

4.5862  +i  4.4771 

8.0436  +i  12.3251 

foA 

-2.6680  -i  0.2732 

-4.4851  -i  5.4008 

-0.4646  -i  4.7289 

-3.4008  -i  5.0893 

-15.8212  +i  7.7273 

/o.6 

0.7068  -i  0.5327 

1.9118  -i  2.8906 

-0.1792  -i  1.5370 

1.9544  -i  1.9449 

-3.4856  +i  4.5756 

fo.8 

0.3036  -i  0.2878 

0.6632  +i  0.0905 

-0.2280  +i  0.5580 

0.4152  +i  0.2162 

-0.0874  +i  0.6210 

9o 

-5.3221 

-16.7992 

4.68 

-17.8351 

5.01 

9o.i 

-18.4821  +i  75.7047 

-32.4973  +i  75.9395 

133.87  +i  93.26 

-40.7346  +i  78.6269 

-29.20  -i  100.29 

90.3 

5.1482  -i  14.6942 

-14.3083  +i  30.8077 

3.72  +i  30.94 

-15.5162  +i  25.3493 

7.14  -i  32.95 
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Table  3:  Complex  Fourier  coefficients  determined  for  nondimensional  fre¬ 
quencies  /  =  0,0.2, 0.4, 0.6, 0.8  of  y(t*)  and  frequencies  g  =  0,0.1, 0.3  of 
y(i*)  at  various  downstream  locations  x/h.  Data  taken  from  cases  shown  in 
Fig.  23.  (i*nt  denotes  the  sampling  interval.) 


case 

Fig.23  (dash, solid) 

Fig.23  (dot, dash-dot) 

x/h 

12.5 

25 

50 

12.5 

tint 

8-38 

10  -  40 

20-40 

7-27 

fo 

1.0959 

2.6197 

0.8649 

-3.7106 

/o.2 

5.0258  -i  0.3302 

-12.3063  -i  15.3748 

4.1233  +i  44.0686 

-0.8199  -i  3.9278 

foA 

-1.2484  +i  5.7662 

-1.0425  -i  16.3421 

-26.9334  +i  3.7302 

2.4076  -i  3.7397 

/o.6 

-0.4002 -i  0.6104 

1.1145  -i  5.6337 

-5.6068  -i  4.2925 

0.8231  +i  0.2940 

/o.8 

1.5359  +i  1.1584 

0.1498  +i  0.2117 

0.2417  +i  1.1853 

-0.2590  +i  1.4226 

00 

0.21 

-12.7 

1.93 

63.36 

00.1 

12.85  +i  536.95 

580.55  +i  657.35 

-391.44  -i  31.04 

226.70  +i  231.71 

00.3 

-9.86  -i  55.09 

25.25  -i  33.98 

-29.35  +i  12.68 

-10.70  +i  19.00 

Table  4:  Complex  Fourier  coefficients  determined  for  nondimensional  fre¬ 
quencies  /  =  0,0.2, 0.4, 0.6, 0.8  of  y(t*)  at  various  downstream  locations 
x/h.  Data  taken  from  cases  shown  in  Fig.  3  and  Fig.  4.  (tint  denotes  the 
sampling  interval.) 


case 

Fig.3 

Fig.4 

12.5 

25 

E&HHHRHl 

12.5 

25 

50 

t’irU 

7-37 

10  -  40 

15.5  -  45.5 

5-25 

9.5  -  24.5 

15-25 

fo 

-1.0487 

-1.1002 

-1.0912 

-1.4044 

-1.0758 

-0.7213 

/o.l 

-0.6770  -  i  1.0303 

-1.8518  -  i  1.9126 

-4.6015  -i  3.4093 

-0.0004  +i  0.0003 

-0.0008-i  0.00115 

0.0045  +i  0.0009 

fo.2 

0.8751  -  i  0.9734 

2.4318  -  i  2.0079 

4.3950  -i  4.8127 

-1.6982  +i  0.3701 

0.8970  -i  2.9529 

1.0881  -i  4.6293 

fo.3 

0.8444  +  i  0.4398 

1.4058  +  i  1.2092 

2.7260+i  0.9146 

0.0026  -i  0.0002 

0.0015  -i  0.00285 

-0.0016  +i  0.0095 

foA 

0.0424  +  i  0.3450 

-0.0645  +  i  0.1913 

-0.0625  -i  0.4539 

-0.4588  +i  0.1487 

0.4805+i  0.2696 

0.2182  +i  0.0205 

/o.5 

0.0512  -  i  0.0836 

0.0737  -  i  0.2080 

0.1080  -i  0.4961 

-0.0027  -i  0.0049 

-0.0041  -i  0.00185 

0.0136  +i  0.0303 

/o.6 

0.0322  +  i  0.0066 

0.0178  +  i  0.1015 

0.1618  -i  0.0052 

0.0465  -i  0.0973 

0.0743  -i  0.1839 

-0.4138  -i  0.3860 
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Figure  captions 

Figure  1:  Schematic  of  Twin-Fluid  Atomizer  Configuration 

Figure  2:  Solution  procedure  for  dilational  mode  film  forcing  (y  =  0,  v  =  0  ;  h  =  y,  u  =  u, 
(j)  =  (j)g.ll  pg  =  pg.1}  r  =  Ti,  y_  =  —  y/2,  pc  =  capillary  pressure  along  lower  interface). 

Figure  3:  Dilationally  distorting  gas-jet-modulated  semi-infinite  planar  film  under  base-case  forc¬ 
ing  conditions  according  to  Eq.  (1)  and  for  base-case  parameter  set:  sjh  —  H/h  =  10, 
pg/ pi  =  0.0012,  Wej  =  10,  Mx  =  My  =  1.  Tp  =  5,  n6] i,2  =  0. 

Figure  4:  Dilationally  distorting  gas-jet-modulated  film  with  base-case  parameter  configuration 
and  modified  pulse  characteristics. 

Figure  5:  Sinusoidally  distorting  gas-jet-modulated  film  under  base-case  forcing  conditions  ac¬ 
cording  to  Eq.  (1)  and  for  base-case  parameter  set:  s/h  =  H/h  —  10,  pg/ pi  =  0.0012, 
Wej  =  10,  Mx  =  My  =  1.  Tp  =  5,  ng-i  =  0,  ng- 1  =  1. 

Figure  6:  Discharging  sinusoidally  distorting  gas-jet-modulated  liquid  film  according  to  Fig.  5 
with  modified  start-up  conditions. 

Figure  T:  Discharging  gas-jet-modulated  liquid  film  according  to  Fig.  5  under  mixed-mode  forc¬ 
ing  conditions  (ng. i  =  0,  ne;2  =  0.5). 

Figure  8:  Effect  of  Weber  number  increase  (Wej  =  25)  on  discharging  sinusoidally  distorting 
gas-jet-modulated  liquid  film. 

Figure  9:  Effect  of  Weber  number  decrease  (Wej  =  5)  on  discharging  sinusoidally  distorting 
gas-jet-modulated  liquid  film. 

Figure  10:  Effect  of  Weber  number  increase  (Wej  =  25)  on  discharging  dilationally  distorting 
gas-jet-modulated  liquid  film. 

Figure  11:  Effect  of  decreased  vertical  gas-jet  momentum  (Mx  =  0.5  on  discharging  sinusoidally 
distorting  gas-jet-modulated  liquid  film.  M,|  =  1). 
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Figure  12:  Effect  of  decreased  gas-jet  injection  angle  a  (a  =  30°,  M  =  4-  M\  =  \/2)  on 

discharging  sinusoidally  distorting  gas-jet-modulated  liquid  film. 

Figure  13:  Effect  of  decreased  gas-jet  injection  angle  (a  =  15°)  on  discharging  sinusoidally  dis¬ 
torting  gas-jet-modulated  liquid  film. 

k 

Figure  14:  Effect  of  increased  vertical  gas-jet  momentum  (Mj_  =  5)  on  discharging  sinusoidally 
distorting  gas-jet-modulated  liquid  film. 

Figure  15:  Effect  of  increased  parallel  gas-jet  momentum  (Mj|  =  5)  on  discharging  sinusoidally 
distorting  gas-jet-modulated  liquid  film. 

Figure  16:  Effect  of  increased  gas-to-liquid  density  ratio  (pg/pt  =  0.006)  on  discharging  sinu¬ 
soidally  distorting  gas-jet-modulated  liquid  film. 

Figure  IT:  Effect  of  decreased  pulse  period  (Tp  =  2.5)  on  discharging  sinusoidally  distorting  gas- 
jet-modulated  liquid  film. 

Figure  18:  Effect  of  non-zero  time-independent  gas-jet  velocity  on  discharging  sinusoidally  dis¬ 
torting  gas-jet-modulated  liquid  film  (Tp  =  5;  u'g.Q  =  0.5tt3;o,ttg;3(i)  =  0.25). 

Figure  19:  Effect  of  increased  pulse  period  (Tp  =  8.75)  on  on  discharging  sinusoidally  distorting 
gas-jet-modulated  liquid  film  with  underlying  non-zero  gas-jet  velocity  ug.s. 

Figure  20:  Effect  of  increased  pulse  period  (Tp  =  20)  on  discharging  sinusoidally  distorting  gas- 
jet-modulated  liquid  film  with  underlying  non-zero  gas-jet  velocity  ug;s. 

Figure  21:  Dimensionless  growth  rate  Im[cj/(&£/o)]  as  function  of  kh  =  2-jrh/X  for  planar  liquid 
film  with  co-flowing  gas  streams  (at  relative  velocity  Uq)  and  We  =  piU^h/cr  =  971 
and  pg/pi  =  0.0012. 

Figure  22:  Influence  of  domain  length  (xmax/h  =  250)  on  discharging  sinusoidally  distorting  gas- 
jet-modulated  liquid  film  with  underlying  non-zero  gas-jet  velocity  ug.s. 
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Figure  23:  Sinusoidally  distorting  liquid-phase-modulated  film  according  to  Eq.  (10),  with  or  with¬ 
out  adjacent  constant  velocity  gas  streams  and  under  base-case  parameter  conditions: 
s/h  =  H/h  =  10,  pg/pi  =  0.0012,  We,  =  10,  Mx  =  Mu  =  1,  Tp  =  5.  ug,3  =  0:  t*  =  9 
(dashed),  19  (solid);  ug.3  =  0.45:  t*  =  9  (dotted),  19  (dash-dot).  Without  continuous 
gas  jets,  energy  input  is  5%  of  that  in  Fig.  5. 

Figure  24:  Sinusoidally  distorting  liquid-phase-modulated  film  according  to  Eq.  (10),  without 
adjacent  constant  velocity  gas  streams  and  under  base-case  parameter  conditions: 
s/h  =  H/h  =  10,  pg/pi  =  0.0012,  Wei  =  10,  Mx  =  M,|  =  1,  Tp  =  5.  Energy 
input  is  twice  that  of  Fig.  23. 

Figure  25:  Dilationally  distorting  liquid-phase-modulated  film  according  to  Eq.  (10)  injected  into 
a  quiescent  ambient  gas  under  base-case  parameter  conditions:  s/h  =  H/h  =  10, 
pg/pl  =  0.0012,  Wei  =  10,  Mj.  =  Mj|  =  1,  Tp  =  5.  Energy  input  is  5%  of  that  in  Fig.  5. 

Figure  26:  Dilationally  distorting  liquid-phase-modulated  film  according  to  Eq.  (10)  with  adjacent 
constant  velocity  gas  streams  and  under  base-case  parameter  conditions:  s/h  =  H/h  = 
10,  pg/pi  —  0.0012,  Wej  =  10,  Mjl  =  My  =  1,  Tp  =  5.  Energy  input  from  direct  film 
modulation  equals  that  in  Fig.  25. 

Figure  27:  Phase  plot  for  y(t*)  at  x/h  =  12.5  according  to  the  case  shown  in  Fig.  11.  (Nondi- 
mensional  time  t*  progresses  in  the  clockwise  direction.) 

Figure  28:  Phase  plot  for  y(t*)  at  x/h  =  12.5  according  to  the  case  shown  in  Fig.  11. 

Figure  29:  Phase  plot  for  y(t *)  at  x/h  —  25  according  to  the  case  shown  in  Fig.  11. 

Figure  30:  Phase  plot  for  y(t*)  at  x/h  =  25  according  to  the  case  shown  in  Fig.  11. 

Figure  31:  Phase  plot  for  y(t*)  at  x/h  =  12.5  according  to  the  case  shown  in  Fig.  23  without 

co-flowing  gas  streams. 
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Figure  32:  Phase  plot  for  y(t*)  at  x/h  =  12.5  according  to  the  case  shown  in  Fig.  23  without 
co-flowing  gas  streams. 
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FIG.  1.  Mehring,  Phys.  Fluids 
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FIG.  2.  Mehring,  Phys.  Fluids 


51 


FIG.  3.  Mehring,  Phys.  Fluids 


FIG.  4.  Mehring,  Phys.  Fluids 


FIG.  5.  Mehring,  Phys.  Fluids 


54 


x/h 


FIG.  6.  Mehring,  Phys.  Fluids 
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FIG.  7.  Mehring,  Phys.  Fluids 
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FIG.  8.  Mehring,  Phys.  Fluids 
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FIG.  9.  Mehring,  Phys.  Fluids 
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FIG.  10.  Mehring,  Phys.  Fluids 
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FIG.  11.  Mehxing,  Phys.  Fluids 
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FIG.  12.  Mehring,  Phys.  Fluids 
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FIG.  13.  Mehring,  Phys.  Fluids 
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FIG.  14.  Mehring,  Phys.  Fluids 
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FIG.  15.  Mehring,  Phys.  Fluids 
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FIG.  16.  Mehring,  Phys.  Fluids 
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FIG.  17.  Mehring,  Phys.  Fluids 
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FIG.  18.  Mehring,  Phys.  Fluids 
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FIG.  19.  Mehring,  Phys.  Fluids 
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FIG.  20.  Mehring,  Phys.  Fluids 
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FIG.  21.  Mehring,  Phys.  Fluids 
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FIG.  22.  Mehring,  Phys.  Fluids 
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FIG.  23.  Mehring,  Phys.  Fluids 
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FIG.  24.  Mehring,  Phys.  Fluids 
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FIG.  25.  Mehring,  Phys.  Fluids 
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FIG.  26.  Mehring,  Phys.  Fluids 
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FIG.  28.  Mehring,  Phys.  Fluids 


77 


3(y/h)/3t 


0.08  - 

0.04  - 

o  - 

-0.04  - 

-0.08  j 

0.95 


I  I  I - 1 - 1 - 1 _ 1—1 _ I _ L-.l  ,  .1 _ I _ I _ I _ L 

0.975  1  1.025 

y/h 


J _ I _ I 

1.05 


FIG.  31.  Mehring,  Phys.  Fluids 
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Abstract 

The  linear  and  nonlinear  distortion  and  breakup  of  three-dimensional  swirling  or 
non-swirling  annular  and  swirling  conical  thin  inviscid  liquid  films  are  analyzed  by 
means  of  a  reduced-dimension  approach.  The  films  discharge  from  an  annular  slit 
nozzle  or  atomizer  into  a  gas  of  negligible  density  at  negligible  gravity  conditions. 
Nonlinear  numerical  simulations  describe  film  distortion  up  to  the  time  when  film  rup¬ 
ture  first  occurs.  Linear  and  nonlinear  solutions  are  presented  and  discussed  for  various 
configurations  and  with  either  dilational  or  sinuous  three-dimensional  modulations  im¬ 
posed  onto  the  films  at  the  nozzle  exit.  Nonlinear  growth  rates  can  be  significantly 
larger  than  predicted  by  linear  theory.  Initially  axisymmetric  disturbances  remain  ax- 
isymmetric  and  fluctuations  in  the  circumferential  direction  generated  by  numerical 
error  are  not  amplified  for  the  considered  cases.  Overall  film  topology  at  the  time 
of  film  rupture  suggests  that  single  dilational  or  sinuous  oblique  waves  will  result  in 


1 


f 


i 


spiraling  filaments  detaching  from  the  continuous  film.  Combination  of  clockwise  and 
counter-clockwise  travelling  dilational  waves  results  in  an  approximately  rectangular 
array  of  larger  ligaments  interspaced  by  thin  fluid  films  indicating  subsequent  cellular- 
type  breakup  for  both  annular  and  conical  configurations.  Results  for  similar  sinuous 
mode  modulations  suggest  film  disintegration  via  shedding  of  continuous  annular  rings 
from  swirling  annular  films  and  filament  shedding  from  swirling  conical  films. 

Introduction 

The  stability  of  liquid  sheets  plays  an  important  role  in  various  technological  applications 
including  spray  combustion  of  liquid  fuels  in  furnaces,  internal  combustion  piston  engines, 
jet  engine  combustion  chambers  and  rocket  motors  [?].  The  combustion  of  liquid  fuels  in  the 
described  applications  is  frequently  achieved  through  the  generation  and  disintegration  of 
swirling  conical  liquid  films  [2].  The  stability  of  annular  films  is  of  relevance  in  applications 
such  as  film  shell  encapsulation  or  shell  formation,  film  blowing,  as  well  as  in  the  use  of 
collapsing  annular  films  as  chemical  reactors  for  toxic  waste  incineration. 

Experimental  observations  of  disintegrating  conical  liquid  sheets  have  already  been  de¬ 
scribed  by  various  authors.  See  Ref.  [3]  for  a  more  detailed  review.  Theoretical  analyses  of  the 
phenomenon  of  liquid-film  breakup  were  in  general  limited  to  linear  planar  or  annular  sheet 
configurations;  the  latter  with  or  without  swirl  [4],  In  particular,  linear  three-dimensional 
analyses  of  planar  or  swirling  annular  sheets  were  only  presented  bv  Ibrahim  and  Akpan  f5l . 
Ponstein  [6]  and  Panchagnula.  Sojka  and  Santangelo  [4].  A  linear  analysis  of  swirling  conical 
sheets  poses  a  variable  coefficient  problem,  and  has  not.  to  the  knowledge  of  the  authors. 
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been  presented  in  the  literature.  The  few  nonlinear  analyses  presented  on  the  subject  were 
also  limited  to  the  planar  geometry  [7,  8,  9,  10,  11],  annular  sheets,  or  liquid  bells  without 
swirl  [12,  13,  14,  15].  A  nonlinear  analysis  of  truly  three-dimensionally  distorting  (swirling 
or  non-swirling)  annular  and  swirling  conical  films  discharging  from  a  nozzle  or  atomizer  and 
subject  to  dilational  or  sinuous  modulations  at  the  nozzle  exit  has  not  been  presented  so 
far,  even  though  the  prescribed  configurations  (in  particular  swirling  conical  liquid  films)  are 
of  importance  in  practical  atomization  systems,  such  as  swirl-cup  atomizers  or  pre-filming 
atomizers. 

Linear  and  nonlinear  analyses  of  discharging,  modulated,  pressure-stabilized,  thin  annu¬ 
lar  liquid  films  were  presented  in  Ref.  [17].  The  linear  analysis  predicted  that,  in  general,  film 
modulation  generates  four  dilational  and  two  sinuous  mode  wavenumbers.  Group  velocity 
analysis  (with  terms  of  0(R~Z)  neglected)  and  consideration  of  the  Sommerfeld  radiation 
condition  downstream  lead  to  the  conclusion  that  only  two  dilational  mode  wavenumbers 
will  appear  on  the  discharging  film  (k2  and  k2  in  Ref.  [17]).  For  sufficiently  large  forcing 
frequencies  u  or  Weber  numbers  We  (based  on  liquid  density,  undisturbed  film  thickness 
and  undisturbed  liquid  velocity  at  the  nozzle  exit),  wavenumber  k3  resulted  in  an  exponen¬ 
tially  growing  dilational  mode  wave  whereas  k2  (its  complex  conjugate  value)  resulted  in 
an  exponentially  decaying  disturbance  downstream.  Linearly  unstable  sinuous  mode  waves 
are  predicted  only  for  Weber  numbers  We  >  2  at  forcing  frequencies  u>  <  (1  IR)Jl  -  2 /We, 
i.e.  where  wavenumber  l2  resulted  in  an  exponentially  growing  wave  and  l\  (the  complex 
conjugate  of  l2)  produced  an  exponentially  decaying  disturbance  downstream.  For  small 
Weber  numbers,  i.e.  We  <  5,  nonlinear  and  linear  numerical  simulations  of  the  initial  and 
boundary-value  problem  were  dominated  by  transient  effects.  In  particular,  film  collapse  and 
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bubble  formation  was  observed  to  occur  near  the  wavefront  as  the  initial  film  disturbance 
propagated  into  the  undisturbed  film.  In  Ref.  [3]  axisymmetric,  swirl-stabilized  annular  film 
were  analyzed.  According  to  linear  analysis  and  in  contrast  to  the  pressure-stabilized  case, 
no  pure  dilational  mode  is  found  on  swirl-stabilized  annular  films.  For  swirl-stabilized  films, 
_  collapse  of  the  undisturbed  annular  configuration  is  prevented  by  the  centrifugal  forces  result¬ 
ing  from  the  swirling  motion;  for  pressure-stabilized  (non-swirling)  annular  film-;  a  constant 
positive  pressure  difference  is  maintained  between  the  inner  core  and  the  surrounding  of 
the  annulus  in  order  to  prevent  the  film  from  collapsing  on  itself  due  to  surface  tension. 
In  the  pressure-stabilized  case,  linear  mode-coupling  occured  solely  from  the  sinuous  to  the 
dilational  mode  via  the  conservation-of-mass  equation.  Nevertheless,  for  large  enough  We¬ 
ber  numbers  and  low  forcing  frequencies,  linear  dilational  or  predominantly  dilational  mode 
waves  behave  qualitatively  the  same  on  both  pressure-  and  swirl-stabilized  films.  However, 
the  instability  range  and  unstable  growth  rate  for  linear  "dilational”  waves  is  significantly 
larger  for  swirl-stabilized  films  in  comparison  to  pressure-stabilized  ones,  i.e.  u>^3  >  ui^p  with 
ui  being  the  forcing  frequency  at  the  stability  limit  and  subscripts  s  and  p  denoting  swirl-  and 
pressure-stabilized  films,  respectively.  A  significantly  larger  growth  rate  of  "dilational”  mode 
waves  on  swirling  axisymmetric  annular  films  has  also  been  observed  in  the  comparison  of  a 
dilationally  modulated  nonlinear  swirl-stabilized  film  and  a  similar  pressure-stabilized  films 
[3].  Qualitative  differences  between  linear  "dilational”  mode  waves  on  swirl-  and  pressure- 
stabilized  annular  films  were  found  at  increased  forcing  frequencies,  e.g.  u >  >  0.6  (for  large 
Weber  numbers),  and  even  at  lower  forcing  frequencies  if  the  Weber  numbers  were  small 
enough.  However,  in  the  parameter  and  comparative  numerical  studies  of  Ref.  [3],  cases 
with  large  modulation  frequencies  were  not  considered  due  to  the  limitations  of  the  em- 
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ployed  thin  film  model.  Small  Weber  number  cases  were  also  excluded  from  the  analysis, 
due  to  the  dominance  of  transient  effects  in  the  corresponding  numerical  simulations.  In 
strong  contrast  to  pressure-stabilized  films,  where  linear  theory  predicts  that  unstable  sinu¬ 
ous  waves  are  only  found  for  Weber  numbers  We  >  2,  predominantly  sinuous-mode  unstable 
waves  on  swirl-stabilized  films  are  only  observed  for  We  <  2.  However,  as  in  Ref.  [17]  and 
due  to  the  dominant  behavior  of  transient  effects  within  the  numerical  analysis,  swirling 
annular  (or  conical)  films  at  small  Weber  numbers  were  not  studied  in  Ref.  [3].  Rather,  the 
parameter  studies  of  Ref.  [3]  (analyzing  the  effect  of  Weber  number,  modulation  amplitude 
and  swirl  number  on  breakup  time  and  breakup  length)  were  conducted  around  four  annular 
(e.g.  swirl  number  k  =  1)  base  cases  with  modulations  of  the  axial  or  transverse  film  ve¬ 
locity  at  the  nozzle  exit  at  We  =  100  and  forcing  frequencies  ui  <  0.42.  The  linear  analysis 
of  Ref.  [3]  was  constrained  to  swirling  axisymmetric  annular  films.  Initial  conditions  for  the 
numerical  simulations  of  axisymmetric  swirling  conical  films  in  Ref.  [3]  were  obtained  by  inte¬ 
grating  the  nonlinear  steady-state  equations  from  the  nozzle  exit  with  appropriate  boundary 
conditions  at  the  nozzle  and  a  4th-order  Runge-Kutta  integrator.  Modulation  of  the  axial 
and/or  transverse  film  velocity  were  of  the  form  v\  T  =  AZ>T  [l  -  e~tm/T‘ ]  sin  (27rt*/Tp)  with 
the  nondimensional  time-period  of  the  harmonic  forcing  Tv  and  the  nondimensional  forcing 
amplitudes  of  the  axial  or  transverse  velocity  modulation  Az<r. 

A  brief  discussion  on  three-dimensionally  modulated  swirling  annular  and  conical  films 
was  included  in  Ref.  [3].  However,  this  discussion  was  limited  to  the  numerical  solution 
of  initially  axisymmetric  film  modulations  or  axisymmetric  modulations  with  superimposed 
small-amplitude  circumferential  disturbances.  Three-dimensional  solutions,  for  sheets  with 
axisymmetric  modulations  only,  remained  axisymmetric  throughout  the  simulations,  indi- 
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eating  that  small  perturbations  due  to  numerical  error  do  not  result  in  a  three-dimensional 
capillary  instability.  For  the  considered  annular  and  conical  films  with  harmonic  forcing  of 
the  axial  or  transverse  velocity,  the  characteristics  of  the  film  distortion  remained  predomi¬ 
nantly  axisymmetric. 

The  present  analysis  extends  the  three-dimensional  discussion  of  Ref.  [3]  by  considering 
film  modulations  which  result  in  truly  three-dimensional  film  distortion  and  rupture,  as 
well  as  by  incorporating  a  detailed,  description  of  the  three-dimensional  model.  In  addition, 
the  present  analysis  includes  a  three-dimensional  linear  spatial  analysis  of  swirling  annular 
films.  In  analogy  to  the  axisymmetric  analysis  of  Ref.  [3],  the  linear  analysis  1)  identifies  the 
appropriate  number  of  boundary  conditions  for  the  numerical  simulations  and  2)  identifies 
the  magnitude  and  importance  of  nonlinear  effects  by  comparison  between  linear  analytical 
and  nonlinear  numerical  results.  Analogous  to  previous  work  by  the  authors  [3,  7,  17,  18, 
19],  group  velocities  of  the  capillary  waves  emanating  from  the  modulated  nozzle  and  from 
downstream  are  obtained  from  the  linear  analysis  and  are  used  to  identify  the  number  of 
boundary  conditions  to  be  specified  at  the  atomizer  exit  and  the  number  at  downstream 
infinity  for  both  linear  and  nonlinear  analyses. 

In  the  analysis  of  thin  liquid  sheets  or  films  presented  here,  only  capillary  and  inertia 
effects  are  considered.  Subsequently,  the  term  ‘annular’  refers  to  thin-walled  cylindrical 
liquid  columns  with  approximately  ring-shaped  cross-sectional  area  of  constant  time-averaged 
radius  along  the  axial  direction.  Sheets  or  films  which  roughly  resemble  hollow  cones,  i.e., 
thin-walled  liquid  columns  with  monotonically  increasing  or  decreasing  annular  radius  of 
the  ring-shaped  cross-sectional  area,  are  named  ‘conical’  sheets.  Clearly,  sheets  or  films 
discharging  from  an  annular  nozzle  and  with  a  monotonically  increasing  or  decreasing  annular 
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radius  within  the  region  near  the  atomizer  are  still  described  as  being  ‘conical’  even  if  the 
annular  radius  undergoes  oscillations  further  downstream.  The  latter  phenomenon  might 
occur  due  to  the  dynamic  exchange  of  translational  or  rotational  kinetic  energy  and  potential 
or  surface  energy. 

Problem  Formulation 

The  present  work  extends  previous  work  by  the  authors  on  non-swirling  annular  and  swirling 
conical  sheets.  [3,  16,  17,  18,  21].  Considered  are  semi-infinitely-long  thin  swirling  (or  non¬ 
swirling)  liquid  films  exiting  from  an  annular  nozzle  or  atomizer  as  shown  in  Fig.  1. 

Three-dimensional  dilational  and  sinuous  disturbances  are  considered.  Liquid  viscosity 
is  neglected  and  the  sheet  is  exiting  into  a  void  under  negligible  gravity.  The  assumption  of 
thin  sheets  allows  the  reduction  of  the  dimensionality  of  the  problem  by  integrating  across 
the  thickness  of  the  sheet.  This  approach,  also  referred  to  as  the  lubrication  approximation, 
has  been  employed  by  Ramos  [14]  and  by  the  authors  [3,  7,  16,  17,  20]  for  the  analyses  of 
thin  planar  and  axisymmetric  annular  sheets  without  swirl,  as  well  as  axisymmetric  swirling 
conical  sheets. 

Only  spatial  film  stability  or  film  distortion  will  be  of  interest  here.  The  spatially  periodic 
temporal  film  distortion  and  film  stability  is  not  relevant  to  the  conical  sheet  with  its  varying 
radius.  The  temporal  instability  can  apply  to  the  annular  cylindrical  sheet  and  for  the 
radially  expanding  sheet. 
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Governing  Equations 


The  governing  equations,  describing  the  unsteady  motion  in  an  incompressible,  inviscid 
three-dimensionally  distorting  liquid  sheet  under  zero  gravity  and  in  a  cylindrical  coordinate 
system,  are  given  by 


dvz  dvr  vr  1  dvg 
dz  dr  +  r  r  86 


dvz  dvz  dvz  1  dvz  1  dp 

W+v‘-aI+v'-dF  +  rVlW  =  -pa; 


dvr  dvT  dvr  1  dvT  vj  1  dp 

dt  JrVzdz  Vt  dr  rVd  d9  r  pdr 


dvg 

~St 


dve  1  dvg  VT  Vg 

+  vr-£-  +  -U9-5T  + - 

dr  r  do  r 


1  dp 
pr  d6 


(1) 

(2) 

(3) 

(4) 


where  vz,  vT  and  vg  are  the  velocity  components  in  the  axial  (z-),  radial  (r-)  and  circumfer¬ 
ential  (0-)  directions,  respectively.  See  Fig.  1.  p  and  pi  denote  the  pressure  and  the  density 
of  the  liquid. 

Indicating  the  outer  and  inner  location  of  the  sheet  by  r+(z,t)  and  r_(z,  t),  we  define 
the  radial  centerline  position  and  the  thickness  of  the  sheet  by  f(z,t )  =  (r+  +  r_)/2  and 
A r(z,t)  =  r+  —  r_.  The  pressure  and  the  radial  velocity  component  at  the  fluid  interfaces 
(p±,  vr,±)  are  given  by  the  following  kinematic  and  dynamic  boundary  conditions, 


dr±  dr± 

Vt±  =  —  +  V:,±  — 


1 

-Vg 


dr± 

~5e 


15) 
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P±  =  Pg,±  ±  CTK±  (g) 

where  the  local  curvature  k±  is  given  by  the  divergence  of  the  unit  normal  vector  at  the 
particular  location  of  the  outer  (+)  or  inner  (— )  sheet  interface,  i.e.  k±  =  V  •  n±.  In  the 
previous  equation  pg  is  the  pressure  of  the  surrounding  gas  and  cr  denotes  the  surface  tension 
coefficient  of  the  liquid.  For  non-swirling  annular  sheets  pg>+  =  0  but  pg_  ^  0  in  order  to 
stabilize  the  annular  sheet  in  its  undisturbed  configuration.  The  latter  is  also  referred  to  as 
‘pressure  stabilization’.  See  Ref.  [17]  in  this  context.  This  contrasts  with  ‘swirl-stabilized’ 
annular  or  conical  sheets,  where  p3i±  =  0  and  stabilization  of  the  undisturbed  (unmodulated) 
sheet  is  obtained  by  balancing  surface  tension  forces  with  the  centrifugal  forces  due  to  swirl. 
Clearly,  gas-phase  effects  are  not  included  in  the  present  analysis.  The  pressure  in  the  gas 
surrounding  the  considered  liquid  films  is  assumed  constant  because  at  atmospheric  pressure, 
the  gas  density  is  negligible  compared  to  the  liquid  density.  The  influence  of  a  gas-phase 
flow  field  surrounding  discharging  planar  liquid  films  has  been  analyzed  by  the  authors  for 
the  case  of  a  two-dimensional  twin-fluid  atomizer  [20]. 

Expanding  the  dependent  variables  in  terms  of  r  -  f(x,  t)  and  employing  a  similar  expan¬ 
sion  for  1/r,  one  obtains  the  following  system  of  equations  for  r,  A r  and  the  leading-order 
expansion  terms  of  the  velocity  components  v2,vT,vg,  i.e.  u0,vQ  and  w0,  respectively. 


dAr  djupAr)  A  r  1  d(A  rw0) 
dt  dz  +  r  V°  '  r  dO 


=  0 


(") 
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(8) 


dr 


dr 


1  df 


v°"ai+'‘l‘di+w°T'ae 


du0  du0  1  du0 
_+Uo_ 


ITS?  .  (—?L  1  dAr\ 

p  dz  P  \  Ar  dz  +  4r  dz  ) 


(9) 


dvp  dvo_  w0  dv0  Wq  _  1 A p 
dt  U°  dz  **"  r  dO  f  p  A r 


(10) 


dw0  dw0  wo  dwo  w0  v0  1  f  dp  A p  dr 

— - h  Un— - 1 - 1 - = - — - - - 

dt  dz  r  dd  r  pr  dd  A r  d6 


(11) 


with 


_  s  p±±t  =  ^[v.ff+_v.s_] 

A p  =  p+  —  p_  =  cr  [V  •  n+  +  V  •  n_] 


whereby  the  divergence  of  the  normal  vector  at  the  outer  and  inner  interfaces  at  a  specific 
location  (0O,  zq)  is  evaluated  from 


V-n  =  (l  -  ^reg  -r]  -  rr„)  \r)  +  fl]  1/2 
+  ^  [r*  +  3/2  2r*rM  +  +  2 r2rzrzB 

+  ^  [r?  +  Q]"3/2  [2 +  2r2r:r J 


(12) 


(13) 
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whereby  ft  =  r2  (1  +  r2)  and  r  =  r  ±  Ar/2. 

Equations  (7)  through  (11)  form  a  closed  system  of  partial  differential  equations  which 
together  with  appropriate  boundary  and  initial  conditions  govern  the  three-dimensional  non¬ 
linear  distortion  of  thin  swirling  liquid  films  exiting  from  a  nozzle  or  atomizer  into  a  void. 
For  the  subsequent  analyses,  the  above  equations  have  been  nondimensionalized  by  using 
the  undisturbed  film  thickness  and  film  discharge  velocity  at  the  nozzle  exit,  i.e.  Ar0  and 
u7  0,  as  characteristic  length  and  velocity,  respectively.  The  Weber  number  in  the  resulting 
nondimensional  equations  is  given  by  We  =  piv^0Ar0/cr. 

Governing  equations  for  the  nondimensional  linearized  problem  have  been  obtained  by 
introducing 

A r*  =  1  +  h,  r*  =  R  +  77,  uj  =  1  +  u*,  ujj  =  V  +  v*,  Wq  =  W  +  w* 

into  the  nondimensional  form  of  Eqns.  (7)  through  (11).  In  the  above  equations,  h  and  77  de¬ 
note  nondimensional  fluctuations  in  the  film  thickness  and  transverse  film  centerline  location, 
respectively,  u *,  v*  and  w”  represent  nondimensional  fluctuations  in  the  axial,  transverse  and 
circumferential  velocity  components.  We  define  e  —  l/v/2We.  Within  the  present  analy¬ 
sis,  V  =  0  and  W  =  2eR  (R2  —  0.25)-0'5  «  2e;  the  latter  represents  the  nondimensional 
swirl- velocity  which  stabilized  the  undisturbed  annular  film  with  nondimensional  radius  R. 

The  linearized  system  of  equations  becomes 


dh  dh  du'  _  1,1  f  dh  t  dw‘) 

dF  +  d?  +  d?  ~  ~~RV  ~R\wde^~dfj 


(14) 
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du'  ^  du*  t  W  du* 
dt *  +  dz*  +  ~R~dO 
dv *  dv *  Wdv *  2„r  ,  W 2 

+  +^r7 
<9«;’  dw*  W  dw* 

~dF  +  ~dF  +  ^'W 

dr]  dr)  W  drj 

dF  +  d?  +  ~Rd6 


where  terms  of  0(R~ 3)  have  been  neglected  and  where  t*  denotes  nondimensional  time 
and  z*  represents  the  nondimensional  axial  coordinate.  Note  that,  the  radiant  9  is  already 
nondimensional.  Eqns.  (14)  to  (18)  apply  in  the  pressure  range  where  gas  density  is  negligibly 
small  and  the  effects  of  gas  inertia  can  be  neglected.  Solutions  to  Eqs.  (14)  -  (18)  are  assume 
to  be  of  the  form 


h  = 


u  = 


TJ  = 


V  — 


w  = 


g  i(kz’  +nB- ut’ ) 
AeHkx‘+n0-ur) 


A^eHkzm+n6-LJtm) 


(19) 


Substitution  of  Eqns.  (19)  into  Eqns.  (14)  -  (18)  provides  the  dispersion  relation  governing 
the  propagation  ot  capillary  waves  on  the  modulated  him.  i.e. 


12 


(a2  -  bg )  [2  ce  —  {g  +  e{f  -e/ 4))  (2/ -  e)] 


+  cnk2  [c1'2  (a2  4-  4<?)  4-  n  a  (g  4-  e  (/  —  e/4))j  (20) 

-  ca  [2fe  +  (2f-e)2  +  4g]  =  0 

with 


a 

6 

c 


/ 

5 


-n-{u-  k) 

fl 

R2 

oj  —  k 


Note  that  fc  and  n  denote  the  wavenumbers  of  capillary  waves  propagating  in  the  axial 
and  circumferential  direction,  respectively,  u  represents  the  angular  frequency  at  which 
modulations  are  forced  onto  the  discharging  film  at  the  atomizer  exit. 

Analogous  to  previous  work  by  the  authors,  evaluation  of  the  group  velocity  for  wavenum¬ 
bers  ki,  i.e.  C{ki)  =  dcj/dRe[A:,]  can  now  be  used  to  determine  which  disturbances  generated 
by  the  nozzle  modulation  at  frequency  ui  will  be  found  downstream.  In  particular,  if  the  group 
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velocity  for  a  particular  axial  wavenumber  is  positive,  then  that  wavenumber  will  be  found 
further  downstream,  if  the  particular  wavenumber  has  negative  group  velocity,  then  the  cor¬ 
responding  disturbance  will  not  be  found  downstream  from  the  nozzle  since  no  modulation 
is  applied  downstream.  The  latter  is  effectively  a  representation  of  the  Sommerfeld  radiation 
condition  which  implies  that  no  energy  is  propagated  upstream  from  infinity  [7,  17].  Note 
that  after  film  rupture  the  discharging  film  will  have  finite  length  which  necessitates  a  recon¬ 
sideration  of  the  downstream  boundary  condition.  In  that  case,  the  rejection  of  wavenumber 
solutions  with  negative  group  velocity  based  on  the  application  of  the  Sommerfeld  radiation 
condition  is  no  longer  valid. 

Equation  (20)  is  the  dispersion  relation  for  linear  three-dimensional  capillary  waves  on 
thin  swirling  annular  films  discharging  into  a  void  and  in  a  reference  frame  fixed  to  the  nozzle 
exit. 

Assuming  exponential  solutions  for  h,  tj,  u*,  v*  and  w*  according  to  Eq.  (19),  the  general 
solution  to  the  boundary- value  problem  for  a  given  forcing  frequency  u>  and  circumferential 
wavenumber  n  is  given  by 


h(z\f ) 

V  (z*,0 

n  l :  .  M 

vm{z\n 


y  A^ei(kzm+n$-uitmwtm) 

i=  i 

5 

J=1 

5 
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=  E^-Je,(fcz‘+nfl““‘’ur) 

J=1 


where  the  index  j  =  1, . . .  5  refers  to  the  5  wavenumber  solutions  to  the  dispersion  relation  for 
.  a  given  forcing  frequency  u  and  circumferential  wavenumber  n  which  have  to  be  considered 
within  the  present  analysis  according  to  the  employed  parameter  range  and  the  discussion 
of  the  previous  paragraph.  These  wavenumber  solutions  are  &2)  £3  and  k5  —  k7,  respectively. 
The  unknown  coefficients  within  the  above  solutions  depend  on  each  other  according  to 
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with  Pj  =  [^n  —  (u;  —  kj)\  and  qj  =  [kjS  +  -^-n2]  for  j  —  1, . . . ,  5.  The  above  system  of 
equations  is  obtained  after  substituting  Eqs.  (21)  into  Eqs.  (14)-(18).  Eqs.  (23)  are  linearly 
dependent  since  det[JW]  =  0  and  can  be  reduced  to 


Ahj  _  fj  +  4 e3Aj/[flpj] 

AvJ  £Pj/2  ~  ^kj/Pj[Qj  ~  R~2} 


= _ L  £V  /2  1  - 
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(27) 


Aw-  j  _  n[g;-  —  R  2][s2fc2  +  en/R(u  -  kj)  -  (ui  -  kj)2/i\  -  2 e(u>  -  k2)jR 
A, jj  ~  en/Rfa  -  R~2]-pj 

with  pj  and  qj  as  defined  above  and  Sj  =  n£2(k2+2R~2)-(e/R) (u—kj)(l—n)—(n/4)(ui—kj)2. 
This  provides  four  independent  equations  for  the  five  unknown  coefficients  of  each  index 
j  —  1, . . . ,  5.  Altogether,  this  yields  twenty  independent  equations  for  twenty-five  unknown 
coefficients.  The  remaining  five  conditions  needed  in  order  to  determine  fully  the  unknown 
coefficients  are  given  by  the  boundary  conditions  specified  at  the  nozzle  exit,  i.e.  one  condition 
for  each  of  the  five  wavenumbers  k 2,3  and  &5_7  with  positive  group  velocities. 

Results  and  Discussion 

Linear  Analysis 

Relevant  Wavenumber  Solutions  /  Wave  Character 

Eq.  (20)  has  been  solved  for  wavenumbers  ki  at  given  values  of  Weber  number  We  or 
s  and  circumferential  wavenumber  n.  Solutions  to  the  dispersion  relation  were  obtained 
numerically  by  using  Mathematical.  Evaluation  of  the  corresponding  group  velocities 
C(ki)  =  dcL//dRe[A;j]  predicts  that  for  u  <  0.6  and  small  n- values  (<  8  for  We  =  5  —  1000 
and  beyond),  information  travels  in  the  same  direction  as  in  the  axisymmetric  case:  Five 
wavenumbers  correspond  to  waves  that  originate  at  the  nozzle  and  send  information  in  the 
downstream  direction  (later  denoted  by  /c2, 3  for  dilational  mode  waves  and  k$-~  for  sinuous 


mode  waves)  and  two  (dilational  mode)  wavenumbers  correspond  to  waves  which  originate  at 
downstream  infinity  that  send  information  upstream  (later  denoted  by  ki,k4),  according  to 
their  positive  and  negative  group  velocities,  respectively.  The  subscript  i  is  an  integer  index 
that  indicates  the  particular  wave  number.  Previous  analyses  of  non-swirling  axisymmetric 
.  discharging  annular  films  indicated  the  existence  of  only  four  wavenumbers  with  positive 
group  velocity  [17].  In  the  present  case,  modulation  of  the  swirl  velocity  allows  for  the 
existence  of  one  additional  wavenumber  solution  with  positive  group  velocity. 

Figs.  2,  3  and  4  illustrate  the  solutions  to  Eq.  (20)  obtained  by  Mathematical^  plotted 
over  n  and  u  for  e  =  1/ ■'/lO  and  R  =  10.  Each  of  the  various  computed  solutions  might 
represent  a  different  wavenumber  in  a  different  parameter  range  (u/,  n).  After  combination  of 
the  various  branches  from  the  different  solutions  in  order  to  represent  the  seven  wavenumber 
solutions  ki(i  =  1, . . . ,  7),  one  observes  that,  wavenumber  plots  k^ui)  for  non-axisymmetric 
cases  with  n  <  5  will  resemble  the  axisymmetric  result  (n  =  0)  qualitatively.  The  described 
collection  of  branches  into  the  different  wavenumber  solutions  was  guided  by  comparison  with 
the  wavenumber  solutions  previously  determined  for  non-swirling  annular  films  (see  Figs.  3 
and  4  of  Ref.  [17]).  Fig.  5  displays  wavenumbers  fct-  as  a  function  of  u  for  the  n  =  3  case. 
For  uj  <  0.7  dashed  and  solid  lines  illustrate  predominantly  dilational  mode  waves  (Jk1_4), 
whereas  dotted  and  dash-dotted  lines  denote  predominantly  sinuous  mode  waves  (Aj5i6  and 
kj).  For  larger  Weber  number  cases  (e.g.  We  =  1000)  and  n  <  8  the  functional  dependence 
of  wavenumbers  A;,  on  forcing  frequency  is  similar  to  the  We  =  5  case  illustrated  here,  even 
though  absolute  values  might  be  quite  different.  Note  that  pure  dilational  or  sinuous  waves 
only  appear  as  the  planar  film  limit  R  — *  co  is  approached.  Prediction  of  the  wave-character 
(i.e.,  dilational  or  sinuous)  is  made  via  evaluation  of  the  amplitude  ratio  h  ° /rft  which  is 
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independent  of  the  imposed  upstream  (or  nozzle)  boundary  conditions,  i.e. 


(gj  -  kj)2' 
4s 2 


(28) 


Evaluation  of  Im(/:j]  and  consideration  of  C(kj )  =  da;/dRe[A:i]  shows  that  for  the  case 
illustrated  in  Fig.  5  no  sinuous  mode  wavenumber  (dotted  and  dash-dotted  lines)  is  to  be 
rejected  (i.e.  C(k 5_7)  >  0)  and  no  unstable  “sinuous”  mode  wavenumber  exists.  The  same 
is  true  for  the  axisymmetric  mode  and  non-axisymmetric  modes  (with  n  <  8)  of  the  similar 
case  (i.e.  R  =  10,  s  =  1/VlO  or  We  =  5)  (see  Figs.  2-4).  Only  one  of  the  four  “dilational” 
mode  wavenumbers,  namely  £3,  produces  exponentially  growing  solutions  in  the  downstream 
direction.  In  Fig.  5,  this  wavenumber  case  is  indicated  by  the  solid  line. 

Fig.  6  illustrates  the  growth  rate  for  the  described  unstable  predominantly  dilational 
wave  ks  as  a  function  of  modulation  frequency  and  circumferential  wavenumber  n  for  the 
particular  case  of  We  =  1000  (or  e  =  l/v^OOO)  and  R  =  10.  From  Fig.  6  and  similar  results 
for  various  other  combinations  of  We  and  R,  we  observe  that  for  large  Weber  numbers,  the 
instability  range  of  the  unstable  dilational  mode  wave  decreases  with  increasing  n-values.  For 
unstable  behavior  and  large  u  values,  the  growth  rate  decreases  with  increasing  n  values. 
In  fact,  for  large  enough  n-values,  no  film  instability  might  exist.  However,  for  unstable 
behavior  at  small  u  values,  the  growth  rate  is  found  to  increase  with  increasing  values  of  n 
until  for  large  enough  n-values,  the  growth  rate  drops  down  to  zero. 

As  r hr  Weber  number  is  decreased  fat  fixed  /?-values'].  the  maximum  growth  rare  of  the 
unstable  dilational-mode  wave  is  increased.  Compare  Fig.  6.  Fig.  7  and  Fig.  8  for  We  = 


1000.50  or  5  and  R  =  10.  respectively.  Fig.  8  also  illustrates  that  for  small  Weber  numbers. 
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and  small  enough  n- values,  the  instability  range  of  non-axisymmetric  modes  (n  ^  0)  is  larger 
than  for  the  n  =  0  case.  However,  as  in  Figs.  6  and  7,  for  a  given  parameter  set  R  and  We, 
the  maximum  growth  rate  is  still  observed  for  the  axisymmetric  (n  =  0)  case. 

For  the  same  case  as  depicted  in  Fig.  7,  Fig.  9  shows  the  dependence  of  the  wave  angle 
_d#,/dz  =  —kx/n  for  the  unstable  predominantly  dilational  wave  (i  —  3)  on  the  circumferential 
wavenumber  n.  The  magnitude  of  the  wave  angle  is  found  to  decrease  with  increasing  n- 
values  and  at  fixed  w- values,  except  for  small  values  of  ui.  Also,  |d0/dz|  increases  essentially 
linearly-  with  u  at  fixed  n,  except  for  small  u  values.  In  the  latter  case,  the  considered 
unstable  dilational-mode  wave  has  negative  wave  velocity  but  positive  group  velocity,  so 
that  the  winding  of  the  helical  structure  generated  by  the  wave  changes  from  clockwise  to 
counter-clockwise.  This  behavior  has  also  been  observed  within  the  numerical  analysis  of 
the  corresponding  initial-and-boundary- value  problem. 

Linear  Field  Solutions 

Solutions  to  the  linear  boundary-value  problem  discussed  earlier  were  obtained  for  the 
case  where  the  five  boundary  conditions  at  z*  =  0  are  given  by  harmonic  variations  imposed 
onto  variables  h,  tj,  u  ,  v*  and  w*  according  to  ,JJt  )  with  complex  forcing  amplitudes 
denoted  by  h®,  and  w® .  For  the  purpose  of  validating  the  numerical  analysis  pre¬ 

sented  below  and  for  the  purpose  of  identifying  the  influence  of  nonlinear  effects  on  the 
growth  of  unstable  waves  (via  comparison  between  linear  analytical  and  nonlinear  numer¬ 
ical  results’),  the  linear  problem  was  solved  for  the  case  where  onlv  the  unstable  dilational 
mode  wavenumber  k^  is  being  generated  at  the  nozzle  through  modulation  of  the  film  thick¬ 
ness.  i.e.  h,(:  =  0.  t  )  =  h°  cos  (nQ  —ut*).  Accordingly,  the  above  equations  are  solved 
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for  the  complex  forcing  amplitudes  t]°,u°,v°  and  w°.  In  particular,  coefficients  pertaining 
to  wavenumbers  other  than  k$,  i.e.  A^^,  A^j^,  Av-^2  and  Av-j^z ,  are  zero  and 

Ah,i  =  h°.  Coefficients  AVt2,  Au- (2,  Av- >2  and  Aw*t2  are  readily  determined  from  Eqns.  (24) 
through  (27).  Table  1  lists  the  various  nonzero  coefficients  (normalized  by  h°)  for  the  cases 
which  have  been  used  to  compare  linear  analytical  and  nonlinear  numerical  solutions. 

Nonlinear  Analysis 

The  governing  equations  for  the  nonlinear  three-dimensional  distortion  of  swirling  annular  or 
conical  liquid  sheets  have  been  solved  numerically  by  using  the  Law-Wendroff  Method  with 
Richtmyer  splitting  [22].  Initial  conditions  for  the  transient  simulations  of  swirling  and  non- 
swirling  annular  films  are  trivial  (A r(z,  t  =  0)  =  Ar0,  f(z,  t  =  0)  =  r0).  The  initial  film  shape 
for  swirling  conical  films  (discharging  from  the  nozzle  with  more  swirl  than  needed  to  stabilize 
the  film  in  its  annular  positions)  was  obtained  by  integrating  the  nonlinear  axisymmetric 
film  equations  in  the  downstream  direction  starting  at  the  nozzle  exit.  See  Ref.  [3]  for  more 
details.  Modulations  of  the  dependent  variables  Ar*,r*r,  Uq,Vq  and  Wq  were  imposed  at  the 
nozzle  exit  (z*  =  0)  in  order  to  generate:  1)  The  unstable  “dilational”  mode  wave  predicted 
by  linear  theory  (annular  case  only)  and  2)  a  system  of  standing  and/or  travelling  waves 
(annular  and  conical  case).  The  former  provides  a  dynamically  simple  system,  which  eases 
the  comparison  between  linear  and  nonlinear  theories.  The  latter  case  has  been  considered  in 
context  with  the  active  control  of  the  film  disintegration  process;  it  represents  a  dynamically 
complex  system. 

For  both  swirl-stabilized  annular  and  ‘conical’  sheets,  the  number  of  boundary  conditions 
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at  the  nozzle  exit  was  chosen  according  to  the  number  of  wavenumbers  hi  with  associated 
positive  group  velocities,  relevant  to  the  linear  boundary-value-problem  analysis  presented 
earlier  for  swirl-stabilized  annular  films,  i.e.  five  boundary  conditions  were  imposed  at  z*  =  0. 
Analogous  to  the  analysis  of  swirling  axisymmetric  films  [3],  boundary  conditions  at  the  noz¬ 
zle  exit  for  nondimensional  sheet  thickness  A r*,  radial  sheet-centerline  location  r*,  and  the 
velocity  components  in  the  axial,  radial  and  circumferential  direction  Uq,Vq  and  Wq  were 
chosen  according  to  the  described  steady-state  solutions,  but  with  additional  harmonic  vari¬ 
ations  of  the  film  thickness  or  the  axial,  circumferential  and  transverse  velocity  components 
generating  axisymmetric  or  three-dimensional  standing  and/or  travelling  sinuous  or  dila¬ 
tions!  mode  waves. 

As  in  Ref.  [3],  additional  numerical  boundary  conditions  (required  to  solve  the  unsteady 
problem)  were  specified  for  d2A r* /dz* 2  and  d2r*/dz*2  with  values  corresponding  to  the 
imposed  steady-state  initial  conditions. 

In  summary,  for  the  analysis  of  the  considered  unsteady  problem,  five  boundary  condi¬ 
tions  and  two  numerical  conditions  were  specified  at  the  nozzle  exit. 

Parameter  Range  and  Forcing  Conditions 

Analogous  to  previous  work  by  the  authors,  and  in  order  to  omit  the  dominance  of  tran¬ 
sient  effects  in  the  majority  of  the  nonlinear  numerical  simulations,  the  present  numerical 
study  was  limited  to  large  Weber  number  flows,  i.e.  We  =  50  or  We  =  1000.  Values  for  the 
radius-to-thickness  ratio  of  the  discharging  films,  as  well  as  values  for  the  imposed  circum¬ 
ferential  wavenumber  and  forcing  frequencies  (determining  downstream  and  circumferential 
disturbances)  were  chosen  small  enough  to  guarantee  film  rupture  and/or  the  appearance  of 
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significant  nonlinear  three-dimensional  effects  within  a  computationally  feasible  time  frame. 
In  addition  the  same  parameters  were  chosen  large  enough  so  that  the  assumptions  of  the 
employed  lubrication  model  are  not  violated.  For  all  the  results  presented  here,  the  nondi- 
mensional  annular  radius  has  been  fixed  at  R  =  10.  Results  computed  for  cases  with  n  =  1 
have  not  been  included  due  to  the  delayed  appearance  of  three-dimensional  effects  even  at 
circumferential  wave  amplitudes  comparable  to  those  of  the  downstream-propagating  waves. 
The  analysis  of  conical  films  was  limited  to  cases  with  swirl  number  k  =  10,  indicating  that 
the  liquid  exits  the  annular  slit  nozzle  with  ten  times  the  amount  needed  to  stabilize  it  in  the 
annular  configuration.  The  steady-state  cone  angle  in  this  case  is  approximately  44  degrees. 

Disturbance  amplitudes  in  the  various  flow  variables  were  chosen  over  a  wide  range.  Small 
amplitude  values  were  used  in  order  to  benchmark  the  numerical  model  by  allowing  direct 
comparison  between  linear  analytical  and  numerical  results  for  a  given  set  of  boundary  or 
forcing  conditions.  Larger  amplitude  values  were  chosen  in  order  to  predict  the  appearance 
of  nonlinear  effects  and  possibly  film  rupture  near  the  nozzle  exit.  As  noted  before,  early 
appearance  of  nonlinear  effects  and  film  rupture  were  essential  in  order  to  perform  the  nu¬ 
merous  nonlinear  quasi-three-dimensional  numerical  simulations  in  a  reasonable  time  frame 
and  at  reasonable  computational  costs.  Note  that,  in  general,  practical  liquid  film  atomizers 
are  designed  (among  other  things)  to  achieve  film  rupture,  disintegration  and  atomization 
at  short  distances  from  the  atomizer  nozzle.  The  latter  being  an  important  requirement  in 
the  design  of  more  compact,  i.e.  shorter  and  lighter,  combustion  systems. 

Here,  we  focus  on  1)  a  comparison  of  linear  and  nonlinear  unstable  dilational  mode  waves 
on  swirling  annular  films  and  2)  a  description  of  film  topology  and  film  rupture  under  a  variety 
of  forcing  conditions  resulting  in  standing  and/or  travelling  sinuous  and/or  dilational  waves 
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via  modulation  of  various  different  flow  variables,  i.e.  (f>  =  77,  ujj,  ujj,  Wg.  The  general  forcing 
condition  for  one  of  the  flow  variables  <f>  employed  here  is  given  by 

*  =  *  +  *'  [l  -  e|-‘‘/T->]  f,(~  +  fMnS  -  j  ^t)  /3(^)  (29) 

Here,  $  represents  the  steady-state  value  of  the  particular  flow  variable,  fi,  /2  and  /3  denote 
either  sine  or  cosine  functions.  The  parameter  j  takes  the  value  0  or  1  in  order  to  model 
a  standing  or  a  travelling  wave,  respectively.  The  exponential  term  [1  —  assures 

a  smooth  start-up  of  the  transient  simulation  preventing  the  film  distortion  process  from 
being  dominated  by  the  dynamics  of  the  initial  wave-front  propagating  into  the  undisturbed 
semi-infinite  film  [23].  Te  denotes  the  characteristic  time  for  the  described  start-up  process. 
The  various  parameter  combinations  and  forcing  conditions  which  have  been  analyzed  are 
summarized  in  Tables  2  through  8,  together  with  data  on  film  break-up  time  and  length  for 
the  various  cases.  If  film  rupture  did  not  occur  within  the  computational  time  frame,  min¬ 
imum  film  thickness  and  corresponding  downstream  location  (at  the  end  of  the  simulation) 
are  tabulated  instead. 

Single  “Dilational”  Wave  Modulation 

Figs.  10  and  11  show  the  instantaneous  film  thickness  as  a  function  of  downstream  dis¬ 
tance  for  a  swirling  annular  film  with  We  =  1000,  R  —  10  and  n  =  0  or  n  =  3,  respectively. 
The  forcing  frequency  in  both  cases  was  u  =  0.26  with  h0  =  0.35.  Boundary  conditions  have 
been  specified  in  order  to  generate  the  previously  discussed  unstable  dilational  wave  only. 
For  the  axisymmetric  case  illustrated  in  Fig.  10  the  nondimensional  film  breakup  length  (i.e. 
length  measured  from  the  nozzle  exit  to  the  first  point  of  film  rupture)  is  rb  =  390.  The 
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rupture  occurs  at  tb  =  422.85  measured  in  nondimensional  terms.  Comparison  with  linear 
analysis  shows  that  the  nonlinear  growth  rate  (at  z*  =  300)  is  about  56  %  larger  than  pre¬ 
dicted  by  linear  theory.  For  the  non-axisymmetric  mode  shown  in  Fig.  11,  breakup  length 
and  time  were  predicted  a s  l*b  =  245  and  t'b  =  275.55,  respectively.  Here,  the  increase  of  the 
film  thickness  A rmai  is  approximately  the  same  as  for  the  n  =  0  case.  On  the  other  hand, 
the  decrease  of  Armtn  in  the  downstream  direction  is  larger  for  the  n  =  3  case  resulting  in 
shorter  film  break-up  time  and  length.  Wave  angle  measurements  obtained  from  the  the 
nonlinear  numerical  simulation  shown  in  Fig.  11  agree  very  well  with  the  wave  angle  pre¬ 
diction  d9/dz  =  —k/n  obtained  from  linear  theory.  In  other  words,  constant  film-thickness 
values  are  found  along  z  =  —{n/k)6  lines. 

Figs.  12  and  13  demonstrate  film  distortion  and  rupture  for  a  a  lower  Weber  number  case, 
i.e.  We  =  50,  at  the  same  annular  radius,  forcing  frequency  and  circumferential  wavenumber 
as  in  Fig.  11  (R  =  10,  u  =  0.26,  n  =  3).  The  disturbance  amplitude  in  the  film  thickness 
for  this  case  was  h0  =  0.15.  The  lower  Weber  number  results  in  an  increased  growth  rate 
of  the  modulated  unstable  “dilational”  wave.  Consequently,  breakup  length  and  time  are 
significantly  smaller  than  for  the  similar  larger  Weber  number  case.  In  particular,  ZJ  =  95 
and  tl  =  115.2  for  the  case  shown  in  Figs.  12  and  13.  Observed  growth  of  A rmax  and  decrease 
of  A rmin  in  the  downstream  direction  correspond  to  the  observation  previously  described  for 
the  larger  Weber  number  case  (Fig.  11).  Also,  distortion  characteristics  are  not  different 
from  the  larger  Weber  number  case.  It  is  noted  here  that,  for  the  similar  axisymmetric  case 
(n  =  0),  film  rupture  does  not  occur  within  the  simulated  time  frame  (0  <  t"  <  150). 

The  present  nonlinear  analysis  considers  an  initial-  and  boundary-value  problem.  Film 
modulation  is  started  at  t’  =  0  after  which  the  front  of  the  film  disturbance  is  propagated 


24 


downstream  into  the  undisturbed  flowing  film.  As  discussed  earlier,  the  temporal  growth  of 
the  propagating  initial  wave-front  might  dominate  over  the  spatial  wave  growth;  the  latter 
resulting  from  continuous  modulation  at  a  prescribed  forcing  frequency.  Comparisons  of  the 
corresponding  film  rupture  lengths  and  times  with  predictions  from  linear  spatial  theory 
are  not  appropriate.  The  dominant  behavior  of  the  described  transient  effect  is  observed  in 
Figs.  14,  15  and  16.  for  n  =  0  or  n  =  3  with  R  =  10,  u  =  0.06,  We  =  50  and  A r*°  =  0.15. 
Even  as  transient  effects  dominate  the  film  distortion  process,  non-axisymmetric  disturbances 
provide  shorter  break-up  length  and  break-up  time.  For  the  axisymmetric  case  in  Fig:  14, 
film  rupture  does  not  occur  within  the  simulated  time  frame,  i.e.  0  <  t*  <  250.  For  the 
non-axisymmetric  n  =  3  case  illustrated  in  Figs.  15  and  16,  film  breakup  (dominated  by 
“start-up”  effects)  is  observed  at  l*b  =  117.5  adn  t*b  =  127.2.  At  the  time  of  rupture,  film 
distortion  is  highly  nonlinear,  particularly  close  to  the  propagating  wave  front. 

Multiple  Dilational  Wave  Modulation  /  Wave  System 

In  order  to  investigate  the  feasibility  of  actively  controlling  the  film  disintegration  process, 
film  distortion  and  film  rupture  were  investigated  for  the  case  where  multiple  superimposed 
waves  are  generated  at  the  nozzle  exit. 

Non-swirling  Annular  Films 

Fig.  17  illustrates  the  case  of  a  dilationally  modulated  non-swirling  (pressure-stabilized) 
annular  film  with  a  harmonic  forcing  imposed  at  z*  =  0  according  to 


uQ{f)  =  l+A 


-r 
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cos  (no)  sin  (— — ) 
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where  t*  denotes  the  nondimensional  time  variable,  and  9  is  measured  in  radians.  For  non¬ 
swirling  annular  sheets,  stabilization  of  the  undisturbed  film  is  achieved  by  a  nonzero  constant 
pressure  difference  across  the  film,  i.e.  pg-  -pg,+  =  cr(l/r0,+  +  l/r0,_)  =  2ar0/(rl  -  Arjj/4), 
where  the  subscript  ‘O’  denotes  undisturbed  values  at  the  nozzle  exit.  The  Weber  number 
for  this  case  is  We  =  1000,  the  undisturbed  annular  film  radius  Fj  =  10,  Tp  =  Te  =  10  and 
A  =  B  =  0.04.  Five  (standing)  waves  have  been  imposed  in  the  circumferential  direction,  i.e. 
n  =  5.  As  observed  from  Fig.  17,  the  imposed  film  modulation  results  in  the  development 
of  a  regular  cellular  structure  on  the  film  in  the  downstream  direction  from  the  nozzle  exit. 
As  the  amplitude  of  the  film  distortion  increases,  a  system  of  fluid  ligaments  forms.  The 
ligaments  are  elongated  in  the  circumferential  direction  and  are  connected  by  a  thinner  film 
of  liquid.  Local  maxima  and  minima  in  the  film  thickness  alternate  in  the  circumferential 
direction  and  also  in  the  downstream  direction  at  constant  0-value.  Points  where  the  film 
breaks  first  (i.e.  where  it  reaches  zero  thickness)  are  found  at  the  same  downstream  location 
z*  and  located  near  the  transitions  between  the  thin  film  regions  and  the  thicker  ligaments 
resulting  in  the  simultaneous  formation  of  larger  stretched  free  liquid  ligaments  inter-spaced 
by  thinner  fluid  films.  Both  larger  ligaments  and  thinner  fluid  fluid  films  are  expected 
to  contract  towards  a  spherical  shape  due  to  surface  tension  (after  being  generated  from 
the  continuous  film)  and  oscillate,  due  to  the  absence  of  viscous  damping,  around  their 
spherical  equilibrium  shapes.  A  simple  theoretical  analysis  considering  the  aspect  of  viscous 
dissipation  in  liquid  atomization  systems  can  be  found  in  Ref.  [24 j.  Note  that  the  current 
implementation  does  not  allow  for  continuation  of  the  simulation  beyond  the  point  when  the 


film  breaks  at  some  location(s). 

Fig.  18  illustrates  the  similar  (pressure-stabilized)  case  as  shown  in  Fig.  17  but  now  with 
a  three-dimensional  sinuous  modulation  enforced  at  the  nozzle.  The  behavior  is  generated 
by  harmonic  variations  of  the  transverse  film  velocity  according  to 

<>o (O  =  C  [l  -  exp  cos  (nd)  sin  (^-)  (32) 

with  C  =  0.1,  Tp  =  Te  =  10  and  n  =  5.  In  this  case,  the  initially  regular  three-dimensional 
sinuous  disturbance  flattens  in  the  downstream  direction,  generating  thicker  fluid  rings  (with 
circumferential  thickness  fluctuations)  which  are  connected  by  thin  films  of  liquid.  The  ob¬ 
served  change  from  a  three-dimensional  sinuous  to  a  more  or  less  axisymmetric  dilational 
structure  results  from  the  development  of  nonlinear  sinuous-dilational  mode  coupling  down¬ 
stream,  but  is  also  related  to  the  presence  of  multiple  downstream  propagating  waves  result¬ 
ing  in  an  envelope  behavior  for  film  thickness  and  film  centerline  location  in  the  downstream 
direction.  The  latter  has  already  been  described  for  non-swirling  axisymmetric  annular  films 
[17]  and  was  observed,  for  example,  also  in  cases  2  and  3  of  Table  6.  As  in  the  dilationally 
modulated  case,  film  rupture  is  observed  first  between  the  larger  fluid  rings.  The  fluid  rings 
themselves  can  be  expected  to  break  due  to  the  dominance  of  the  most  unstable  wavelength 
analogously  to  the  Rayleigh  instability  mechanism  for  straight  liquid  jets  [25,  26]. 

Swirling  Annular  Films 

Fig.  19  illustrates  the  film  thickness  distribution  for  a  non-swirling  pressure-stabilized 
annular  liquid  film,  when  the  front  of  the  disturbance  generated  by  the  modulation  has 
propagated  to  about  z'  —  35  in  the  downstream  direction.  Problem  parameters  for  this  case 
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axe  the  same  as  for  the  case  illustrated  in  Fig.  18.  The  film  was.  modulated  according  to  the 
forcing  conditions  prescribed  in  Eq.  (32).  Fig.  20  shows  the  instantaneous  thickness  distri¬ 
bution  at  the  same  time  for  the  corresponding  swirl-stabilized  case.  Comparison  of  the  film 
distortion  for  pressure-stabilized  and  swirl-stabilized  films  shows  no  significant  differences  in 
the  particular  film  thickness  distribution.  However,  as  already  observed  for  axisymmetric 
films  [3],  disturbance  amplitudes  are  significantly  larger  in  the  swirl-stabilized  case,  reducing 
the  length  between  nozzle  exit  and  downstream  position  where  film  rupture  first  occurs. 

Swirling  Conical  Films 

Fig.  21  shows  the  outer  film  surface  location  for  the  case  of  a  swirling  conical  sheet 
generated  by  liquid  which  exits  the  annular  slit  nozzle  with  ten  times  the  amount  of  swirl 
needed  to  stabilize  the  film  in  its  annular  configuration  when  exiting  the  nozzle.  In  other 
words,  the  swirl  velocity  for  this  case  is  ten  times  larger  than  for  the  case  illustrated  in 
Fig.  20. 

Initial  conditions  for  this  case  were  obtained  by  solving  the  nonlinear  steady-state  ax¬ 
isymmetric  equations  as  described  in  Ref.  [3].  The  liquid  film  is  modulated  sinusoidally  at 
the  nozzle  (z*  =  0)  according  to  Eq.  (32).  The  Weber  number  for  this  case  is  We  =  1000, 
R  =  10  and  Tv  =  Te  =  10.  As  in  the  previous  case,  5  standing  waves  were  imposed  at  the 
nozzle  exit,  i.e.  n  =  5,  with  C  =  0.1. 

For  the  same  case  and  at  the  same  time,  Fig.  22  shows  the  instantaneous  film  thickness 
distribution.  Tnterestinglv.  areas  of  the  film,  with  a  local  maximum  of  the  film  thickness  at  the 
nozzle  exit  develop  into  areas  with  two  local  maxima  in  the  film  thickness.  The  latter  is  not 
found  for  the  similar  swirling  annular  sheet  of  Fig.  20  and  results  from  the  film  divergence  in 
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the  conical  case.  Film  divergence  due  to  excess  swirl  also  causes  the  regular  cellular  structure 
observed  at  the  nozzle  to  stretch  forming  a  net-like  or  web-like  configuration  of  thicker  fluid 
ligaments  imposed  onto  a  thinner  layer  of  liquid  (see  Fig.  23).  As  previously  described  for 
non-swirling  annular  films,  film  rupture  first  occurs  in  the  transition  regions  between  the 
thinner  liquid  layers  and  the  thicker  ligaments.  For  both  non-swirling  and  swirling  annular 
films  described  earlier,  initial  film  break-up  points  are  located  at  the  same  downstream 
position  and  more  or  less  symmetrically  on  both  sides  of  the  thin  films  or  thicker  ligaments 
In  the  swirling  conical  case,  this  symmetry  is  lost  and  initial  rupture  points  are  found  only 
on  one  side  of  the  thin  film  regions,  whereas  the  other  side  still  remains  connected  to  the 
neighboring  larger  liquid  mass.  Depending  on  the  flow  conditions,  the  fluid  within  the  thin 
films  might  be  subsequently  re-absorbed  into  the  larger  ligaments.  Note  that,  for  the  annular 
films  of  Figs.  19  and  20,  the  continuous  liquid  film  is  expected  to  disintegrate  initially  into 
thinner  liquid  layers  located  between  larger  fluid  ligaments.  However,  recombination  of  the 
larger  fluid  masses  with  the  smaller  ones  is  still  possible  due  to  different  relative  velocities 
of  their  center  of  gravity.  The  latter  is  analogous  to  the  recombination  of  main  and  satellite 
droplets  observed  in  liquid  jet  atomization  processes  such  as  ink-jet  printing  [27]. 

For  larger  disturbance  amplitudes,  film  breakup  occurs  close  to  the  nozzle  exit  (see 
Fig.  24)  with  film  distortion  being  highly  nonlinear  prior  to  breakup.  For  the  case  shown 
in  Fig.24  fluid  accumulates  into  spike-like  structures  with  the  possibility  of  drop  generation 
before  film  breakup  into  larger  ligaments  occurs.  It  is  conceivable  that  droplets  ejected 
into  the  inner  region  of  the  conical  film  will  impact  onto  the  film  causing  film  rupture  and 
atomization  at  the  downstream  location  where  they  have  impacted.  The  relevance  of  the 
impingement  of  droplets  onto  discharging  liquid  films  near  the  nozzle  exit  has  already  been 
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considered,  by  other  authors  [28]. 

Clearly,  in  order  to  investigate  liquid-film  dynamics  after  initial  break-up  of  the  continu¬ 
ous  film  is  predicted,  the  current  model  has  to  be  extended  to  describe  the  dynamics  of  the 
newly  generated  free  film  edges.  At  the  film  edges,  lubrication  equations  cannot  be  employed. 
Rather,  the  authors  envision  a  combination  of  a  dynamic-film-edge  model  in  combination 
with  the  current  lubrication  model  away  from  the  edges.  Limitations  of  this  combined  model 
are  reached  when  the  thickness-to-length  ratio  of  the  retracting  free  liquid  ligaments  is  of 
0(1).  Also,  as  noted  earlier,  rejection  of  boundary  conditions  at  the  nozzle  exit  on  the  basis 
of  the  Sommerfeld  radiation  condition  can  no  longer  be  employed  after  rupture  of  the  semi¬ 
infinite  film  generates  free  ligaments  and  a  finite-length  continuous  film  is  discharging  from 
the  nozzle  exit. 


Summary 

Linear  and  nonlinear  analyses  of  modulated  three-dimensionally  distorting  thin  inviscid  free 
liquid  films  discharging  into  a  gas  of  negligible  density  are  presented.  The  nonlinear  numer¬ 
ical  analysis  uses  a  lubrication  model  reducing  the  three-dimensional  problem  to  a  system 
of  two-dimensional  unsteady  equations.  Linear  theory  for  swirling  annular  films  predicts 
that  for  We  >  2  only  one  unstable  wave  is  generated  on  the  film  due  to  its  modulation  at 
the  nozzle  exit.  Depending  on  Weber  number  and  film  radius,  linear  growth  rates  for  this 
unstable  dilacional  mode  wave  might  lie  larger  for  iioii-axisymmetne  modes  than  for  the 
corresponding  axisymmetnc  case.  However,  for  the  considered  parameter  range,  maximum 
growth  rates  for  a  given  Weber  number  and  annular  radius  are  still  observed  for  the  ax- 
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isymmetric  case.  Modulation  of  multiple  dilational  mode  waves  at  the  nozzle  exit  allows  for 
uniform  breakdown  of  swirling  and  non-swirling  annular  films  indicated  by  the  generation 
of  uniformly  sized  liquid  volumes  which  are  uniformly  spaced  and  connected  by  thinner  liq¬ 
uid  layers.  Film  topology  and  break-up  of  three-dimensionally  modulated  conical  films  are 
discussed.  Comparison  with  the  corresponding  swirling  annular  film  shows  that  for  sinuous 
film  modulation,  film  divergence  causes  the  separation  of  initially  formed  larger  fluid  blobs 
into  a  pair  of  smaller  fluid  volumes  connected  by  an  even  thinner  liquid  layer. 

Within  the  considered  large  Weber  number  and  low  forcing  frequency  range,  wave  an¬ 
gle  measurements  from  the  nonlinear  numerical  simulations  agree  very  well  with  predictions 
&Q/&Z  =  —k/n  obtained  from  linear  theory  for  both  sinuous  and  dilational  waves  on  swirling 
annular  films.  On  clockwise  swirling  conical  films,  spiraling  dilational  and  sinuous  waves 
moving  in  the  clock-wise  direction  increase  in  slope  as  the  film  thins  out  in  the  downstream 
direction.  The  slope  of  counter-clockwise  propagating  waves  decreases  with  downstream  dis¬ 
tance.  The  described  changes  in  slope  can  be  attributed  to  a  decrease  in  swirl  velocity  with 
increase  in  annular  film  radius  downstream  that  results  from  the  conservation-of-angular- 
momentum  principle.  Based  on  the  location  of  the  initial  film  rupture  points  and  the  thick¬ 
ness  distribution  at  the  time  of  film  rupture,  Figs.  25  through  29  illustrate  projected  break-up 
patterns  characteristic  for  the  cases  with  pure  standing  or  travelling  dilational  or  sinuous 
mode  waves  modulated  at  the  nozzle.  For  mixed  standing/travelling  and  sinuous/dilational 
wave  modulation,  the  resulting  break-up  pattern  are  modified  or  mixed  versions  of  the  pat¬ 
terns  illustrated  in  Figs.  25  -  29.  Here,  filled  circles  indicate  locations  of  local  maxima 
in  the  film  thickness,  solid  lines  are  perpendicular  to  the  directions  along  which  the  various 
modulated  waves  (two  in  the  standing-wave  cases,  one  in  the  travelling- wave  cases)  are  prop- 
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agating.  Dashed  lines  represent  predictions  for  the  free  film  edges  generated  after  the  initial 
rupture  of  the  film.  The  illustrated  break-up  pattern  (dashed  lines)  are  derived  from  the 
assumption  that  after  the  initial  local  rupture  of  the  film  ,  the  free  liquid  edge(s)  propagates 
mainly  along  the  direction  where  the  film  thickness  is  the  smallest.  This  assumption  is  based 
on  the  observation  that  the  velocity  of  a  free  liquid-film  edge  is  given  by  ue  =  \j2a/(pt), 
where  a,  p  and  t  denote  surface  tension  coefficient,  density  and  thickness  of  the  uniform  film 
[29].  For  the  investigated  cases  with  travelling  dilational  or  sinuous  mode  circumferential 
waves  (wavenumber  n),  superposition  with  the  imposed  axial  wave  (wavenumber  k)  results 
in  an  oblique  wave  spiraling  clock-wise  downstream  on  the  annular  or  conical  film.  Here,  ini¬ 
tial  film  breakup  occurs  simultaneously  at  various  z  =  constant  points.  The  perforations  are 
expected  to  expand  along  the  line  of  minimum  film  thickness  (see  dashed  lines  in  Fig.  25  for 
the  annular  case)  which  results  in  the  formation  of  liquid  filaments  spiraling  downstream  in 
the  counter-clockwise  direction  (“filament  beak-up”).  On  swirling  annular  films  with  super¬ 
imposed  dilational  standing-wave  modulations  in  the  circumferential  direction,  film  rupture 
first  occurs  after  the  formation  of  larger  fluid  blobs  and  upstream  just  behind  these  larger 
fluid  masses.  Again,  initial  film  rupture  will  take  place  simultaneously  at  various  z  ^constant 
positions.  Subsequent  expansion  of  the  film  perforations  can  be  expected  to  generate  a  more 
or  less  circular  pattern  within  each  cell  formed  by  neighboring  fluid  blobs.  This  type  of  film 
rupture  can  be  characterized  as  “cellular  break-up”  and  is  illustrated  in  Fig.  26.  The  same 
break-up  pattern  can  also  be  observed  for  non-swirling  (pressure-stabilized)  annular  films 
and  for  swirling  conical  films  with  similar  forcing  conditions  (see  Fig.  27).  However,  due 
to  film  divergence  the  “cells”  will  stretch  in  the  conical  case  and  film  thickness  within  the 
cells  will  be  reduced.  Therefore,  expansion  of  the  film  perforation  is  expected  to  occur  sig- 
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nificantly  faster  than  in  the  corresponding  annular  case.  Note,  however,  that  initial  rupture 
of  the  film  in  the  conical  case  is  delayed.  Fig.  28  illustrates  the  break-up  pattern  expected 
for  the  considered  swirling  annular  cases  with  a  standing  sinuous  circumferential  wave  im¬ 
posed  at  the  nozzle  exit.  As  for  the  dilational  case  discussed  earlier,  initial  film  perforation 
_  takes  place  just  behind  the  larger  fluid  blobs.  However,  in  this  case  the  perforations  are  ex¬ 
pected  to  propagate  (more  or  less)  along  z  —  constant  lines  in  the  circumferential  direction 
(see  Fig.  28),  ultimately  resulting  in  the  detachment  of  fluid  rings  with  pronounced  thick¬ 
ness  fluctuations  in  the  circumferential  direction  (“ring  break-up”).  The  similar  standing 
wave  sinuous  mode  modulation  imposed  onto  the  conical  film  geometry  yields  the  “filament 
break-up”  pattern  shown  in  Fig.  29.  For  the  investigated  cases,  the  filaments  are  spiraling 
clock-wise  in  the  downstream  direction.  In  contrast  to  the  single  wave  modulation  of  Fig.  25, 
breakup  of  the  filaments  in  Fig.  29  will  be  greatly  influenced  by  the  existing  non-uniform 
mass  distribution  along  the  filaments. 
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Table  1:  Nonzero  normalized  coefficients  within  Eqs.  21  if  only  the  unstable  dilational  mode 
wave  with  wavenumber  kz  is  modulated  at  the  nozzle  exit. 


case 

We  = 1000 

We  =  50,  u 

=  0.26 

n 

=  0 

n 

=  3 

n  =  0 

n  —  3 

uj  =  0.06 

u  =  0.26 

uj  =  0.06 

w  =  0.26 

Re[&3] 

0.259998 

0.247686 

0.259957 

0.205644 

Im[Ar3] 

-0.00336706 

-0.00241779 

-0.0150983 

-0.0122502 

R e[AVt2/h°} 

4.07431 

1.07193 

-2.35187 

1.31047 

1.06706 

1.85022 

In^A^/h0] 

0.158797 

0.0202575 

1.91089 

0.102595 

0.0903446 

0.819221 

Re[At-,2/h°] 

-0.00383658 

-0.000203372 

-0.0102061 

0.00577779 

-0.00406128 

-0.000287701 

0.151768 

0.0143359 

0.00764552 

0.0169967 

0.0640431 

0.107195 

Re[-Av*^/h°] 

-0.0814861 

-0.00479382 

-0.117893 

-0.00604028 

-0.0213412 

-0.0679461 

Im^-,2  /h°] 

-0.003167593 

-0.0000905945 

-0.13906 

0.0121581 

-0.00180689 

0.0833378 

RefA^/h0] 

0.0262353 

0.0036093 

-0.0459605 

0.00305534 

0.0161147 

0.0180419 

Im  ^.,2/h0] 

-0.00143516 

0.0000660645 

0.00966504 

0.00169273 

0.00131817 

0.0204783 

Table  2:  Dilational  travelling  wave  modulation  on  swirling  annular  films  generated  via  har¬ 
monic  oscillations  of  axial  film  velocity  u*  [it!  =  10,  Tp  =  15,  We  =  1000,  fx=fz  =  sin, 
f2  =  cos,  j  =  1,  u*0  =  0.025,  tj  -  v*  =  u*  -  w*  =  0]. 
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3 
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2tt/3 
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4 

3 

20tt/3 

0.1 

1 

117.6 

97.5 

5 

3 

2tt 

0.1 

10 

114.6 

94.8 
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Table  3:  Dilational  travelling  or  standing  wave  modulation  on  swirling  annular  films  gen¬ 
erated  via  harmonic  oscillations  of  axial  film  velocity  u*  [J?  =  10,  Tp  =  15,  We  =  1000, 
/i  =  sin,  fz  =  1 ,  h  —  cos,  n  =  3,  Ttr  =  2i r,  Uq  =  0.025,  tj  =  v*  =  u*  =  w*  =  0]. 


case 

uh 

i 

k 

k 

0 

0.1 

l 

1 

118.1 

97.8 

1 

1.0 

l 

1 

105.7 

87.1 

2 

0.1 

0 

1 

119.4 

99.5 

3 

0.5 

0 

10 

116.9 

97.8 

Table  4:  Sinuous  travelling  wave  modulation  on  swirling  annular  films  generated  via  harmonic 
oscillations  of  transverse  film  velocity  v*  [i?  =  10,  Tp  =  25,  We  =  1000,  fi  =  fz  =  cos, 
h  =  COS,  Ttr  =  2tT,  Vq  =  0.1,  tj  =  u*  =  w*  =  0]. 


case 

n 

vh 

k 

£&,m 

0 

3 

0.1 

i 

142.6 

113.7 

1 

3 

0.5 

i 

131.7 

102.6 

2 

0 

0.1 

i 

145.1 

116.2 

3 

3 

0.5 

10 

181.7 

153 

Table  5:  Sinuous  standing  and  travelling  wave  modulation  on  swirling  annular  films  generated 
via  harmonic  oscillations  of  transverse  film  velocity  v*  [f?  =  10,  Tp  =  25,  We  =  1000,  f\  =  cos, 
fz  —  —  cos,  n  =  3,  Ttr  =  2-tt,  uj  =  0.1,  tj  =  u*  =  w*  =  0]. 


case 

u3  d 

3 

k 

^6,771 

^6,m 

0 

0.5 

1 

1 

139  2 

85.2 

1 

0.5 

0 

i 

107.1 

78.9 

2 

0.5 

0 

10 

145.6 

118.1 
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Table  6:  Sinuous  or  mixed  sinuous/dilational  standing  wave  modulation  on  non-swirling 
annular  films  generated  via  harmonic  oscillations  of  transverse  film  velocity  v*  or  both  film 
thickness  77  and  v\  [R  =  10,  Tp  =  10,  We  =  1000,  k  =  1,  /3  =  1,  j  =  0,  n  =  5,  u*  =  w*  =  0]. 


case 

h 

h 

v  0,  n 

0 

sin 

cos 

0.1” 

75 

39.8  (0.4) 

1 

sin 

1 

0.5” 

65 

50.3  (0.15) 

0 

cos 

0.2’ 

2 

sin 

sin 

0.5” 

100 

48.2  (0.23) 

sin 

cos 

0.5’ 

3 

sin 

sin 

0.4” 

75 

46.8  (0.05) 

sin 

cos 

0.4’ 

Table  7.  Sinuous  or  dilational  travelling  or  standing  wave  modulation  on  swirling  films 
generated  via  harmonic  oscillations  of  either  u*,  v'  or  w*.  [iZ  =  10,  Tv  —  Tx  =  10,  We  =  1000, 
/2  =  cos,  /1  =  0, 77  =  0]. 


case 

h 

3 

« (’),  0*  n,  u*  ( 

)  n 

k 

h,(m) 

0 

sin 

0 

0.25” 

5 

1 

37.5 

23.1  (0.14) 

1 

sin 

0 

0.25” 

5 

10 

36.4 

23.4 

2 

sin 

0 

0.15” 

5 

10 

56.3 

33.1 

3 

sin 

0 

0.1” 

5 

10 

78 

57.5  (0.1) 

4 

1 

1 

2 

2 

10 

70 

52.4  (0.015) 

5 

1 

1 

0.05’ 

2 

10 

77.5 

55  (0.019) 

6 

sin 

0 

0.05’ 

2 

10 

75 

52.5 

7 

sin 

0 

0.5 

4 

10 

27.6 

18.5 
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Table  8:  Modulation  of  single  unstable  dilational  mode  wave  onl  Dilational  travelling  wave 
modulation  on  swirling  annular  films  generated  via  harmonic  oscillations  of  axial  film  velocity 
u*  [i?  =  10,  Tp  =  15,  We  =  1000,  A  =  /3  =  sin,  f2  =  cos,  j  =  1,  =  0.025,  t)  =  v*  =  u*  = 

w*  =  0]. 


case 

n 

Ttr 

^34 

A; 

*6, (ml 

0 

i 

2tt 

0 

1 

120 

100 

1 

i 

2tt 

0.1 

1 

100 

80.1  (0.2) 

2 

3 

2t r 

0.1 

1 

117.9 

98 

3 

3 

2tt/3 

0.1 

1 

122.1 

102 

4 

3 

20ir/3 

0.1 

1 

117.6 

97.5 

5 

3 

27T 

0.1 

10 

114.6 

94.8 
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Figure  9: 


Figure  10: 


Figure  Captions 

Schematic  depictions  of  the  investigated  semi-infinite  annular  and  conical  swirling  sheet  con¬ 
figurations:  (a)  dilational  modulation,  (b)  sinuous  modulation. 

Solutions  1-3  to  Eq.  (20)  obtained  via  Mathematica™. 

Solutions  4  and  5  to  Eq.  (20)  obtained  via  Mathematical^. 

Solutions  6  and  7  to  Eq.  (20)  obtained  via  Mathematica™. 

Solutions  to  Eq.  (20)  for  wavenumbers  ki(i  =  1, . . . ,  7)  [We  =  5,  ft  =  10,  n  —  3].  For  u  <  0.7, 
fci_4  are  predominantly  dilational  mode  waves  and  £5-7  are  predominantly  sinuous  mode 
waves. 

Growth  rates  for  unstable  dilational  mode  wave  £3  [We  =  1000,  ft  =  10]:  o  :  n  =  0,  □  :  n  =  1, 
O  :  n  =  2,  A  :  n  =  3,  V  :  n  =  4. 

Growth  rates  for  unstable  dilational  mode  wave  k$  [We  =  50,  R  =  10]:  o  :  n  =  0,  □  :  n  =  1, 
0:n  =  2,  A:n  =  3,  V:ti  =  4. 

Growth  rates  for  unstable  dilational  mode  wave  k^  [We  =  5,  R  =  10]:  o  :  n  —  0,  □  :  n  =  1, 
O  :  n  =  2,  A  :  n  =  3,  V  :  n  =  4. 

Wave  angle  dd/dz*  =  k/n  for  unstable  dilational  mode  wave  k$  [We  =  50,  R  =  10]:  □  :  n  =  1, 
O  :  n  =  2,  A  :  n  =  3,  V  :  n  —  4. 

Instantaneous  film  thickness  distribution  for  swirling  annular  axisymmetric  liquid  film  (n  =  0) 
with  modulation  of  the  film  thickness  at  the  nozzle  exit  in  order  to  generate  the  unstable 
dilational  mode  wave  predicted  by  linear  theory  [ft  =  10,  u  =  2tt Tp  =  0.26,  We  =  1000, 
AAr*  =  0.35;  i;  =  390,  tmb  =  422.85]. 
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Figure  11:  Instantaneous  three-dimensional  film  thickness  distribution  for  swirling  annular  liquid  film 
(n  =  3)  with  modulation  of  the  film  thickness  at  the  nozzle  exit  in  order  to  generate  the 
unstable  dilational  mode  wave  predicted  by  linear  theory  [i?  =  10,  u  =  2irTp  =  0.26,  We  = 
1000,  AAr-  =  0.35;  l’b  =  245,  t J  =  275.55]. 

Figure  12:  Instantaneous  film  three-dimensional  thickness  distribution  for  swirling  annular  liquid  film 
(n  =  3)  with  modulation  of  the  film  thickness  at  the  nozzle  exit  in  order  to  generate  the 
unstable  dilational  mode  wave  predicted  by  linear  theory  [R  =  10,  u  =  2ttTp  =  0.26,  We  =  50, 
AAr*  =  0.15]. 

Figure  13:  Instantaneous  film  interface  locations  for  the  swirling  annular  liquid  film  of  Fig.  12.  Film 
break-up  length  and  time  are  1%  =  95  and  tb  =  115.2,  respectively. 

Figure  14:  Instantaneous  film  thickness  distribution  for  swirling  annular  axisymmetric  liquid  film  (n  =  0) 
with  modulation  of  the  film  thickness  at  the  nozzle  exit  in  order  to  generate  the  unstable 
dilational  mode  wave  predicted  by  linear  theory  [/?  =  10,  uj  =  2xTp  =  0.06,  We  =  50, 
AAr-  =  0.15;  tj  >  225]. 

Figure  15:  Instantaneous  three-dimensional  film  thickness  distribution  for  swirling  annular  liquid  film 
(n  =  3)  with  modulation  of  the  film  thickness  at  the  nozzle  exit  in  order  to  generate  the 
unstable  dilational  mode  wave  predicted  by  linear  theory  [72  =  10,  oj  =  2ttTp  =  0.06,  We  =  50, 
AAr-  =  0.15;  i;  =  117.5,  ij  =  127.2]. 

Figure  16:  Slices  through  the  distorting  liquid  film  of  Fig.  15  at  various  downstream  locations.  Nozzle 
exit  located  at  r*  =  0. 

Figure  17:  Instantaneous  film-interface  locations  for  non-swirling  (pressure-stabilized)  annular  liquid 
film  dilationally  modulated  at  the  nozzle  exit  z“  =  0  according  to  Eqns.  (30)  and  (31) 
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[R  =  io,  We  =  1000,  Tp  =  Te  =  10,  n  =  5,  A  =  B  =  0.04]. 

Figure  18:  Instantaneous  film-interface  locations  for  non-swirling  (pressure-stabilized)  annular  liquid 
film  sinusoidally  modulated  at  the  nozzle  exit  z*  =  0  according  to  Eq.  (32)  [R  —  10,  We  = 
1000,  Tp  =  Te  =  10,  n  =  5,  C  =  0.1]. 

Figure  19:  Instantaneous  non-dimensional  film  thickness  distribution  for  the  non-swirling  (pressure- 
stabilized)  annular  liquid  film  of  Fig.  18,  subject  to  three-dimensional  sinuous  forcing  at 
the  nozzle  exit  z*  =  0  according  to  Eq.  (32)  [R  —  10,  We  =  1000,  Tv  =  Te  =  10,  n  =  5, 
C  =  0.1]. 

Figure  20:  Instantaneous  non-dimensional  film  thickness  distribution  for  swirling  annular  liquid  film 
subject  to  three-dimensional  sinuous  forcing  at  the  nozzle  exit  z*  =  0  according  to  Eq.  (32) 
[R  =  10,  We  =  1000,  Tp  =  Te  =  10,  n  =  5,  C  =  0.1]. 

Figure  21:  Instantaneous  outer  interface  topology  for  swirling  conical  liquid  film  with  three-dimensional 
sinuous  modulation  enforced  at  the  nozzle  exit  z*  =  0  according  to  Eq.  (32)  [iZ  =  10, 
We  =  1000,  Tp  =  Te  =  10,  n  =  5,  C  =  0.1,  swirl  number  k  —  10]. 

Figure  22:  Instantaneous  non-dimensional  film  thickness  distribution  for  swirling  conical  liquid  film  of 
Fig.  21. 

Figure  23:  Disturbed  and  undisturbed  outer  film  interface  for  the  swirling  conical  liquid  film  of  Fig.  21. 

Figure  24:  Instantaneous  film-interface  locations  for  swirling  conical  liquid  film  with  three-dimensional 
sinuous  modulation  enforced  at  the  nozzle  exit  z*  =  0  according  to  Eq.  (32)  [R  =  10, 
We  =s  1000,  Tp  =  Te  =  10,  n  =  5,  C  =  0.25,  swirl  number  k  =  10]. 

Figure  25:  Sketch  of  general  film  break-up  pattern  for  swirling  annular  liquid  film  in  (z,  9)  plane  at 
instant  in  time  with  one  dilational  or  sinuous  travelling  wave  modulated  at  the  nozzle.  Solid 
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lines:  maximum  thickness  contour  lines  (local  line  thickness  =  measure  of  maximum  local 
film  thickness).  Dashed  lines:  film  perforation  pattern. 

Figure  26:  Sketch  of  general  film  break-up  pattern  for  swirling  annular  liquid  film  in  ( z,9 )  plane  at 
instant  in  time  with  standing  dilational  circumferential  wave  modulated  at  the  nozzle.  Solid 
line:  maximum  thickness  contour  line  for  single  oblique  wave.  Dashed  line:  instantaneous 
free  film  edge.  Shaded  area:  void.  Solid  circle:  local  film-thickness  maximum. 


Figure  27:  Sketch  of  general  film  break-up  pattern  for  swirling  conical  liquid  film  in  (z,  9)  plane  at  instant 

$ 

in  time  with  standing  dilational  circumferential  wave  modulated  at  the  nozzle.  Solid  line: 
maximum  thickness  contour  line  for  single  oblique  wave.  Dashed  line:  instantaneous  free  film 
edge.  Shaded  area:  void.  Solid  circle:  local  film- thickness  maximum. 

Figure  28:  Sketch  of  general  film  break-up  pattern  for  swirling  annular  liquid  film  in  (z,  9)  plane  at 
instant  in  time  with  standing  sinuous  circumferential  wave  modulated  at  the  nozzle.  Solid 
line:  maximum  thickness  contour  line  for  single  oblique  wave.  Dashed  line:  instantaneous 
free  film  edge.  Solid  circle:  local  film-thickness  maximum. 

Figure  29:  Sketch  of  general  film  break-up  pattern  for  swirling  conical  liquid  film  in  ( z,0 )  plane  at 
instant  in  time  with  standing  sinuous  circumferential  wave  modulated  at  the  nozzle.  Solid 
line:  maximum  thickness  contour  line  for  single  oblique  wave.  Dashed  line:  instantaneous 
free  film  edge.  Solid  circle:  local  film-thickness  maximum. 
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swirling  conical  sheet 


swirling  conical  sheet 
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Solutions  to  Eq.(20) 


Figure  3: 


47 


Solutions  to  Eq.(20) 
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Solutions  to  Eq.(20) 


Figure  5: 
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Dynamic  Stretching  of  A  Planar  Liquid  Bridge 

C.  Mehringa)  and  W.  A.  Sirignano0^ 

Department  of  Mechanical  and  Aerospace  Engineering ,  University  of  California ,  Irvine 

(February  18,  2003) 

Abstract 

A  thin  incompressible  viscous  planar  free  liquid  film  in  a  void  and  under  zero 
gravity  is  analyzed  by  means  of  a  reduced-dimension  (lubrication)  approach. 
Linear  analysis  focuses  on  films  with  harmonic  modulations  in  the  axial  film 
velocity  enforced  at  the  ends  of  the  planar  bridge.  Effect  of  changes  in  the 
problem  parameters  on  the  overall  distortion  characteristics  of  the  film  are 
discussed.  Nonlinear  film  distortion  and  break-up  is  investigated  for  the  case 
of  temporally  increasing  velocity  at  the  end  of  the  film  resulting  in  continuous 
film  stretching  eventually  leading  to  film  rupture.  Implementation  of  the 
employed  numerical  model  is  validated  for  the  linear  limit  by  comparison 
with  the  analytical  linear  solutions  and  for  harmonically  modulated  film-end 
velocities. 

Within  the  nonlinear  analysis  of  the  continuously  stretched  film  bridge, 
several  distinct  film  topologies  are  identified  depending  on  liquid  Weber  num¬ 
ber  and  Reynolds  number,  i.e.,  the  magnitude  of  the  stretching  rate  (end 
velocity)  compared  to  signal  propagation  rates  through  the  liquid  via  capil¬ 
lary  waves  and  viscous  action.  That  is,  the  Weber  number  is  the  square  of 
the  ratio  of  stretching  rate  to  capillary  wave  velocity  while  Reynolds  number 
is  the  ratio  of  stretching  rate  to  the  characteristic  viscous  velocity.  Here,  film 
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topology  is  typically  characterized  by  three  distinct  regions,  i.e.  a  film  wedge 
forming  at  the  pulling  end(s),  the  film  center  region  and  a  transition  region. 
The  size  and  shape  of  these  regions  greatly  depends  on  the  particular  case 
under  investigation. 

Film  distortion  characteristics  observed  for  continuously  compressed  pla¬ 
nar  films  conform  with  observations  made  by  other  authors  for  the  similar 
case  of  contracting  free  liquid  films. 


Typeset  using  REVTgX 
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I.  INTRODUCTION 


Written  reports  on  the  analysis  of  thin  films  originate  as  far  back  as  Leonardo  da  Vinci1 
who  described  the  formation  of  a  liquid  film  when  pulling  a  reed  out  of  soapy  water.  Today, 
liquid  films  are  found  in  a  wide  variety  of  engineering  applications  including  novel  concepts 
such  as  liquid  film  space  radiators,2,3  liquid  lithium  films  for  inertial  fusion  reactors,4  liquid 
bells  employed  as  toxic  waste  incinerators5  and  modulated  injected  fuel  films  for  actively 
controlled  fuel  spray  formation  in  the  combustion  chambers  of  advanced  gas  turbine  engines.6 

Previous  work  by  the  authors  has  focused  on  the  initial  distortion  and  disintegration 
of  continuously  discharging  films  (or  periodically  disturbed  infinite  films)  with  relevance  to 
the  described  actively  controlled  fuel  atomization  process.7-10  In  an  effort  to  understand  the 
dynamics  of  the  larger  fluid  ligaments  initially  shed  from  the  discharging  continuous  film,  a 
finite  length  continuously  stretched  or  compressed  two-dimensional  planar  liquid  film  in  a 
void  (also  referred  to  as  liquid  bridge)  has  been  analyzed.  The  avoidance  of  considering  a 
finite-length  film  with  free  ends  allows  the  use  of  lubrication  equations  (previously  derived  by 
the  authors)  without  additional  modelling  needs  at  the  film  ends.  Those  equations  become 
invalid  for  a  freely  contracting  or  stretching  film  whose  ends  are  characterized  by  a  radius 
of  curvature  of  the  film  interface  that  becomes  as  small  as  the  film  thickness. 

Thin  planar  films  in  a  surrounding  void  and  subject  to  two-dimensional  symmetric  (or 
antisymmetric)  disturbances  are  well  known  to  be  linearly  stable.11  However,  nonlinear  the¬ 
ory  has  shown  that  the  superposition  of  multiple  symmetric  disturbances  on  a  periodically 
disturbed  infinite  film  can  produce  nonlinear  film  instability.12,13  Furthermore,  for  very  thin 
films,  long-range  molecular  forces  (acting  between  the  two  film  interfaces)  can  also  result  in 
film  instability,  even  by  linear  theory.  There  have  been  various  linear  and  nonlinear  analyses 
on  this  subject.  For  a  detailed  review  the  reader  is  referred  to  Ref.  11  and  14. 

Previous  work  focused  on  periodically  disturbed  films  with  zero  or  constant  “base  flow” 
velocity.  Previous  analyses  of  finite  length  liquid  bridges  focused  on  cylindrical  fluid  columns 
or  annular  liquid/soap  films  rather  than  two-dimensional  films.15’16  The  proposed  analysis  a 
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dynamically  stretched  or  compressed  planar  liquid  bridge  has  not,  to  the  knowledge  of  the 
authors,  been  reported  in  the  literature. 

A  brief  nonlinear  analysis  of  a  contracting  two-dimensional  free  liquid  film  was  presented 
in  Ref.  17.  The  authors  described  topological  changes  in  the  interface  distortion  with  changes 
in  Reynolds  number  and  Weber  number.  The  present  work  focuses  on  the  interface  dynamics 
of  the  stretching  film  for  various  Weber  numbers,  Reynolds  numbers,  film-thickness-to-length 
ratios  and  the  acceleration  period  (time  over  which  the  film  ends  are  accelerated  from  zero 
velocity  to  the  maximum  pull  velocity). 

In  this  context,  it  is  worthwhile  to  note  that  in  the  spray  formation  application  described 
earlier,  stretching  of  the  initially  formed  liquid  ligaments  proceeds  the  contraction  process, 
the  latter  taking  place  once  a  sufficiently  stretched  ligament  has  ruptured. 


II.  GOVERNING  EQUATIONS  AND  PROBLEM  SET-UP 


The  governing  equations  within  the  liquid  phase  of  an  incompressible  viscous  planar  free 
liquid  film  at  zero  gravity  are  given  by: 


du 

dx 


dw 

+  a7  =  0 


(i) 


du  du  du 

m+udi  +  wTz 


1  d(p  +  $)  /  d2u  d2u\ 

p  dx  y&r2  +  dz2) 


(2) 


dw  dw 

W  +  nfc 


dw  1  d(p  +  $) 

+  w—  = - - - 

dz  p  dz 


+  v 


d2w  d2w\ 

dx2  dz2  J 


(3) 


where  x  is  the  direction  along  the  film,  z  is  the  transverse  direction,  u  and  w  are  the 
respective  velocity  components,  p  is  pressure,  v  is  kinematic  viscosity,  and  $  denotes  the 
potential  energy  function  per  unit  volume  which  accounts  for  the  presence  of  long-range 
molecular  forces  and  which  depends  on  the  thickness  of  the  free  film.18’19  When  the  film 
thickness  is  large,  the  effect  of  the  Van  der  Waals  (long-range  molecular)  forces  is  well 
represented  through  the  Navier-Stokes  equations  with  the  use  of  surface  tension  in  the 
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interface  conditions  and  the  <h-term  is  negligible.  However,  $  becomes  important  and  can 
even  dominate  film  dynamics  as  the  film  thickness  drops  below  100-1000  Angstroms.14  The 
present  analysis  will  focus  on  film  dynamics  at  fairly  large  film  thicknesses  so  that  here,  Van 
der  Waals  forces  are  included  in  the  governing  equations  only  to  illustrate  their  potential 
relevance  as  film  thinning  progresses,  i.e.  before  film  rupture  occurs.  The  majority  of  results 
presented  here  neglect  the  effects  through  the  $-term  of  the  described  intermolecular  forces; 
rather,  their  integrated  effect  appears  in  the  continuum  representation. 

Boundary  conditions  at  the  film  interfaces  z  =  h(x,  t)±  are  given  by  the  following  kine¬ 
matic  and  dynamic  conditions. 

Kinematic  condition: 


w  — 


dh±  dh± 
dt  ^U±  dx 


(4) 


Dynamic  condition: 


T±  •  n-t  =  —K±an±  (5) 

where  +  and  —  refer  to  the  upper  and  lower  interface,  respectively.  In  Eq.  (5),  T  denotes  the 
stress  tensor,  k  represents  the  interface  curvature  and  n  is  the  outward  pointing  unit  normal 
vector  at  the  interface.  Eq.  (5)  has  a  tangential  and  a  normal  component.  The  tangential 
component,  i.e.  the  total  shear  stress  on  the  interface,  is  identical  to  zero  since  externally 
imposed  surface  forces,  as  well  as  surface  gradients  of  surface  tension  a  have  been  neglected 
here. 

The  stress  tensor  is  given  by 

-P  + 

f=  +  (6) 

; p  +  2  p£ 

and 
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(7) 


dh±/dx ,  ±1) 
yjl  +  Jdh^Jdx)2 

_v  a  i  d-^l^ 

[i  +  (dh±/dx)2\w 


(8) 


Accordingly,  the  components  of  the  dynamic  condition  Eq.  (5)  in  the  axial  and  transverse 
directions  can  be  written  as: 


1 

CN 

-H 

1 

1 _ 

du±  dw± 

- —  +  -7~ 

\dx)  _ 

dz  dx 

-  4 


-P±  + 


2  p 

Cl± 


e2±- 


dw±  dh±  ( du± 
dz  dx  1  dz 


+ 


dh±  du±  _ 

(9) 

dx  dx 

dw±\ 

dx)\=aK± 

(10) 

with  ei±  =  [1  +  ( dh±/dx )2]  and  e2±  =  [1  -  ( dh±/dx )2].  Since  the  film  thickness  is  small 
compared  to  the  length  of  the  film,  we  employ  a  polynomial  expansion  of  the  dependent 
variables  6  =  u,  w,p,  or  $  in  terms  of  z,  i.e. 


6  =  90(x,  t )  +  9i(x,  t)z  +  02(x,  t)z2  +  ... 


(11) 


analogous  to  the  one  presented  in  Ref.  [5]  for  inviscid  films. 

Considering  zero-  and  first-order  terms  in  the  series  expansion  for  u  and  w  within  the 
kinematic  condition  (4)  yields, 


dh+ 

w*-~§r+u° 

a/i_ 

w-  =  ~dT  +  u  o 


dh+ 

dx 
dh _ 


dx 


+  U\ 


U'2 


(12) 


and 


<5 

W+=Wq+W  1- 

w.  =  w0  -  wi  ^  (13) 

where  5  denotes  the  film  thickness  6(x,t)  =  h+  —  h-  and  h±  =  ±6/2,  respectively.  The 
summation  +  w _  and  subtraction  w+  —  w _  then  provides 
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Note  that,  the  forces  driving  the  transverse  motion  of  the  fluid  are  of  first  order  (i.e.  pi  and 
Therefore,  the  transverse  fluid  motion  in  the  considered  dilational  case,  the  transverse 
fluid  motion  can  be  neglected  in  this  analysis.  Nevertheless,  substituting  w0  from  Eqs.  (14) 
into  Eq.  (18),  one  obtains  an  equation  relating  the  first-order  expansion  coefficients  u\,pi 
and  $i. 

Eqs.  (16)  and  (17)  represent  the  leading-order  equations  considered  within  the  present 
analysis.  In  order  to  close  the  lowest-order  system  of  equations  and  proceed  with  its  (nu¬ 
merical)  solution,  p0,  $0  and  u?  remain  to  be  expressed  in  terms  of  known  quantities  or 
quantities  for  which  governing  equations  are  given,  i.e.  8  and  u0. 
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An  approximation  for  rt2  is  obtained  by  employing  the  second-degree  polynomial  expan¬ 
sion  of  u  in  the  evaluation  of  d2u[dz2.  This  yields  d2u/dz*  —  2u2  or  after  integration  from 
the  lower  to  the  upper  interface  boundaries, 

1 


U2  =  25 


' du 
,  dz , 


'du' 
,  dz , 


(19) 


The  right-hand-side  of  Eq.  (19)  can  be  expressed  by  means  of  the  shear-stress  boundary 
condition  Eq.  (9)  solved  for  ( du/dz)± ,  i.e. 


'du\  _  ( dw\ 


\®z) 


dx) 


■1  +4 
± 


-l 


dh±  du± 
dx  dx 


(20) 


where  an  expression  for  (du>/dx)±  is  obtained  by  taking  the  derivative  of  Eq.  (1)  in 
x— direction  and  subsequent  integration  in  z—  direction.  The  result  of  this  procedure  yields 


'dw\ 

,dx). 


'  dw 
,  dx 


)  i =rf^d „-f  / 

J  _J  Jh-  dx 2  dx  Jh 


d  fh+  du  un+ 
dx  Jh .  dx  Z  dx  dx 


dh+  du+  dh _  du- 

+ 


d_ 

dx 

d 


'  d_  rh+ 

dx  Jh. 


u  d  .j.-,  _d5 

=  Tx  Tx{uS)~uTx 


dh+  dh- 

udz  — — u+  +  u- 

dx  dx 

d5du  _  d_ 

dx  dx  dx 


dx  dx 
d5_du 
dx  dx 


.du 
1  dx 


dSdu  _  ;d2u 
dxdx  dx 2  1  ' 


where  u  —  J^+  udz/5  —  uq  +  <52/2u2  +  h.o.t.  Introducing  Eq.  (21)  into  Eq.  (20)  yields  the 
integral  form  of  the  shear  stress  boundary  condition, 


'du 
.  dz 


-  ©  J 


rd2uo  .d5  du0 


_fi«V 

\2  dx  J 


+  0^ 1> 


(22) 


if  only  the  two  lowest-order  terms  in  the  expansion  of  u  are  considered.  Therefore, 


ld2n0  ,  2  d5  du0  ^(x2  x^  \ 
U2  =  2  -d^+6di-^  +  0{S'6di] 


(23) 


The  lowest-order  approximation  for  the  interface  pressures  p±  is  obtained  by  considering 
only  terms  up  to  order  one  in  the  series  expansion  for  p,  i.e. 


Po  ±Pi- 


(24) 


This  yields  po  =  (p+  +p_)/ 2  whereby  p+  and  p_  are  taken  from  the  normal  stress 
boundary  condition  Eq.  (10).  The  result  is 
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2  dx 2 


.  .  on  -3/2 

,  1  d6\2  n  du0  ,  1  /0<J 

1  +  j(&  1  +  2l& 


,-Ip 

4  V0x 


+  0(<52,5^) 


where  again  only  the  two  lowest-order  terms  in  the  approximation  to  u  have  been  considered 
in  Eq.  (10). 

Note  that,  when  evaluating  tt2  and  p0  in  Eq.  (17),  the  employed  lower-order  approxima¬ 
tion  to  the  normal  stress  and  shear  stress  boundary  conditions  Eqs.  (9)  and  (10)  have  not 
employed  any  simplifications  with  respect  to  the  interface  geometry;  in  particular,  the  full 
curvature  term  is  retained  in  the  evaluation  of  the  surface  tension  within  Eq.  (10). 

Using  Eqs.  (22)  and  Eq.  (25),  Eqs.  (16)  and  (17)  can  be  rewritten  to  obtain  the  leading- 
order  equations  for  symmetric  film  distortion  governing  the  film  thickness  6  =  (h+  —  h_) 


and  the  axial  velocity  Uq 


8(u5) 
dt  dx 


du  du  d2u  ,t  ■. ,  .  1 85  du 

^  +  ^  =  ^[4  +  2(/-l),]  +  4,---p 


dx2  1  '  w 

.  du  86  d25 


n  du  dS  d25  2  cr  838  ,_3/2 

4~2u - u  H - f 

dx  dx  dx 2  2  p  dx 3 

_3 ad5_  (&8\2  5/2  3 Ad6_ 

8  p  dx  l  Ox2 )  p8 4  dx 


where  the  subscript  ‘0’  in  u0  has  been  omitted  and 


.  ,  fidsy  ,  (ids 

+  V2dx)  ,9~  \2dx 


Furthermore,  the  last  term  in  Eq.  (27)  represents  —  (d$o/dx)/p  where  $o  has  been  expressed 


—  .<4  6 


according  to  Ref.  18,  with  A  =  A°/(6tt)  where  A°  is  the  Hamaker  constant  ^4°  =  10  13  erg, 
but  without  an  additive  constant.20,21  If  the  denominators  in  Eqs.  (7 , 8)  are  approximated 
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by  1  (i.e.  considering  ( dh±/dx )2  <C  1),  then  Eqs.  (26 , 27)  reduce  to  the  equations  presented 
by  Erneux  and  Davis  in  Ref.  21  for  long  wavelength  disturbances  on  thin  films,  i.e. 


96  d(u5)  _ 

dt  dx 

du  9u_  d2u  v  35  du  a  3*6  3A  35 
dt  U  dx  V  dx2  ^  5  dx  dx  2p  dx 3  p5 4  dx 


(29) 


Here  we  have  a  parameter  4  in  the  viscous  term  to  which  viscous  forces  inside  the  liquid 
contribute  one,  shear  stress  at  the  interface  contributes  one,  and  normal  stress  contributes 
two.  The  second  term  on  the  right-hand-side  is  a  viscosity  correction  term  which  results 
from  the  shear  stress  at  the  interface. 

In  the  present  analysis,  no  approximations  are  made  in  the  formulation  of  the  interface 
conditions.  In  particular,  the  full  curvature  term  is  retained  in  Eq.  (10).  This  “ad-hoc” 
assumption  in  an  otherwise  lower-order  one-dimensional  model  follows  the  analyses  and 
observations  by  Ruckenstein  and  Jain,18  Eggers22  and  Eggers  and  Dupont,23  who  studied  the 
dynamics  of  drop  formation  from  cylindrical  liquid  jets  or  columns.  Similar  considerations 
were  made  by  Mehring  and  Sirignano7,8  who  analyzed  capillary  waves  on  thin  planar  and 
annular  free  liquid  films. 

Eqs.  (26)  and  (27)  have  been  used  to  study  the  dynamic  stretching  of  a  thin  liquid  film 
subject  to  equal  but  opposite  pulling  velocities  at  both  ends.  Due  to  the  symmetry,  it 
is  sufficient  to  analyze  only  one  half  of  the  film.  The  configuration  considered  within  the 
present  work  is  shown  in  Fig.  1,  with  the  plane  of  symmetry  located  at  x  =  0.  The  film  is 
assumed  to  be  initially  undisturbed  with  5(x,  0)  =  5{  and  u(x,  0)  =  0  for  0  <  x  <  l{  where 
(2/j)  is  the  initial  total  length  of  the  film.  Boundary  conditions  in  the  plane  of  symmetry 
at  x  =  0  for  t  >  0  are  given  by  ( 35/dx )  =  0  and  it  =  0.  The  pulling  velocity  up(t)  at 
x  =  l(t)  =  Jq  up(t')  dt1  +  was  assumed  to  increase  smoothly  from  0  to  a  maximum  steady 
value  U  according  to 


u[l(t),t]  =  up(t)  =  £7/2  [1  —  cos  (27Tf/Tp)]  for  0  <  t  <  Tv/2  and 

up(t)  =  U  for  t  >  Tp/2  (30) 
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The  film  thickness  S  at  the  pulling  end  was  assumed  to  be  fixed,  i.e.  <5[/(f),t]  =  <5,. 

With  the  prescribed  initial  and  boundary  conditions,  Eqs.  (26)  and  (27)  have  been 
solved  numerically  after  nondimensionalization  and  after  introducing  the  following  coor¬ 
dinate  transformation: 


V  =  W)  '  T  =  t  (31) 

where  2l(t)  denotes  the  instantaneous  total  film  length. 

With  this  change  of  variables,  the  problem  becomes  a  fixed-boundary-value  problem, 
which  greatly  simplifies  its  numerical  analysis.  Note  that,  with  increasing  film  length  l(t), 
the  fixed  spatial  resolution  in  77  used  within  the  subsequent  numerical  analysis,  corresponds 
to  a  continuously  decreasing  resolution  in  x,  i.e.  physical  space.  Therefore,  predictions  for 
film  lengths  at  breakup  and  time  until  breakup  have  to  be  considered  with  some  caution.  It 
is  noted  here  that,  the  present  work  is  focused  on  a  description  of  film  dynamics  in  the  course 
of  film  stretching  (or  compression)  prior  to  film  rupture  and  not  on  an  exact  prediction  of 
film  length  at  rupture  or  time  until  film  rupture. 

For  the  non-dimensionalization,  the  following  non-dimensional  variables  have  been  em¬ 


ployed, 


=  i  u*  =  -  r  =  -  r*  =  —  C  -  — 
Si  ’  U'  u'  k/uc  ’  C  k 


where  6i,  li  and  U  denote  the  initial  film  thickness,  initial  film  length  and  maximum  pull 
velocity,  as  introduced  earlier.  The  capillary  velocity  uc  is  given  by  uc  =  ^/cr/ (/)£,).  Intro¬ 
ducing  the  prescribed  transformation  and  non-dimensionalization,  Eqs.  (26)  and  (27)  are 
given  by: 


dS*  RTT  u‘  ~  Wp  dS*  r—5'  du *  n 

v—  +  vWe- - - - - — h  VWe—  = 

3r*  l*  drj  l*  dr] 


du *  \u*-r]u'du * 

—  +  VWe± - - - ——  = 

dr*  l*  dr] 
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3A'  dSm 


2?  _ 

+  /*4  ^  drj  dp  drj1  +  S‘ 4  lm\/We  dp 

,  C2  t*  —3/2  W  _  3C4  5/2  ^  ,W\2 

2l*3VWeJ  dp3  8 l‘*y/WeJ  dp  \  dp2  j 


with 

\2l *  y 

and  Ohnesorge  number,  Weber  number  and  nondimensional  ‘modified’  Hamaker  constant 
A*  given  by 


r  =  i  + 


'  C  dS'' 
2 /•  dr? 


0  = 


Oh  = 


Vfe 

Re 


where  the  Reynolds  number  Re  is  defined  a s  Re  =  U  Si/ v. 


(35) 


Linear  Analysis 


We  now  consider  a  linear  analysis  of  the  described  planar  liquid  film  with  its  ends  har¬ 
monically  modulated  180°  out  of  phase  keeping  the  symmetry  .  The  average  mean  axial 
velocity  of  the  film  is  zero  in  this  case  and  the  average  half-length  is  L.  Therefore,  here  the 
capillary  velocity  uc  =  \jojpS i  has  been  employed  in  order  to  nondimensionalize  Eqs.  (26) 
and  (27). 

Replacing  variables  S',  u*  and  l*  (now  nondimensionalized  using  uc  and  Si)  by  their  mean 
and  fluctuating  values,  i.e. 

S'  =  1  +  S'  ,  u  =  u'  ,  r  =  L  +  l' 


linearization  of  Eqs.  (33)  and  (34)  yields 

dS/_  l_dv/_ 
dr”  +  L  dr] 

du'  4  1  d2u'  n  Al  1  1  dS'.  1  1  d35' 
dr  Re  I?  dr, 2  ■  Re2  I.  dr,  21?  drf 


(36) 

(37) 
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with  A'  =  A/(pv26i)  and  now  77  =  x/L.  Note  that,  in  contrast  to  the  nonlinear  analysis 
of  the  continuously  stretched  or  compressed  liquid  bridge,  the  Reynolds  number  within  this 
linear  analysis  is  based  on  the  initial  (undisturbed)  film  thickness  6i  and  the  capillary  velocity 
uc  =  \JaJpdi.  A  modal  analysis  of  the  above  equations  is  conducted  by  assuming  wavelike 


solutions  of  the  form 


5' =  a- ei{uT'~kT,) 
v!  =  b-ei(uir'-kn) 


where  to  is  the  non-dimensional  oscillation  frequency  at  a  particular  location  along  the  film 
and  k  =  27r<5j/A  is  the  wavenumber  characterizing  the  disturbances  appearing  on  the  film. 

Introducing  the  above  equations  for  5'  and  u'  into  the  linear  equations  Eq.  (36)  and 
Eq.  (37)  yields  the  characteristic  equation  which  determines  the  dependence  of  the  distur¬ 
bance  wavenumbers  k  on  the  forcing  frequency  u. 

-  (41) 

with  K  —  k/L  and  u  =  JoJ{p5f)  Q  where  Q  denotes  the  dimensional  forcing  frequency. 

In  order  to  clearly  illustrate  the  effects  of  viscosity,  surface  tension  and  intermolecu- 
lar  forces  on  the  stability  of  the  liquid  film,  Eq.  (41)  is  solved  for  forcing  frequency  u>  in 
dependence  of  real  wavenumber  k,  i.e. 


oil  2  =  -  <  i  ~~k2  ±i  +  12 A*  k2  —  2  «4 

’  2  Re  Re2  ' — v — '  ' — ■* — " 

v - - *  v - ^  molecular  surface  tension 


where  .4*  =  A1 1 Re2  —  A/(o62)  as  already  introduced  earlier.  In  Eq.  (42),  the  first  term  in 
curly  brackets  and  the  first  term  in  square  brackets  account  for  viscous  effects.  Note  that, 
the  combined  effect  of  both  terms  will  always  act  stabilizing  on  the  distorting  film,  i.e.  result 
in  real  values  for  u.  The  second  and  third  terms  on  the  right-hand-side  of  Eq.  (42)  account 
for  intermolecular  forces  and  surface  tension,  respectively.  While  the  former  might  result  in 
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film  instability  due  to  the  production  of  a  complex  w-value  with  negative  imaginary  part; 
the  latter,  by  itself,  will  always  produce  real  oz-values  indicative  of  a  stable  film.  While  the 
surface  tension  represents  the  effect  of  intermolecular  forces  for  larger  thickness,  it  acts  in 
opposite  manner  to  the  $-term  which  is  intended  to  represent  intermolecular  forces  in  the 
very  thin  film  situation.  That  is,  when  the  intermolecular  forces  act  primarily  on  neighboring 
molecules  in  the  same  surface  layer  as  happens  for  thicker  films,  the  effect  is  stabilizing  and 
is  well  represented  by  the  surface  tension  formulation.  However,  when  the  intermolecular 
forces  act  significantly  upon  molecules  in  the  other  surface  layer  of  a  very  thin  film,  the 
effect  is  destabilizing  and  is  represented  through  the  $-term. 

Eq.  (42)  show  that,  for  k 2  <  6A*,  u>  will  take  only  positive  or  negative  imaginary  values 
lji  and  u>2-  The  latter  will  result  in  an  exponentially  growing  solution  while  the  former  will 
produce  a  solution  which  decays  exponentially  in  time.  Note  that  the  unstable  growth  rate 
u)2  will  be  smaller  than  the  decay  rate  u>i.  For  6 A*  <  k2  <  6A*  /(I  —  S/Re2),  the  square-root 
term  in  Eq.  (42)  will  be  positive  but  smaller  than  4 k2 /Re.  In  this  case,  Eq.  (42)  will  produce 
solely  exponentially  decaying  solutions.  If  k2  >  6.4* /(I  —  S/Re2),  the  square-root  term  in 
Eq.  (42)  will  be  complex  and  the  resulting  frequencies  uii#  will  contain  an  exponentially 
decaying  part  (due  to  viscous  damping)  and  a  travelling-wave  part. 

If  the  effect  of  intermolecular  forces  is  neglected  in  the  analysis,  no  film  instability  exists 
and  u>  has  a  double  root  at  Re  =  2^2,  i.e.  u  1,2  =  2iK2/Re  with  exponentially  decaying 
solutions  only.  Exponentially  decaying  solutions  are  also  found  for  Re  <  2y/2,  however 
here  the  decay  rates  associated  with  each  solution  of  oj  are  different.  For  Reynolds  numbers 
above  2\f2  each  solution  for  oj  is  associated  with  an  exponentially  decaying  solution  and  a 
travelling-wave-type  solution.  Consider  the  linear  boundary-value  problem  of  a  planar  liquid 
film  bridge  (i.e.,  a  free  liquid  film  bounded  at  both  ends  by  a  solid  vertical  surface)  with 
harmonic  forcing  of  its  ends  at  frequency  oj.  Here,  Eq.  (41)  has  to  be  solved  for  wavenumbers 
k  dependent  upon  oj.  The  prescribed  modulation  will  result  in  a  film  distortion  composed 
of  four  wavenumbers  since  Eqs.  (36)  and  (37)  form  a  fourth-order  system.  If  intermolecular 
forces  are  neglected,  kx  through  /c4  are  given  by 
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(43) 


^1,2  =  ki^/L  =  ±2  —iu>  ^Re  1  —  \J Re  2  —  1/8^ 
«3,4  =  kz^/L  =  ±2  —  iui  ^i?e-1  +  \j Re~2  —  1/8^ 


where  the  ‘  +  ’  or  ‘  -  ’  sign  relates  to  k\  and  k^  or  kz  and  k 4,  respectively.  Note  that  k{ 
(i  =  1  through  4)  are  always  complex.  No  film  instability  exists  if  intermolecular  forces 
can  be  neglected.  Using  the  four  solutions  obtained  for  wavenumbers  kt  in  (43),  the  general 
solution  to  the  described  linear  boundary-value  problem  is  given  by: 

S'  =  £  a,-  •  e*‘(WT*-fcM)  ,  ,  (44) 

i=l  *=1 

where,  according  to  Eq.  (36),  the  eight  coefficients  a{  and  6,-  are  linearly  dependent  on  each 
other  according  to  four  relations 


dj  k, 

bi  uj  L 


(45) 


The  remaining  four  unknown  coefficients  are  determined  by  four  boundary  conditions 
imposed  at  the  ends  of  the  liquid  bridge.  For  the  present  analysis,  these  boundary  conditions 
have  been  specified  as 


u'(t*,T)  =  ±L )  =  ±U0  cos  (27 r/r*) 

5'{t*,t]  =  ±L)  =  0  (46) 


The  prescribed  boundary-value  problem  has  been  solved  for  a  variety  of  Reynolds  numbers 
Re,  forcing  frequencies  /  =  and  thickness-to-length  ratios,  i.e.  L- values. 


III.  RESULTS 

A.  Linear  Analysis 

We  now  describe  the  analytical  results  obtained  for  the  linear  boundary-value  problem 
specified  in  the  previous  section.  Some  figures  also  contain  the  corresponding  numerical 
solutions  obtained  by  using  the  numerical  scheme  employed  to  solve  the  nonlinear  distortion 
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of  a  continuously  stretched  or  compressed  film  discussed  later.  The  excellent  agreement 
between  linear  analytical  and  numerical  solutions  illustrates  the  proper  implementation  of 
the  lubrication  equations  Eqs.  (26)  and  (27),  at  least  in  the  linear  limit. 

Effect  of  Changes  in  Reynolds  number. 

For  large  Reynolds  numbers  Re,  viscous  damping  of  the  capillary  waves  generated  by  the 
modulation  at  the  film  ends  is  very  small  and  the  instantaneous  envelope  in  the  film  thick¬ 
ness  distribution  is  essentially  uniform  along  the  film.  See  Fig.  2.  As  the  Reynolds  number 
decreases,  capillary  waves  are  subject  to  viscous  damping  and  the  maximum  disturbance 
amplitude  in  the  film  thickness  is  decreased  and  the  amplitude  of  the  instantaneous  en¬ 
velope  in  S'  along  the  film  decreases  towards  the  film  center;  see  Fig.  3.  As  the  Reynolds 
number  is  decreased  further,  the  decay  rate  of  the  described  envelope  towards  the  film  center 
increases.  At  low  enough  Re-values,  no  wave  phenomena  is  observed.  However,  for  the  case 
shown  in  Fig.  4,  information  on  the  film  modulation  still  reaches  the  film  center  resulting  in 
a  time- varying  local  extremum  in  S'.  Variations  in  S'  are  smallest  at  the  film  center.  At  very 
low  Reynolds  numbers,  information  propagation  is  dominated  by  viscous  effects;  see  Fig.  5. 
We  observed  that,  whereas  the  overall  maximum  disturbance  amplitude  in  S'  (at  a  particular 
time)  decreases,  its  amplitude  at  the  film  center  is  increased.  Spatial  variations  in  S'  are 
moderate  at  very  low  Reynolds  numbers,  approaching  uniform  instantaneous  film  thickness 
along  the  film,  with  exception  close  to  the  film  ends  where  S'  =  0  has  been  enforced.  For 
the  same  case,  the  axial  velocity  profile  along  the  film  is  nearly  linear  as  the  velocity  at  the 
film  ends  reaches  its  maximum  values  (positive  or  negative). 

Effect  of  Changes  in  Thickness-to- Length  Ratio  l/L: 

Comparison  amongst  Figs.  3,  4,  5  and  Figs.  6,  7,  8  illustrates  that,  the  characteristics  of  film 
topology  and  film  dynamics  (including  instantaneous  maximum  disturbance  amplitudes)  at 
the  film  ends  do  not  change  with  decreasing  film  thickness-to-length  ratio.  However,  at 
larger  L-values  the  decay  rate  of  the  envelope  in  S'  along  the  film  towards  the  film  center  is 
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increased  so  that  the  center  region  remains  essentially  undisturbed.  For  very  larger  L- values, 
decay  rates  of  the  envelope  in  the  film  thickness  disturbance  do  not  change  significantly  with 
changes  in  L.  This  can  be  observed  by  comparing  Fig.  7  and  9,  for  example. 

Effect  of  Changes  in  Forcing  Frequency  f: 

Changes  in  film  distortion  characteristics  with  increases  in  modulation  frequency  /  (at  con¬ 
stant  Re-  and  L- values),  such  as  those  observed  in  Figs.  4  and  10  or  Figs.  5  and  11,  resemble 
those  observed  for  the  similar  case  but  larger  L- value  if  results  are  plotted  over  an  axial 
coordinate  normalized  by  the  nondimensional  film  length  L.  See  Figs.  10  and  7  or  Figs.  11 
and  8.  Note  however  that  maximum  thickness  disturbances  are  significantly  smaller  for  the 
case  with  increased  forcing  frequency.  Also,  recall  that  maximum  disturbance  amplitudes 
did  not  notably  change  with  changes  in  the  thickness-to-length  ratio  (i.e.  1/L- value)  of  the 
film. 


B.  Nonlinear  Analysis 

The  nonlinear  analysis  has  focused  on  films  which  are  continuously  stretched  or  com¬ 
pressed  at  their  ends.  Numerical  solutions  of  the  nondimensional  nonlinear  evolution  equa¬ 
tions  were  obtained  explicitly  by  employing  central  differencing  for  the  spatial  derivatives  in 
combination  with  one-sided  spatial  differencing  for  the  evaluation  of  the  highest-order  spatial 
derivative  at  the  pulling  end.  Time  step  At*  and  mesh  size  Arj  were  fixed  throughout  each 
simulation  and  have  been  chosen  in  order  to  guarantee  stability  of  the  numerical  solution  and 
in  order  to  provide  an  accurate  prediction  of  film  topology  and  film  dynamics  in  the  course  of 
film  stretching.  The  latter  has  been  demonstrated  by  comparison  with  results  using  reduced 
time-steps  and/or  mesh  sizes.  It  is  noted  here  once  more  that,  predictions  of  film  break-up 
length  and  break-up  time  for  significantly  stretched  films  have  to  be  considered  with  some 
caution  due  to  the  deteriorating  spatial  resolution  resulting  from  the  transformation  of  the 
stretching-film  problem  to  a  fixed  boundary-value  problem.  Figs.  12  through  17  illustrate 
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results  for  various  combinations  of  Reynolds  number  Re  and  Weber  number  We  at  fixed 
thickness- to-length  ratio  (  and  constant  Tp-value,  i.e.  C  —  0-1  and  Tp  =  1,  respectively. 

Figs.  12  and  13  show  the  variation  of  film  thickness  and  axial  film  velocity  with  time  for 
the  case  with  We  —  1  and  Re  =  100.  Note  that  We  =  1  implies  that  the  maximum  pulling 
velocity  U  equals  the  capillary  wave  velocity  uc  =  JaJ[p5 *).  Precisely,  uc  represents  the 
capillary  wave  speed  for  sinuous  mode  waves  and  not  dilational  waves  as  considered  here; 
the  latter  have  been  shown  to  be  dispersive  with  wave  velocity  inversely  proportional  to 
wavelength.  For  dilational  waves  of  wavelength  greater  than  Si,  the  wave  speed  is  less 
than  ucJ  The  nondimensional  length  =  k/h  at  which  the  film  ruptures  is  2.42,  with 
rupture  occuring  after  nondimensional  time  rb*  =  Tb/(li/U)  =  1.92.  In  this  analysis,  film 
rupture  is  assumed  to  take  place  if  the  value  of  the  film  thickness  5  reaches  zero  or  drops 
below  zero  value  locally. 

From  Fig.  12,  we  observe  the  formation  of  a  liquid  wedge  at  the  pulling  end.  Also, 
capillary  waves  initially  generated  at  the  pulling  end  propagate  into  the  film  leading  to  an 
absolute  maximum  in  the  film  thickness  at  tj  =  x/l  =  0  before  and  at  the  time  of  film  break¬ 
up.  For  the  same  case  as  shown  in  Fig.  12,  Fig.  13  illustrates  the  formation  of  a  singularity 
in  the  axial  velocity  at  the  point  of  film  rupture,  the  latter  being  located  close  to  the  pulling 
end  dividing  the  film  into  a  wedge-shaped  region  and  a  stretched  film  region. 

Fig.  14  shows  the  similar  case  as  in  Fig.  12  but  with  a  decrease  in  Weber  number  from 
We  =  1  to  0.01  (at  fixed  Reynolds  number).  Since  We  =  ( pU26i)/a  =  U2/ul,  where  uc 
denotes  the  capillary  velocity,  the  lower  Weber  number  indicates  that  the  pulling  velocity 
is  significantly  smaller  than  the  velocity  at  which  the  longer  wavelength  dilational  capillary 
waves  propagate  along  the  film.  Accordingly,  the  time  scale  for  the  exchange  of  information 
along  the  liquid  film  is  significantly  smaller  than  the  time  scale  associated  with  the  generation 
of  film  disturbances  due  to  the  accelerating  film.  Consequently,  a  distinct  liquid  wedge  at 
the  pulling  end  is  no  longer  formed.  As  the  film  stretches  over  a  longer  period  of  time 
without  rupture,  the  film  thickness  at  77  =  0  continuously  decreases  at  later  times  and  its 
maximum  value  is  reached  at  the  pulling  end.  The  overall  film  length  at  breakup  is  increased, 
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i.e.  §  =  4.36  with  film  rupture  taking  place  at  a  significantly  later  time  r6*  =  35.75  than  for 
the  case  shown  in  Figs.  12  and  13. 

Fig.  15  displays  numerical  results  obtained  for  the  case  with  Re  =  100  and  We  =  100. 
In  this  case,  the  maximum  pulling  velocity  is  significantly  larger  than  the  capillary  wave 
speed.  Accordingly,  the  effects  of  film  thinning  by  film  stretching  are  not  communicated 
fast  enough  towards  the  film  center  such  that,  at  the  point  of  film  rupture,  the  center  of 
the  film  remains  essentially  undisturbed  with  film  thickness  6*(r]  =  0,  rfc*)  close  to  1.  Due  to 
the  significant  increase  in  stretching/pulling  velocity,  the  film  breaks  very  early  after  only 
r6*  =  0.79.  However,  due  to  the  large  value  of  U  the  film  has  already  stretched  significantly 
at  the  time  of  rupture,  i.e.  II  =  3.98.  Note  that  this  value  is  larger  than  the  one  predicted  for 
the  case  with  decreased  Weber  number,  whereas  the  opposite  is  true  for  the  observed  break¬ 
up  time  Tj.  Fig.  15  illustrates  that  an  increase  in  maximum  pulling  velocity  also  results 
in  a  narrower  film  wedge  (both  in  the  transformed  and  the  original  spatial  coordinates), 
i.e.  thinner  thickness  boundary-layer  at  the  pulling  end,  and  film  breakup  closer  to  the  pulling 
end.  It  is  also  interesting  to  note  that  for  the  case  shown  in  Fig.  15,  a  singular  behavior  in 
the  axial  velocity  at  the  point  of  film  rupture  is  no  longer  observed  (not  illustrated).  Notable 
differences  between  the  results  illustrated  in  Fig.  15  and  the  similar  case  with  long-range 
molecular  forces  considered,  only  appeared  for  film  thicknesses  below  5*  =  0.01. 

Fig.  16  illustrates  a  case  for  small  Weber  number  and  Reynolds  number,  i.e.  We  —  0.01 
and  Re  =  1.  This  can  reflect  a  hundred-fold  decrease  in  velocity  from  the  case  of  Fig.  15. 
Comparison  with  Fig.  14  shows  that  break-up  time  and  film  length  at  rupture  increase  as 
the  Reynolds  number  is  decreased  from  100  to  1  (i.e.  I ;  =  16.8  at  r6*  =  162.64).  The  break¬ 
up  time  and  length  are  also  much  greater  here  than  for  the  case  of  Fig.  15.  A  long  thin 
liquid  filament  or  band  of  “uniform”  thickness  is  formed  with  a  film  wedge  located  at  the 
pulling  end.  Filament  pinching  occurs  at  its  transition  into  the  film  wedge.  Furthermore,  in 
comparison  to  the  case  shown  in  Fig.  12  and  Fig.  13,  the  amplitude  of  the  velocity  singularity 
which  is  forming  at  the  pinch-point  is  greatly  reduced. 

The  influence  of  a  reduction  in  Reynolds  number  from  Re  =  100  to  1  for  a  larger  Weber 
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number  case,  i.e.  We  =  100,  can  be  observed  by  comparing  Fig.  15  and  Fig.  17.  For  the 
larger  Re  number  case  in  Fig.  15,  information  is  propagated  into  the  film  mainly  due  to 
capillary  waves,  however,  since  the  pulling  velocity  in  this  case  is  significantly  larger  than 
the  capillary  wave  velocity,  the  center  of  the  film  at  77  =  0  remained  virtually  undisturbed 
by  the  pulling  of  its  ends.  On  the  other  hand,  for  the  same  value  of  Weber  number  but 
smaller  Reynolds  number,  propagation  of  information  from  the  pulling  end  into  the  film 
is  greatly  determined  by  viscous  diffusion.  Consequently,  the  film  thickness  in  Fig.  17  has 
a  local  maximum  at  the  centerline  whose  value  continuously  decreases  as  time  progresses. 
At  the  same  time,  film  length  at  breakup  and  breakup  time  (i.e.  =  3.77  and  rt*  =  0.77) 

remain  close  to  the  values  predicted  for  the  larger  Reynolds  number  case  of  Fig.  15.  This 
result,  and  the  observation  that  the  thickness  boundary  layer  in  Fig.  17  at  the  pulling  end  is 
slightly  smaller  than  that  displayed  in  Fig.  15,  shows  that  viscous  diffusion  results  in  a  more 
uniform  film  thickness  in  the  stretched  film  region  outside  of  the  thickness  boundary-layer, 
i.e.  away  from  the  film  wedge  at  the  pulling  end.  Also,  whereas  for  the  smaller  Reynolds 
number  case  film  thickness  is  more  uniform  away  from  the  thickness  boundary  layer,  film 
pinching  in  the  vicinity  of  the  liquid  wedge  is  more  pronounced  than  for  the  corresponding 
larger  Reynolds  number  case  at  the  same  time  (e.g.  r*  =  0.7).  In  other  words,  the  local 
minimum  in  film  thickness  at  the  transition  from  film  wedge  to  the  outer  stretched-film 
region  is  more  distinct  in  Fig.  17  than  in  Fig.  15.  In  fact,  for  the  larger  Reynolds  number 
case  film  distortion  might  be  characterized  as  ’film  tearing’  rather  than  ‘film  pinching.’ 

Comparison  of  Fig.  17  with  Fig.  16  also  illustrates  that  for  small  Reynolds  numbers,  an 
increase  in  Weber  number  results  in  a  decrease  in  the  size  of  the  film  thickness  boundary- 
layer.  Note  that  such  a  layer  w-as  not  observed  for  the  case  in  Fig.  14  with  smaller  viscous 
diffusion  in  comparison  to  the  case  of  Fig.  16.  This  indicates,  of  course,  that  viscous  diffusion 
not  only  serves  the  propagation  of  information,  but  it  also  is  effective  in  damping  wave 
amplitudes  (via  viscous  dissipation)  and  in  delaying  film  rupture,  the  latter  occuring  at 
Tj  =  35.75  and  l*b  =  4.36  in  Fig.  14,  whereas  r6*  =  16.80  and  r6*  =  162.64  for  the  case  in 
Fig.  16. 
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Fig.  18  illustrates  a  case  with  Re  —  1  and  We  =  1.  Analogous  to  the  similar  case  with 
larger  Weber  number  shown  in  Fig.  17,  a  distinct  maximum  in  the  film  thickness  develops 
at  the  film  center  with  its  maximum  value  continuously  decreasing.  However,  the  stabilizing 
effect  of  surface  tension  prevents  film  rupture  in  the  lower  Weber  number  case  so  that 
film  stretching  progresses  over  a  much  longer  time-period  resulting  in  more  significant  film 
stretching  before  film  rupture  ( l b  =  8.63,  r6*  =  8.13).  In  this  context,  the  results  of  Fig.  18 
for  We  =  1  resemble  those  of  Fig.  16  for  We  =  0.01  (both  for  Re  =  1),  even  though  in  the 
latter  case  the  film  stretches  over  a  longer  distance  resulting  in  a  nearly  uniform  very  thin 
film  region  near  the  film  center. 

Comparison  of  Figs.  12, 18  and  19  illustrates  the  effect  of  changes  in  Reynolds  number  at 
low  Weber  number  (i.e.  We  =  1).  As  the  Reynolds  number  is  decreased  from  Re  =  100  to 
1,  capillary  waves  are  damped  by  viscous  action,  film  rupture  is  delayed  and  film  stretching 
is  enhanced.  Fig.  19  shows  that  for  a  very  viscous  film  (i.e.  very  small  Re- values),  the  film 
thickness  essentially  decreases  uniformly  over  the  entire  film  as  film  stretching  progresses. 
Film  rupture  eventually  occurs  in  the  immediate  vicinity  of  the  pulling  ends  where  8*  =  1 
is  enforced.  The  prescribed  observation  has  also  been  made  for  the  case  of  a  compressed 
highly  viscous  film  (see  below).  It  remains  to  be  noted  that,  break-up  length  and  time  for 
the  case  shown  in  Fig.  19  (i.e.,  l*b  =  5.52  and  r6*  =  5.22)  are  not  very  different  from  those 
observed  for  the  case  shown  in  Fig.  18. 

Effect  of  Changes  in  (  or  Tp 

The  various  cases  discussed  in  the  previous  section  have  been  analyzed  with  respect  to 
changes  in  thickness-to-length  ratio  of  the  film  (  as  well  as  with  respect  to  changes  in  the 
time  interval  Tp  over  which  the  films  ends  are  accelerated  from  zero  velocity  to  the  maximum 
pull  velocity  U,  in  particular,  C  was  decreased  from  0.1  to  0.01  and  Tp  was  increased  from  1 
to  10,  respectively.  Note  that,  a  smaller  C  value  corresponds  to  a  longer  film  assuming  the 
film  thickness  is  kept  constant. 

Figs.  20  through  24  illustrate  instantaneous  film  thickness  distributions  for  the  various 
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cases  (characterized  by  We  and  Re  values)  near  the  time  of  film  rupture.  Solid  lines  represent 
results  for  the  “base”  cases  discussed  earlier,  dashed  and  dot-dashed  lines  denote  the  results 
for  the  corresponding  cases  (i.e.,  with  same  We-  and  Re- values)  but  with  decreased  £- 
value  or  increased  Tp-value,  respectively.  Fig.  23  contains  an  additional  result  for  the  case  of 
£  =  0.025.  Note  that  the  results  are  plotted  versus  x/l,  so  that  for  direct  comparison  between 
the  £  =  0.1  results  (solid  and  dash-dotted  lines)  and  the  £  =  0.01  result  (dashed  lines),  the 
latter  has  to  be  stretched  by  a  factor  of  ten  to  give  the  dimensional  length  comparison, 
assuming  that  the  film  thickness  is  the  same  in  all  three  cases.  Also,  the  characteristic  time 
used  in  the  nondimensionalization  is  rTej  =  l{/uc.  This  implies  that,  for  direct  comparison 
between  the  £  =  0.01  and  the  £  =  0.1  results  in  one  particular  figure,  the  time  associated 
with  the  £  =  0.01  result  has  to  be  multiplied  by  a  factor  of  ten  to  yield  the  same  dimensional 
time,  again  assuming  the  initial  film  thickness  is  the  same  in  all  3  cases. 

The  figures  illustrate  that,  for  a  particular  combination  of  We  and  Re,  nondimensional 
film  length,  i.e.  l/U  at  the  point  of  rupture  is  significantly  smaller  for  the  £  =  0.01  case. 
That  is,  the  initially  thinner  sheet  requires  a  smaller  fractional  increase  in  length  in  order 
to  break.  However,  dimensional  time  until  break-up  and  dimensional  break-up  length  will 
be  significantly  larger  for  the  lower  Rvalue.  Recall  in  this  context,  that  Trej  =  li/uc.  The 
prescribed  observation  is  particularly  obvious  in  Figs.  20  and  22  through  24.  For  £  =  0.01 
with  Re  =100  and  We  =  1  or  100  (Figs.  20  and  22),  the  center  region  of  the  film  remains 
essentially  undisturbed  throughout  the  simulation.  However,  for  small  Weber  numbers  (i.e. 
We=0.01  at  Re  =  1  or  100),  capillary  waves  from  one  end  of  the  film  are  able  to  reach  the 
film  center  prior  to  film  rupture  resulting  in  a  disturbed  film  center  region.  (See  Figs.  21 
and  23.)  The  similar  observation  is  made  for  smaller  Reynolds  numbers  but  larger  Weber 
numbers  (i.e.  Re=l  at  We  =  100).  Here,  the  viscous  forces  cause  the  film  thickness  to 
decrease  in  the  center  region  even  for  the  smaller  (-value  case. 

If  the  acceleration  of  the  film  ends  is  reduced  by  increasing  Tp  from  1  to  10,  the  films 
stretch  further  before  rupture  occurs.  This  is  particularly  obvious  from  Figs.  20,  22  and  24. 
As  the  films  are  stretched  over  a  longer  time  interval  (prior  to  rupture),  the  film  thickness  in 
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the  film  center  region  will  be  significantly  smaller  at  the  time  of  film  rupture  (in  comparison 
to  the  similar  cases  but  smaller  Tv  value).  See  Figs.  20,  22,  24.  The  prescribed  observation 
is  less  pronounced  in  the  low  Weber  number,  low  Reynolds  number  case,  where  the  film 
stretches  significantly  even  in  the  Tp  =  1, ,  £  =  0.1  “base”  case.  Due  to  the  delay  in  film 
rupture  the  maximum  film  thickness  in  the  film  center  at  the  time  of  rupture  will  be  smaller 
for  larger  Tv  value.  However,  due  to  the  decreased  film  acceleration  in  this  case,  the  film 
thickness  in  the  film  center  region  at  a  given  nondimensional  time  will  be  larger  for  smaller 
Tp  values.  See  Fig.  21  in  this  context. 

Table  1  summarizes  the  various  parameter  configurations  analyzed  within  the  present 
study  including  predictions  of  break-up  length  Tb,  break-up  time  r6*/C  and  film  thickness  in 
the  film  center  at  the  time  of  film  rupture  6*(rb,r]  =  0)  for  the  various  cases.  Note  that, 
Tb  /C  is  proportional  to  the  dimensional  break-up  time  and  5*(rb  ,rj  =  0)  indicates  the  extent 
to  which  the  propagation  of  disturbances  towards  the  center  has  proceeded  at  the  time  of 
film  rupture. 

Film  Compression 

In  order  to  illustrate  how  the  dynamics  of  a  planar  film  bridge  can  relate  to  the  dynamics  of 
a  free  planar  film,  the  present  analysis  was  extended  to  include  negative  film  end-velocities, 
effectively  leading  to  film  compression.  Fig.  25a  illustrates  a  film  compression  case  with 
velocity  boundary  conditions  according  to  Eq.  (30)  but  with  -U/2  instead  of  Uj 2.  The 
Weber  number  and  Reynolds  numbers  in  this  case  were  1  and  100,  respectively.  A  fluid 
blob  (i.e.  a  liquid  cylinder)  forms  at  the  compression  end(s)  with  capillary  waves  proceeding 
the  cylinder  and  propagating  towards  the  film  center.  Figure  25b  shows  the  predictions 
obtained  for  the  similar  case  but  Re  =  1  instead  of  100.  Here,  a  liquid  rim  is  still  formed  at 
the  film  ends.  Due  to  increased  viscous  damping,  however,  no  capillary  waves  proceed  the 
rim.  As  the  Reynolds  number  is  further  decreased,  information  from  the  film  ends  is  rapidly 
propagated  to  the  film  center  via  viscous  forces.  Consequently,  the  film  thickness  increases 
more  or  less  uniformly  across  the  entire  film  except  near  the  film  ends  due  to  the  prescribed 


23 


boundary  condition  S  *  (r*)  =  1. 

Figures  26a-c  illustrate  computational  results  presented  by  Brenner  and  Gueyffier17  for  a 
freely  contracting  film  and  for  Weber  numbers  and  Reynolds  numbers  comparable  to  those 
shown  in  Figs.  25a-c,  respectively.  The  observations  made  in  Ref.  17  for  the  freely  contract¬ 
ing  film  are  essentially  the  same  as  those  described  earlier  for  the  compressed  film  bridge. 

Neglect  of  Intermolecular  Force  Term 

Figure  27  illustrates  results  for  the  same  case  as  in  Fig.  12  but  with  the  effect  on  intermolec¬ 
ular  forces  (modelled  via  the  potential  $)  neglected.  Comparison  between  Figs.  12  and  27 
shows  that  omission  of  the  destabilizing  effect  of  long-range  intermolecular  forces  can  prevent 
film  rupture  and  produce  non-physical  global  flow  features,  such  as  the  continuous  shedding 
of  smaller  ligaments  from  the  pulling  end.  The  ligament  shedding  has  not  been  observed  for 
any  of  the  other  cases  discussed  in  this  work.  However,  it  illustrates  that  the  consideration 
of  a  continuum  model  alone  for  the  analysis  of  film  rupture  can  produce  misleading  results. 


IV.  CONCLUDING  REMARKS 

The  dynamic  stretching  (and  compression)  of  a  thin  viscous  liquid  film  in  a  void  has 
been  analyzed  by  means  of  reduced-dimension  or  lubrication  approach.  The  numerical  im¬ 
plementation  of  the  nonlinear  governing  equations  has  been  validated  for  the  linear  limit  by 
comparison  with  analytical  linear  results  for  a  film  bridge  with  harmonic  modulations  of  the 
axial  film-end  velocities.  Nonlinear  solutions  for  the  continuously  stretching  film  have  been 
presented  for  various  Weber  number  and  Reynolds  number  combinations,  for  two  different 
initial  film-thickness-to-film-length  ratios  and  at  two  different  acceleration  conditions  for  the 
film  ends.  Consistent  results  are  found  from  both  the  linear  and  nonlinear  analyses.  The 
importance  of  the  film  stretching  rate  in  comparison  to  either  the  capillary  wave  velocity 
or  the  characteristic  viscous  velocity  has  been  demonstrated.  If  the  stretching  rate  is  large, 
the  distortion  signals  are  slow  in  reaching  large  distances  from  the  forced  end;  so  distortion 
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is  confined  to  a  smaller  portion  of  the  sheet.  For  small  stretching  rates,  the  distortion  is 
significant  over  a  larger  portion  of  the  film. 

For  films  with  an  initial  length-to-thickness  ratio  of  0(10)  and  with  a  time-scale  for  film- 
end  acceleration  (from  zero  to  maximum  pull  velocity)  comparable  to  the  propagation  time 
of  capillary  waves  from  film  end  through  film  center,  various  distinct  film  topologies  were 
observed,  depending  on  Weber  number  and  Reynolds  number.  For  cases  with  length-to- 
thickness  ratio  of  0(100)  information  on  the  acceleration  of  the  film  ends  does  not  reach  the 
film  center  region.  Here,  evolution  of  film  topology  is  similar  for  all  the  considered  Weber 
number  and  Reynolds  number  combinations;  with  a  more  or  less  undisturbed  film  center 
region  smoothly  transitioning  into  a  narrow  liquid  wedge  at  the  pulling  end.  For  cases  with 
low  film-end  acceleration,  film  rupture  is  significantly  delayed  resulting  (in  all  the  considered 
cases)  in  a  significantly  stretched  film  prior  to  rupture  with  or  without  significant  amounts 
of  fluid  remaining  in  the  film  center  region. 

Nonlinear  film  distortion  has  also  been  investigated  for  the  continuously  compressed 
planar  film  bridge,  illustrating  the  relevance  of  the  film  bridge  analysis  for  the  contracting 
and  stretching  of  a"  free  planar  film  such  as  those  found  in  the  atomization  process  of  liquid 
fuels  in  typical  gas  turbine  combustors. 
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Tables 


We 

Re 

c 

EHI 

i; 

Tb  K 

S"(n,n  =  o) 

Figure  No. 

1 

100 

0.1 

mm 

2.42 

19.2 

1.25 

12,  20  (solid) 

1 

100 

0.01 

1 

1.52 

102 

1 

20  (dashed) 

1 

100 

0.1 

10 

6.45 

104.4 

0.44 

20  (dash-dotted) 

100 

100 

0.1 

1 

3.98 

7.9 

0.98 

15,  22  (solid) 

100 

100 

0.01 

1 

2.04 

53 

1 

22  (dashed) 

100 

100 

0.1 

EH 

12.32 

54.2 

0.68 

22  (dash-dotted) 

0.01 

1 

0.1 

1 

16.80 

1626.4 

0.01 

16,  23  (solid) 

0.01 

1 

0.025 

mm 

6.14 

2074.4 

0.35 

23  (dot-dot-dashed) 

0.01 

1 

0.01 

1 

3.03 

1.14 

23  (dashed) 

0.01 

1 

0.1 

EH 

16.31 

0.02 

23  (dash-dotted) 

0.01 

100 

1 

4.36  ! 

357.5 

0.03 

14,  21  (solid) 

0.01 

100 

0.01 

mm 

1.50 

552 

1.07 

21  (dashed) 

0.01 

100 

0.1 

EH 

4.11 

373.8 

21  (dash-dotted) 

100 

1 

0.1 

mm 

3.77 

7.7 

0.49 

17,  24  (solid) 

100 

1 

0.01 

1 

2.08 

53 

0.96 

24  (dashed) 

pliM 

1 

EH 

7.65 

44.7 

0.22 

24  (dash-dotted) 

■ 

1 

1 

8.63 

0.29 

18 

1 

0.01 

— 

5.52 

52.2 

0.17 

19 

1 

100 

0.01 

m 

0.83 

176 

1 

26  (a) 

1 

1 

0.01 

5 

0.70 

215 

1 

26  (b) 

1 

0.01 

0.01 

5 

0.56 

247 

1.65 

26  (c) 

Table  1:  Summary  of  parameter  survey  including  results  for  nondimensional  values  of  break-up  length  / j , 
break-up  time  r6*/C  and  film  thicknes  at  the  film  center  during  film  rupture  =  0). 
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Figure  captions 


Figure  1:  Schematic  of  investigated  free  stretching  film  configuration  (h  —  5). 

Figure  2:  Solid:  Linear  analytical  solution  for  Uq  =  0.01,  Re  =  100,/  =  0.1,  L  =  10.  Dashed. 
Numerical  solution  for  the  same  case. 

Figure  3:  Solid:  Linear  analytical  solution  for  Uq  =  0.01,  Re  =  10,/  =  0.1,  L  —  10.  Dashed: 
Numerical  solution  for  the  same  case. 

Figure  4:  Solid:  Linear  analytical  solution  for  U0  =  0.01,  Re  =  1,/  =  0.1,  L  =  10.  Dashed: 
Numerical  solution  for  the  same  case. 

Figure  5:  Linear  analytical  solution  for  Uo  —  0.01,  Re  =  0.2,  /  =  0.1,  L  =  10. 

Figure  6:  Linear  analytical  solution  for  Uo  =  0.01,  Re  —  —  0.1,  L  =  100. 

Figure  7:  Linear  analytical  solution  for  Uq  —  0.01,  Re  =  l,f  =  0.1,  L  =  100. 

Figure  8:  Linear  analytical  solution  for  Uo  =  0.01,  Re  =  0.2,  /  =  0.1,  L  =  100. 

Figure  9:  Linear  analytical  solution  for  Uq  =  0.01,  Re  =  1,  /  =  0.1,  L  =  250. 

Figure  10:  Linear  analytical  solution  for  Uo  —  0.01,  Re  =  1,  /  =  0.5,  L  =  10. 

Figure  11:  Linear  analytical  solution  for  Uo  =  0.01,  Re  =  0.2,  /  =  0.5,  L  =  10. 

Figure  12:  Film  thickness  5*  as  a  function  of  time  r*  and  normalized  spatial  coordinate  77  for 
We  =  1  and  Re  =  100  [l*b  =  2.42,  r6*  =  1.92], 

Figure  13:  Axial  film  velocity  u*  as  a  function  of  time  r*  and  normalized  spatial  coordinate  77  for 
We  =  1  and  Re  =  100  [i;  =  2.42,  r*b  =  1.92]. 

Figure  14:  Film  thickness  5*  as  a  function  of  time  r*  and  normalized  spatial  coordinate  77  for 
We  =  0.01  and  Re  =  100  [ZJ  =  4.36,  r6*  =  35.75]. 


29 


Figure  15:  Film  thickness  8*  as  a  function  of  time  r*  and  normalized  spatial  coordinate  77  for 
We  =  100  and  Re  =  100  [Z6*  =  3.98,  =  0.79]. 

Figure  16:  Film  thickness  5*  as  a  function  of  time  r*  and  normalized  spatial  coordinate  77  for 
We  =  0.01  and  Re  =  1  [i;  =  16.80,  r4*  =  162.64], 

Figure  17:  Film  thickness  8*  as  a  function  of  time  r*  and  normalized  spatial  coordinate  rj  for 
We  =  100  and  Re  =  1  [l*b  =  3.77,  rft*  =  0.77]. 

Figure  18:  Film  thickness  8*  as  a  function  of  time  r*  and  normalized  spatial  coordinate  77  for 
We  =  1  and  Re  =  1  [l*b  =  8.63,  rb*  =  8.13].  . 

Figure  19:  Film  thickness  8*  as  a  function  of  time  r*  and  normalized  spatial  coordinate  77  for 
We  =  1  and  Re  =  0.01  [l*b  =  5.52,  r6*  =  5.22], 

Figure  20:  Film  thickness  8*  at  various  times  r*  plotted  versus  normalized  spatial  coordinate  x/k 
for  We  =  1  and  Re  =  100.  Solid:  <^=0.1,  Tp  =  1,  dashed:  •Ji//m-=0.01,  Tp  =  1, 
dash-dotted:  <5i/Zi=0.1,  Tp  =  10. 

Figure  21:  Film  thickness  <5*  at  various  times  r*  plotted  versus  normalized  spatial  coordinate  x/U 
for  We  =  0.01  and  Re  =  100.  Solid:  8^= 0.1,  Tp  =  1,  dashed:  <^=0.01,  Tp  =  1, 
dash-dotted:  8j/li= 0.1,  Tp  =  10, 

Figure  22:  Film  thickness  8*  at  various  times  r*  plotted  versus  normalized  spatial  coordinate  x/li 
for  We  =  100  and  Re  =  100.  Solid:  8i/U= 0.1,  Tp  =  1,  dashed:  5i/Zi=0.01,  Tp  =  1, 
dash-dotted:  <5,/Z,=0.1,  Tp  =  10. 

Figure  23:  Film  thickness  8 *  at  various  times  r*  plotted  versus  normalized  spatial  coordinate  x/lt 
for  We  =  0.01  and  Re  =  1.  Solid:  £,//, =0.1,  Tp  =  1,  dashed:  8t/lt=  0.01,  Tp  =  1, 
dash-dotted:  8i/k=  0.1,  Tp  =  10,  dash-dot-dotted:  5,//, =0.025,  Tp  =  1. 

Figure  24:  Film  thickness  8*  at  various  times  r*  plotted  versus  normalized  spatial  coordinate  x/li 
for  We  =  100  and  Re  =  1.  Solid:  £//<={).  1,  Tp  =  1,  dashed:  d,//,=0.01.  Tp  =  1, 
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dash-dotted:  <5i/Zi=0.1,  Tv  =  10, 


Figure  25:  Film  compression  for  We  =  1  and  Re  =  100  (a),  1  (b)  and  0.01  (c). 

Figure  26:  Film  compression  for  We  =  1,  Re  =  1.  Results  presented  in  Ref.  17  for  contracting  free 
planar  liquid  film  and  Weber  and  Reynolds  numbers  comparable  to  those  presented  in 
Fig.  25. 

Figure  27:  Film  thickness  5*  as  a  function  of  time  r*  and  normalized  spatial  coordinate  77  for  We  = 
1  and  Re  =  100  and  without  the  consideration  of  intermolecular  forces,  i.e.  A*  =  0. 
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FIGURES 


x=0:  dh/dx=0  ,  u=0  x=l(t)=uP(t)  t:  h=ho,u=uP(t) 


FIG.  1.  Mehring,  Phys.  Fluids 


FIG.  2.  Mehring,  Phys.  Fluids 
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FIG.  10.  Mehring,  Phys.  Fluids 
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FIG.  18.  Mehring,  Phys.  Fluids 


■n 


FIG.  19.  Mehring,  Phys.  Fluids 
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FIG.  21.  Mehring,  Phys.  Fluids 
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FIG.  22.  Mehring,  Phys.  Fluids 
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FIG.  24.  Mehring,  Phys.  Fluids 
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FIG.  26.  Mehring,  Phys.  Fluids 
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FIG.  27.  Mehring,  Phys.  Fluids 
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